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Abstract 

The first LHC results at 7-8 TeV, with the discovery of a candidate Higgs boson 
and the non observation of new particles or exotic phenomena, have made a big step 
towards completing the experimental confirmation of the Standard Model (SM) of 
fundamental particle interactions. It is thus a good moment for me to collect, update 
and improve my graduate lecture notes on Quantum Chromodynamics (QCD) and 
the theory of Electroweak (EW) Interactions, with main focus on Collider Physics. I 
hope that these lectures can provide an introduction to the subject for the interested 
reader, assumed to be already familiar with quantum field theory and some basic 
facts in elementary particle physics as taught in undergraduate courses. 
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1 Gauge Theories and the Standard Model 



1.1 An Overview of the Fundamental Interactions 

A possible goal of fundamental physics is to reduce all natural phenomena to a set of 
basic laws and theories that, at least in principle, can quantitatively reproduce and predict 
the experimental observations. At microscopic level all the phenomenology of matter and 
radiation, including molecular, atomic, nuclear and subnuclear physics, can be understood 
in terms of three classes of fundamental interactions: strong, electromagnetic and weak 
interactions. For all material bodies on the Earth and in all geological, astrophysical and 
cosmological phenomena a fourth interaction, the gravitational force, plays a dominant 
role, while it is instead negligible in atomic and nuclear physics. In atoms the electrons are 
bound to nuclei by electromagnetic forces and the properties of electron clouds explain the 
complex phenomenology of atoms and molecules. Light is a particular vibration of electric 
and magnetic fields (an electromagnetic wave). Strong interactions bind the protons and 
neutrons together in nuclei, being so intensively attractive at short distances that they 
prevail over the electric repulsion due to the equal sign charges of protons. Protons and 
neutrons, in turn, are composites of three quarks held together by strong interactions to 
which quarks and gluons are subject (hence these particles are called "hadrons" from the 
Greek word for "strong"). To the weak interactions are due the beta radioactivity that 
makes some nuclei unstable as well as the nuclear reactions that produce the enormous 
energy radiated by the stars and by our Sun in particular. The weak interactions also 
cause the disintegration of the neutron, the charged pions, the lightest hadronic particles 
with strangeness, charm, and beauty (which are "flavour" quantum numbers) as well as 
the decay of the quark top and of the heavy charged leptons (the muon fj,~ and the tau 
t~). In addition all observed neutrino interactions are due to these weak forces. 

All these interactions (with the possible exception of gravity) are described within the 
framework of quantum mechanics and relativity, more precisely by a local relativistic quan- 
tum field theory. To each particle, described as pointlike, is associated a field with suitable 
(depending on the particle spin) transformation properties under the Lorentz group (the 
relativistic space-time coordinate transformations). It is remarkable that the description 
of all these particle interactions is based on a common principle: "gauge" invariance. A 
"gauge" symmetry is invariance under transformations that rotate the basic internal de- 
grees of freedom but with rotation angles that depend on the space-time point. At the 
classical level gauge invariance is a property of the Maxwell equations of electrodynamics 
and it is in this context that the notion and the name of gauge invariance were intro- 
duced. The prototype of all quantum gauge field theories, with a single gauged charge, is 
QED, Quantum Electro- Dynamics, developed in the years from 1926 until about 1950, 
which indeed is the quantum version of Maxwell theory. Theories with gauge symmetry 
in 4 space-time dimensions are renormalizable and are completely determined given the 
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symmetry group and the representations of the interacting fields. The whole set of strong, 
electromagnetic and weak interactions is described by a gauge theory, with 12 gauged non- 
commuting charges, which is called "the Standard Model" of particle interactions (SM). 
Actually only a subgroup of the SM symmetry is directly reflected in the spectrum of 
physical states. A part of the electroweak symmetry is hidden by the Higgs mechanism for 
spontaneous symmetry breaking of the gauge symmetry. 

The theory of general relativity is a classic description (in the sense of non quantum 
mechanical) of gravity that goes beyond the static approximation described by Newton law 
and includes dynamical phenomena like, for example, gravitational waves. The problem 
of formulating a quantum theory of gravitational interactions is one of the central prob- 
lems of contemporary theoretical physics. But quantum effects in gravity become only 
important for energy concentrations in space-time which are not in practice accessible to 
experimentation in the laboratory. Thus the search for the correct theory can only be done 
by a purely speculative approach. All attempts at a description of quantum gravity in 
terms of a well defined and computable local field theory along similar lines as for the SM 
have so far failed to lead to a satisfactory framework. Rather, at present the most com- 
plete and plausible description of quantum gravity is a theory formulated in terms of non 
pointlike basic objects, the so called "strings", extended over distances much shorter than 
those experimentally accessible, that live in a space-time with 10 or 11 dimensions. The 
additional dimensions beyond the familiar 4 are, typically, compactified which means that 
they are curled up with a curvature radius of the order of the string dimensions. Present 
string theory is an all-comprehensive framework that suggests a unified description of all 
interactions together with gravity of which the SM would be only a low energy or large 
distance approximation. 

A fundamental principle of quantum mechanics, the Heisenberg indetermination princi- 
ple, implies that, for studying particles with spatial dimensions of order Ax or interactions 
taking place at distances of order Ax, one needs as a probe a beam of particles (typi- 
cally produced by an accelerator) with impulse p > h/Ax, where h is the reduced Planck 
constant (h = hj1-n\ Accelerators presently in operation, like the Large Hadron Collider 
(LHC) at CERN near Geneva, allow to study collisions between two particles with total 
center of mass energy up to 2E ~ 2pc < 7 — 14 TeV. These machines, in principle, can 
allow to study physics down to distances Ax > 10~ 18 cm. Thus, on the basis of results 
from experiments at existing accelerators, we can confirm that, down to distances of that 
order of magnitude, indeed electrons, quarks and all the fundamental SM particles do not 
show an appreciable internal structure and look elementary and pointlike. We expect that 
quantum effects in gravity will certainly become important at distances Ax < 10~ 33 cm 
corresponding to energies up to E ~ Mp/c 2 ~ 10 19 GeV, where Mpi is the Planck mass, 
related to Newton constant by Gjy = hc/M^. At such short distances the particles that 
so far appeared as pointlike could well reveal an extended structure, like for strings, and 
be described by a more detailed theoretical framework of which the local quantum field 
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theory description of the SM would be just a low energy/large distance limit. 

From the first few moments of the Universe, after the Big Bang, the temperature 
of the cosmic background went down gradually, starting from kT ~ Mpic 2 , where k = 
8.617...10 -5 eV °K~ 1 is the Boltzmann constant, down to the present situation where 
T ~ 2.725°K. Then all stages of high energy physics from string theory, which is a purely 
speculative framework, down to the SM phenomenology, which is directly accessible to 
experiment and well tested, are essential for the reconstruction of the evolution of the 
Universe starting from the Big Bang. This is the basis for the ever increasing relation 
between high energy physics and cosmology. 

1.2 The Architecture of the Standard Model 

The SM is a gauge field theory based on the symmetry group SU (3) ® SU (2) U(l). 
The transformations of the group act on the basic fields. This group has 8+3+1= 12 gen- 
erators with a non trivial commutator algebra (if all generators commute the gauge theory 
is said to be "abelian", while the SM is a "non abelian" gauge theory). SU(2) §§U(1) de- 
scribes the electroweak (EW) interactions [1]- [3] and the electric charge Q, the generator 
of the QED gauge group U (1)q, is the sum of T 3 , one of the SU{2) generators and of Y/2, 
where Y is the U(l) generator: Q = T 3 + Y/2. SU(3) is the "colour" group of the theory 
of strong interactions (QCD: Quantum Chromo- Dynamics [4]- [6]). 

In a gauge theory 1 to each generator T is associated a vector boson (also said gauge 
boson) with the same quantum numbers as T, and, if the gauge symmetry is unbroken, 
this boson is of vanishing mass. These vector (i.e. of spin 1) bosons act as mediators 
of the corresponding interactions. For example, in QED the vector boson associated to 
the generator Q is the photon 7. The interaction between two charged particles in QED, 
for example two electrons, is mediated by the exchange of one (or seldom more than one) 
photon emitted by one electron and reabsorbed by the other one. Similarly in the SM there 
are 8 gluons associated to the £77(3) colour generators, while for SU{2) (g) £7(1) there are 4 
gauge bosons W + , W~ , Z° and 7. Of these, only the gluons and the photon 7 are massless 
because the symmetry induced by the other 3 generators is actually spontaneously broken. 
The masses of W + , W~ and Z° are quite large indeed on the scale of elementary particles: 
mw ~ 80.4 GeV, rriz ~ 91.2 GeV are as heavy as atoms of intermediate size like rubidium 
and molibdenum, respectively. In the electroweak theory the breaking of the symmetry is 
of a particular type, denoted as spontaneous symmetry breaking. In this case charges and 
currents are as dictated by the symmetry but the fundamental state of minimum energy, 
the vacuum, is not unique and there is a continuum of degenerate states that all together 
respect the symmetry (in the sense that the whole vacuum orbit is spanned by applying 
the symmetry transformations). The symmetry breaking is due to the fact that the system 
(with infinite volume and infinite number of degrees of freedom) is found in one particular 

1 Much of the material in this Chapter is a revision and update of ref. [7] 
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vacuum state, and this choice, which for the SM occurred in the first instants of the Universe 
life, makes the symmetry violated in the spectrum of states. In a gauge theory like the SM 
the spontaneous symmetry breaking is realized by the Higgs mechanism [8-11] (described 
in detail in Sect. (1.7)): there are a number of scalar (i.e. of zero spin) Higgs bosons with a 
potential that produces an orbit of degenerate vacuum states. One or more of these scalar 
Higgs particles must necessarily be present in the spectrum of physical states with masses 
very close to the range so far explored. It now seems that that the Higgs particle has been 
found at the LHC with uih ~ 126 GeV [12, 13] thus making a big step towards completing 
the experimental verification of the SM. 

The fermionic (all of spin 1/2) matter fields of the SM are quarks and leptons. Each 
type of quark is a colour triplet (i.e. each quark flavour comes in three colours) and also 
carries electroweak charges, in particular electric charges +2/3 for up-type quarks and -1/3 
for down-type quarks. So quarks are subject to all SM interactions. Leptons are colourless 
and thus do not interact strongly (they are not hadrons) but have electroweak charges, in 
particular electric charges -1 for charged leptons (e~, and r~) while it is for neutrinos 
(v e , and v r ). Quarks and leptons are grouped in 3 "families" or "generations" with 
equal quantum numbers but different masses. At present we do not have an explanation 
for this triple repetition of fermion families: 

u u u v e 
d d d e 

The QCD sector of the SM (see Chapter 2) has a simple structure but a very rich dynam- 
ical content, including the observed complex spectroscopy with a large number of hadrons. 
The most prominent properties of QCD are asymptotic freedom and confinement. In field 
theory the effective coupling of a given interaction vertex is modified by the interaction. 
As a result, the measured intensity of the force depends on the tranferred (four)momentum 
squared, Q 2 , among the participants. In QCD the relevant coupling parameter that ap- 
pears in physical processes is a s = e 2 /4n where e s is the coupling constant of the basic 
interaction vertices of quark and gluons: qqg or ggg (see eqs.(28,29,30,31)). Asymptotic 
freedom means that the effective coupling becomes a function of Q 2 : a s (Q 2 ) decreases 
for increasing Q 2 and vanishes asymptotically. Thus, the QCD interaction becomes very 
weak in processes with large Q 2 , called hard processes or deep inelastic processes (i.e. with 
a final state distribution of momenta and a particle content very different than those in 
the initial state). One can prove that in 4 space-time dimensions all pure-gauge theories 
based on a non commuting group of symmetry are asymptotically free and conversely. The 
effective coupling decreases very slowly at large momenta with the inverse logarithm of 
Q 2 : a s (Q 2 ) = l/61ogQ 2 /A 2 where b is a known constant and A is an energy of order a 
few hundred MeV. Since in quantum mechanics large momenta imply short wavelenghts, 
the result is that at short distances (or Q > A) the potential between two colour charges 
is similar to the Coulomb potential, i.e. proportional to a s (r)/r, with an effective colour 
charge which is small at short distances. On the contrary the interaction strenght becomes 
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large at large distances or small transferred momenta, of order Q < A. In fact all observed 
hadrons are tightly bound composite states of quarks (baryons are made of qqq and mesons 
of qq), with compensating colour charges so that they are overall neutral in colour. In fact, 
the property of confinement is the impossibility of separating colour charges, like individual 
quarks and gluons or any other coloured state. This is because in QCD the interaction 
potential between colour charges increases at long distances linearly in r. When we try 
to separate the quark and the antiquark that form a colour neutral meson the interaction 
energy grows until pairs of quarks and antiquarks are created from the vacuum and new 
neutral mesons are coalesced and observed in the final state instead of free quarks. For 
example, consider the process e + e~ — > qq at large center of mass energies. The final state 
quark and antiquark have large energies, so they separate in opposite directions very fast. 
But the colour confinement forces create new pairs in between them. What is observed 
is two back-to-back jets of colourless hadrons with a number of slow pions that make the 
exact separation of the two jets impossible. In some cases a third well separated jet of 
hadrons is also observed: these events correspond to the radiation of an energetic gluon 
from the parent quark- antiquark pair. 

In the EW sector the SM (see Chapter 3) inherits the phenomenological successes of 
the old (V — A) <g> (V — A) four-fermion low-energy description of weak interactions, and 
provides a well-defined and consistent theoretical framework including weak interactions 
and quantum electrodynamics in a unified picture. The weak interactions derive their name 
from their intensity. At low energy the strength of the effective four-fermion interaction 
of charged currents is determined by the Fermi coupling constant Gf- For example, the 
effective interaction for muon decay is given by 

£ cff = (G F /V2) [iyy a (l - 75 )/i] [ei a (l - 7sK] , (2) 

with [14] 

Gf = 1.1663787(6) x 1(T 5 GeV" 2 . (3) 

In natural units h — c — 1, Gp has dimensions of (mass) -2 . As a result, the intensity of 
weak interactions at low energy is characterized by GpE 2 , where E is the energy scale for 
a given process (E ~ m M for muon decay). Since 

G F E 2 = G F m 2 p (E/m p ) 2 ~ 10 -5 (£/m p ) 2 , (4) 

where m p is the proton mass, the weak interactions are indeed weak at low energies (up to 
energies of order a few ten's of GeV). Effective four fermion couplings for neutral current 
interactions have comparable intensity and energy behaviour. The quadratic increase with 
energy cannot continue for ever, because it would lead to a violation of unitarity. In fact, 
at large energies the propagator effects can no longer be neglected, and the current-current 
interaction is resolved into current-^ gauge boson vertices connected by a W propagator. 
The strength of the weak interactions at high energies is then measured by gw, the W — fi- 
coupling, or, even better, by aw = Qw/^ analogous to the fine-structure constant a 
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of QED (in Chapter 3, gw is simply denoted by g or g 2 ). In the standard EW theory we 
have 

a w = V2G F m 2 w /n = 1/30 . (5) 

That is, at high energies the weak interactions are no longer so weak. 

The range ry/ of weak interactions is very short: it is only with the experimental 
discovery of the W and Z gauge bosons that it could be demonstrated that r w is non- 
vanishing. Now we know that 

h _, fi 
r w = ~ 2.5 x 10 cm , 

m w c 

corresponding to m-w — 80.4 GeV. This very large value for the W (or the Z) mass makes 
a drastic difference, compared with the massless photon and the infinite range of the QED 
force. The direct experimental limit on the photon mass is [14] m 7 < 1 10~ 18 eV. Thus, 
on the one hand, there is very good evidence that the photon is massless. On the other 
hand, the weak bosons are very heavy. A unified theory of EW interactions has to face 
this striking difference. 

Another apparent obstacle in the way of EW unification is the chiral structure of weak 
interactions: in the massless limit for fermions, only left-handed quarks and leptons (and 
right-handed antiquarks and antileptons) are coupled to Ws. This clearly implies parity 
and charge-conjugation violation in weak interactions. 

The universality of weak interactions and the algebraic properties of the electromagnetic 
and weak currents [the conservation of vector currents (CVC), the partial conservation of 
axial currents (PCAC), the algebra of currents, etc.] have been crucial in pointing to a 
symmetric role of electromagnetism and weak interactions at a more fundamental level. 
The old Cabibbo universality [15] for the weak charged current: 

jweak _ ^ 7a (i _ 75 )^ _|_ j? e7a (l - 75 ) e + COS 9 C U"f a (l - 7 5 )rf + 

+ sin6» c uj a (l - 7 5 )s + ... , (7) 

suitably extended, is naturally implied by the standard EW theory. In this theory the 
weak gauge bosons couple to all particles with couplings that are proportional to their 
weak charges, in the same way as the photon couples to all particles in proportion to their 
electric charges [in Eq. (7), d! = cos8 c d + sin0 c s is the weak-isospin partner of u in a 
doublet. The (u,d') doublet has the same couplings as the (v e ,£) and doublets]. 

Another crucial feature is that the charged weak interactions are the only known in- 
teractions that can change flavour: charged leptons into neutrinos or up-type quarks into 
down-type quarks. On the contrary, there are no flavour-changing neutral currents at tree 
level. This is a remarkable property of the weak neutral current, which is explained by the 
introduction of the Glashow-Iliopoulos-Maiani (GIM) mechanism [16] and has led to the 
successful prediction of charm. 
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The natural suppression of flavour-changing neutral currents, the separate conservation 
of e, \i and r leptonic flavours that is only broken by the small neutrino masses, the mech- 
anism of CP violation through the phase in the quark- mixing matrix [17], are all crucial 
features of the SM. Many examples of new physics tend to break the selection rules of the 
standard theory. Thus the experimental study of rare flavour-changing transitions is an 
important window on possible new physics. 

The SM is a renormalizable field theory which means that the ultra-violet divergences 
that appear in loop diagrams can be eliminated by a suitable redefinition of the parameters 
already appearing in the bare lagrangian: masses, couplings and field normalizations. As 
it will be discussed later, a necessary condition for a theory to be renormalizable is that 
only operator vertices of dimension not larger than 4 (that is m 4 where m is some mass 
scale) appear in the lagrangian density C (itself of dimension 4, because the action S is 
given by the integral of C over d A x and is dimensionless in natural units: H — c — 1). Once 
this condition is added to the specification of a gauge group and of the matter field content 
the gauge theory lagrangian density is completely specified. We shall see the precise rules 
to write down the lagrangian of a gauge theory in the next Section. 

1.3 The Formalism of Gauge Theories 

In this Section we summarize the definition and the structure of a gauge Yang-Mills 
theory [18]. We will list here the general rules for constructing such a theory. Then these 
results will be applied to the SM. 

Consider a lagrangian density C[(j>, <9 M 0] which is invariant under a D dimensional con- 
tinuous group T of transformations: 

<j ) '(x) = U(6 A )<j ) (x) (A =1,2, ...,£>) . (8) 

with: 

U(9 A ) = exp [ig J2 A T A ] ~ l + igJ2 eATA + ■ ■ ■ . ( 9 ) 

A A 

The quantities 9 A are numerical parameters, like angles in the particular case of a rotation 
group in some internal space. The approximate expression on the right is valid for 6 A 
infinitesimal. Then, g is the coupling constant and T A are the generators of the group 
T of transformations (8) in the (in general reducible) representation of the fields 0. Here 
we restrict ourselves to the case of internal symmetries, so that T A are matrices that are 
independent of the space-time coordinates and the arguments of the fields <fi and 0' in 
Eq.(8) are the same. If U is unitary, then the generators T A are hermitian, but this need 
not be the case in general (though it is true for the SM). Similarly if U is a group of 
matrices with unit determinant, then the traces of T A vanish: tr(T ) = 0. In general, the 
generators satisfy the commutation relations 

[T A ,T B ]=iC ABC T c . (10) 
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For A, B, C... up or down indices make no difference: T A = T4 etc. The structure constants 
Cabc are completely antisymmetric in their indices, as can be easily seen. Recall that if 
all generators commute the gauge theory is said to be "abelian" (in this case all the 
structure constants Cabc vanish), while the SM is a "non abelian" gauge theory. We 
choose to normalize the generators T A in such a way that for the lowest dimensional non- 
trivial representation of the group T (we use t A to denote the generators in this particular 
representation) we have 

tr(t A t B ) = l -5 AB . (11) 

A normalization convention is needed to fix the normalization of the coupling g and of the 
structure constants Cabc- In the following, for each quantity f A we define 

f = TA f A ■ (12) 

A 

For example, we can rewrite eq. (9) in the form: 

U(9 A ) = exp [igO] ~ 1 + igO + ..., (13) 

If we now make the parameters 9 A depend on the space-time coordinates 6 A = A (x IJ ), 
£[0,(9^0] is in general no longer invariant under the gauge transformations U[9 A (x fl )], 
because of the derivative terms: indeed d^cj)' = d^(U(f>) 7^ Ud^cp. Gauge invariance is 
recovered if the ordinary derivative is replaced by the covariant derivative: 

D^ = d, + igV^ , (14) 

where V A are a set of D gauge vector fields (in one-to-one correspondence with the group 
generators) with the transformation law 

V; = UV.U- 1 - {l/ig^d^U- 1 . (15) 

For constant 9 A , V reduces to a tensor of the adjoint (or regular) representation of the 
group: 

V; = UV.U- 1 ~ + ig[0, VJ . . . , (16) 

which implies that 

V'f = VZ-gC ABC 9 A V* (17) 

where repeated indices are summed up. 

As a consequence of Eqs. (14) and (15), D^ifi has the same transformation properties 
as 0: 

(D^)' = U(D^) . (18) 

In fact 

(D^)' = (d. + tgV'^' = (d,U)<P + Ud,<P + igUV,<f ) -(d,U)<P = U(D,<p) . (19) 
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Thus £[0, D^cj)] is indeed invariant under gauge transformations. But, at this stage, 
the gauge fields appear as external fields that do not propagate. In order to construct 
a gauge-invariant kinetic energy term for the gauge fields V^ 4 , we consider 

[Dp, D v ]<t> = ig{d^ u - d v V^ + ig[V M V v }}<f> = igF^ , (20) 

which is equivalent to 

K = W - d »K - 9CabcV*VC . (21) 

From Eqs. (8), (18) and (20) it follows that the transformation properties of F^ u are those 
of a tensor of the adjoint representation 

= UF^U- 1 ■ (22) 

The complete Yang-Mills lagrangian, which is invariant under gauge transformations, can 
be written in the form 

£ym = - 1 -TtF, v F^ + £[0, Drf] = - \ £ F*F A »» + £[0, Drf] . (23) 

A A 

Note that the kinetic energy term is an operator of dimension 4. Thus if £ is renormalizable, 
also £ym is renormalizable. If we give up renormalizability then more gauge invariant 
higher dimension terms could be added. It is already clear at this stage that no mass term 
for gauge bosons of the form m 2 V^V tl is allowed by gauge invariance. 



1.4 Application to QED and QCD 

For an abelian theory, as for example QED, the gauge transformation reduces to 
U[6(x)] = exp[ieQ6(x)], where Q is the charge generator (for more commuting genera- 
tors one simply has a product of similar factors). The associated gauge field (the photon), 
according to Eq. (15), transforms as 

K = v * ~ d A x ) ■ ( 24 ) 

and the familiar gauge transformation by addition of a 4-gradient of a scalar function is 
recovered. The QED lagrangian density is given by: 

£ = ~\f^F, u + ^(i#-m^ . (25) 

Here p = -D M 7 M , where 7 M are the Dirac matrices and the covariant derivative is given in 
terms of the photon field and the charge operator Q by: 

D^ = d^ + xeA^Q (26) 
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and 

= d^A u -d u A^ (27) 

Note that in QED one usually takes the e~ to be the particle, so that Q — — 1 and the 
covariant derivative is = — ieA^ when acting on the electron field. In the abelian 
case, the tensor is linear in the gauge field so that in the absence of matter fields 
the theory is free. On the other hand, in the non abelian case the F A V tensor contains both 
linear and quadratic terms in V A , so that the theory is non-trivial even in the absence of 
matter fields. 

According to the formalism of the previous section, the statement that QCD is a renor- 
malizable gauge theory based on the group £77(3) with colour triplet quark matter fields 
fixes the QCD lagrangian density to be 

8 nj 
A=l j=l 

Here q 3 - are the quark fields (of n/ different flavours) with mass m 3 - and is the covariant 
derivative: 

D„ = 0„ + ie s ^ (29) 
e s is the gauge coupling and later we will mostly use, in analogy with QED 

a, = ±. (30) 

Also, g M = J2a t A 9^ where g A , A = 1, 8, are the gluon fields and t A are the SU(3) group 
generators in the triplet representation of quarks (i.e. tA are 3x3 matrices acting on q); 
the generators obey the commutation relations = %Cabc^ C where Cabc are the 

complete antisymmetric structure constants of £77(3) (the normalization of Cabc an d of 
e s is specified by that of the generators t A \ Tr[t A t B ] = 5 AB /2, see eq.(ll)). Finally we 
have: 

K = - d v gt - e s C AB cg*g C u (31) 

Chapter 2 is devoted to a detailed description of the QCD as the theory of strong inter- 
actions. The physical vertices in QCD include the gluon-quark-antiquark vertex, analogous 
to the QED photon-fermion-antifermion coupling, but also the 3-gluon and 4-gluon ver- 
tices, of order e s and e 2 s respectively, which have no analogue in an abelian theory like QED. 
In QED the photon is coupled to all electrically charged particles but itself is neutral. In 
QCD the gluons are coloured hence self-coupled. This is reflected in the fact that in QED 
F^y is linear in the gauge field, so that the term F? v in the lagrangian is a pure kinetic 
term, while in QCD F A is quadratic in the gauge field so that in F A2 we find cubic and 
quartic vertices beyond the kinetic term. Also instructive is to consider the case of scalar 
QED: 

C = -\f^F^ + (7J M 0)t( J D^)-m 2 (0V) (32) 
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For Q = 1 we have: 

(D^(D^) = (d^y(dy) + ieA^d^cj) - tf(dy)] + e 2 A^^<P (33) 

We see that for a charged boson in QED, given that the kinetic term for bosons is quadratic 
in the derivative, there is a gauge-gauge-scalar-scalar vertex of order e 2 . We understand 
that in QCD the 3-gluon vertex is there because the gluon is coloured and the 4-gluon 
vertex because the gluon is a boson. 

1.5 Chirality 

We recall here the notion of chirality and related issues which is crucial for the formu- 
lation of the EW Theory. The fermion fields can be described through their Right Handed 
(RH) (chirality +1) and Left Handed (LH) (chirality -1) components: 

*Pl,r = [(1 T 7 5 )/2]^, fej = $[(1± 7s)/2] , (34) 

with 75 and other Dirac matrices defined as in the book by Bjorken-Drell [19] . In particular, 
75 = 1,75 = 75- Note that, as follows from Eq. (34), one has: 

4>l = Vho = ^[(1 - 7 5 )/2]7o = ^7o[(l - 7 5 )/2]7o = ^[(1 + 7s)/2] . 

The matrices P± = (1 ± 7s)/2 are projectors. They satisfy the relations P±P± = P±, 
P±Pzf = 0, P + + P = 1. The P projectors project on fermions of definite chirality. For 
massless particles, chirality coincides with helicity. For massive particles, a chirality +1 
state only coincides with a +1 helicity state up to terms suppressed by powers of m/E. 

The sixteen linearly independent Dirac matrices (r) can be divided into 75-even (r^) 
and 75-odd (To) according to whether they commute or anticommute with 75. For the 
75-even, we have 

4>T E tp = 4> L F E ipR + 4>rF e iPl (r E = 1, 175, o>) , (35) 
whilst for the 75-odd, 

4>T ^ = ipL^otpL + 4>rT ^r (To = 7/,, 7^7s) • (36) 

We see that in a gauge lagrangian fermion kinetic terms and interactions of gauge bosons 
with vector and axial vector fermion currents all conserve chirality while fermion mass terms 
flip chirality. For example, in QED if an electron emits a photon, the electron chirality is 
unchanged. In the ultrarelativistic limit, when the electron mass can be neglected, chirality 
and helicity are approximately the same and we can state that the helicity of the electron 
is unchanged by the photon emission. In a massless theory the LH and the RH fermion 
components are uncoupled and can be transformed separately. If in a gauge theory the LH 
and RH components transform as different representations of the gauge group one speaks 
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of a chiral gauge theory, while if they have the same gauge transformations one has a 
vector gauge theory Thus, QED and QCD are vector gauge theories because, for each 
given fermion, ipL and ipR have the same electric charge and the same colour. Instead, the 
standard EW theory is a chiral theory, in the sense that ipi and tpR behave differently under 
the gauge group (so that parity and charge conjugation non conservation are made possible 
in principle). Thus, mass terms for fermions (of the form ipLipR + h.c.) are forbidden in 
the EW gauge-symmetric limit. In particular, in the Minimal Standard Model (MSM: i.e. 
the model that only includes all observed particles plus a single Higgs doublet), all ipL are 
SU (2) doublets while all ipR are singlets. 

1.6 Quantization of a Gauge Theory 

The lagrangian density Cym in eq.(23) fully describes the theory at the classical level. 
The formulation of the theory at the quantum level requires that procedures of quanti- 
zation, of regularization and, finally, of renormalization are also specified. To start with, 
the formulation of Feynman rules is not straightforward. A first problem, common to all 
gauge theories, including the abelian case of QED, can be realized by observing that the 
free equations of motion for V^ 4 , as obtained from eqs. (21,23), are given by 

[&g llv -d li d v ]V A " = (37) 

Normally the propagator of the gauge field should be determined by the inverse of the 
operator [d 2 g^ u — <9 M <9„] which, however, has no inverse, being a projector over the transverse 
gauge vector states. This difficulty is removed by fixing a particular gauge. If one chooses 
a covariant gauge condition d^V^ = then a gauge fixing term of the form 

AC GF = -^J2\ d "<\ 2 ( 38 ) 

has to be added to the lagrangian (1/A acts as a lagrangian multiplier). The free equations 
of motion are then modified as follows: 

[4,-(l-l/W]^ = 0. (39) 

This operator now has an inverse whose Fourier transform is given by: 

^) = ?Tfc[- 9) - + (l-A)^]«" («) 

which is the propagator in this class of gauges. The parameter A can take any value and it 
disappears from the final expression of any gauge invariant, physical quantity. Commonly 
used particular cases are A = 1 (Feynman gauge) and A = (Landau gauge). 

While in an abelian theory the gauge fixing term is all that is needed for a correct 
quantization, in a non abelian theory the formulation of complete Feynman rules involves 
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a further subtlety. This is formally taken into account by introducing a set of D fictitious 
ghost fields that must be included as internal lines in closed loops (Faddeev-Popov ghosts 
[20]). Given that gauge fields connected by a gauge transformation describe the same 
physics, clearly there are less physical degrees of freedom than gauge field components. 
Ghosts appear, in the form of a transformation Jacobian in the functional integral, in the 
process of elimination of the redundant variables associated with fields on the same gauge 
orbit [21]. By performing some path integral acrobatics the correct ghost contributions 
can be translated into an additional term in the lagrangian density. For each choice of 
the gauge fixing term the ghost langrangian is obtained by considering the effect of an 
infinitesimal gauge transformation V^ c = Vf^ — gC A bc@ A V^ ~ 9^9° on the gauge fixing 
condition. For d^V^ = one obtains: 

d»vf = &>vF - gC ABC d»(6 A Vf) - d?9 c = -[d 2 5 A c + gC ABC V*dy A (41) 

where the gauge condition d^V^ = has been taken into account in the last step. The 
ghost lagrangian is then given by: 

^Cchost = f[d 2 5 AC + gC ABC V*d% A (42) 

where r\ A is the ghost field (one for each index A) which has to be treated as a scalar field 
except that a factor (—1) for each closed loop has to be included as for fermion fields. 

Starting from non covariant gauges one can construct ghost-free gauges. An example, 
also important in other respects, is provided by the set of "axial" gauges: n^V A = where 
n M is a fixed reference 4- vector (actually for n M spacelike one has an axial gauge proper, for 
n 2 = one speaks of a light-like gauge and for n M timelike one has a Coulomb or temporal 
gauge). The gauge fixing term is of the form: 

^gf = ~Y,\^V a \ 2 (43) 

With a procedure that can be found in QED textbooks [19] the corresponding propagator, 
in Fourier space, is found to be: 

D AB (a) = — \-a + n ^ q+7lvq ^ - w Vg" i S ab (u) 

In this case there are no ghost interactions because n^V!^, obtained by a gauge transfor- 
mation from n^V^, contains no gauge fields, once the the gauge condition n^V A = has 
been taken into account. Thus the ghosts are decoupled and can be ignored. 

The introduction of a suitable regularization method that preserves gauge invariance is 
essential for the definition and the calculation of loop diagrams and for the renormalization 
programme of the theory. The method that is by now currently adopted is dimensional 
regularization [22] which consists in the formulation of the theory in n dimensions. All 
loop integrals have an analytic expression that is actually valid also for non integer values 
of n. Writing the results for n = 4 — e the loops are ultraviolet finite for e > and the 
divergences reappear in the form of poles at e = 0. 
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Figure 1: The potential 7 = 1/2 /i 2 M 2 + 
1/4 A(M 2 ) 2 for positive (a) or negative fi 2 (b) 
(for simplicity, M is a 2-dimensional vector). 
The small sphere indicates a possible choice 
for the direction of M. 



1.7 Spontaneous Symmetry Breaking in Gauge Theories 

The gauge symmetry of the SM was difficult to discover because it is well hidden in 
nature. The only observed gauge boson that is massless is the photon. The gluons are 
presumed massless but cannot be directly observed because of confinement, and the W 
and Z weak bosons carry a heavy mass. Indeed a major difficulty in unifying the weak 
and electromagnetic interactions was the fact that e.m. interactions have infinite range 
(m 7 = 0), whilst the weak forces have a very short range, owing to mw,z ^ 0. The solution 
of this problem is in the concept of spontaneous symmetry breaking, which was borrowed 
from condensed matter physics. 

Consider a ferromagnet at zero magnetic field in the Landau-Ginzburg approximation. 
The free energy in terms of the temperature T and the magnetization M can be written 
as 

F(M,T) ~ F (T) + 1/2 /i 2 (T)M 2 + 1/4 A(T)(M 2 ) 2 + ... . (45) 

This is an expansion which is valid at small magnetization. The neglect of terms of higher 
order in M 2 is the analogue in this context of the renormalizability criterion. Also, A(T) > 
is assumed for stability; F is invariant under rotations, i.e. all directions of M in space are 
equivalent. The minimum condition for F reads 

dF/dMi = 0, [/i 2 (T) + A(T)M 2 ]M = . (46) 

There are two cases, shown in Fig. 1. If fi 2 > 0, then the only solution is M = 0, there is 
no magnetization, and the rotation symmetry is respected. In this case the lowest energy 
state (in a quantum theory the vacuum) is unique and invariant under rotations. If /i 2 < 0, 
then another solution appears, which is 

|M | 2 = -fi 2 /X . (47) 

In this case there is a continuous orbit of lowest energy states, all with the same value of 
|M| but different orientations. A particular direction chosen by the vector Mo leads to a 
breaking of the rotation symmetry. 

For a piece of iron we can imagine to bring it to high temperature and to let it melt in 
an external magnetic field B. The presence of B is an explicit breaking of the rotational 
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symmetry and it induces a non zero magnetization M along its direction. Now we lower 
the temperature while keeping B fixed. Both A and [j? depend on the temperature. With 
lowering T, /z 2 goes from positive to negative values. The critical temperature T crit (Curie 
temperature) is where (J?(T) changes sign: /z 2 (T crit ) = 0. For pure iron T crit is below the 
melting temperature. So at T = T crit iron is a solid. Below T crit we remove the magnetic 
field. In a solid the mobility of the magnetic domains is limited and a non vanishing M 
remains. The form of the free energy is again rotationally invariant as in eq.(45). But now 
the system allows a minimum energy state with non vanishing M in the direction where 
B was. As a consequence the symmetry is broken by this choice of one particular vacuum 
state out of a continuum of them. 

We now prove the Goldstone theorem [23]. It states that when spontaneous symmetry 
breaking takes place, there is always a zero-mass mode in the spectrum. In a classical 
context this can be proven as follows. Consider a lagrangian 

C=\\d^-V{<P). (48) 

The potential V(4>) can be kept generic at this stage but, in the following, we will be mostly 
interested in a renormalizable potential of the form (with no more than quartic terms): 

^(0) = -^ 2 2 + ^A0 4 . (49) 

Here by we mean a column vector with real components 0j (1=1, 2. ..N) (complex fields 
can always be decomposed into a pair of real fields), so that, for example, 2 = ^0 2 . 
This particular potential is symmetric under a NxN orthogonal matrix rotation 0' = 0(f), 
where O is a SO(N) transformation. For simplicity, we have omitted odd powers of 0, 
which means that we assumed an extra discrete symmetry under -h- —0. Note that, for 
positive /i 2 , the mass term in the potential has the "wrong" sign: according to the previous 
discussion this is the condition for the existence of a non unique lowest energy state. More 
in general, we only assume here that the potential is symmetric under the infinitesimal 
transformations 

_> ; = + 50, = iseHf^ . (so) 

where 59 A are infinitesimal parameters and t& are the matrices that represent the symmetry 
group on the representation of the fields 0j (a sum over A is understood). The minimum 
condition on V that identifies the equilibrium position (or the vacuum state in quantum 
field theory language) is 

(^/90 4 )(0 l = 0°) = O. (51) 

The symmetry of V implies that 

SV = {dV/d^Sfa = i56 A (dV/d<p i )tf J <p j = . (52) 
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By taking a second derivative at the minimum fa = 0° of both sides of the previous 
equation, we obtain that, for each A: 

The second term vanishes owing to the minimum condition, Eq. (51). We then find 

dV ^- ^HM = . (54) 



dfadfa^ 1 

The second derivatives = (d 2 V/d(pkdfa)(fa = faf) define the squared mass matrix. 
Thus the above equation in matrix notation can be written as 

M 2 t A (j) = . (55) 

In the case of no spontaneous symmetry breaking the ground state is unique, all symmetry 
transformations leave it invariant, so that, for all A, t A <p° = 0. On the contrary, if, for 
some values of A, the vectors (t A (p°) are non-vanishing, i.e. there is some generator that 
shifts the ground state into some other state with the same energy (hence the vacuum is 
not unique), then each t A <p ^ is an eigenstate of the squared mass matrix with zero 
eigenvalue. Therefore, a massless mode is associated with each broken generator. The 
charges of the massless modes (their quantum numbers in quantum language) differ from 
those of the vacuum (usually taken as all zero) by the values of the t A charges: one says 
that the massless modes have the same quantum numbers of the broken generators, i.e. 
those that do not annihilate the vacuum. 

The previous proof of the Goldstone theorem has been given in the classical case. In 
the quantum case the classical potential corresponds to the tree level approximation of 
the quantum potential. Higher order diagrams with loops introduce quantum corrections. 
The functional integral formulation of quantum field theory [21], [24] is the most appropri- 
ate framework to define and compute, in a loop expansion, the quantum potential which 
specifies, exactly as described above, the vacuum properties of the quantum theory. If the 
theory is weakly coupled, e.g. if A is small, the tree level expression for the potential is 
not too far from the truth, and the classical situation is a good approximation. We shall 
see that this is the situation that occurs in the electroweak theory with a moderately light 
Higgs (see Chapter 3, Sec. 3.5). 

We note that for a quantum system with a finite number of degrees of freedom, for 
example one described by the Schroedinger equation, there are no degenerate vacua: the 
vacuum is always unique. For example, in the one dimensional Schroedinger problem with 
a potential: 

V(x) = -/i 2 /2 x 2 + X x 4 /4 , (56) 

there are two degenerate minima at x = ±Xq = (fjp/X) 1 ' 2 which we denote by |+) and | — }. 
But the potential is not diagonal in this basis: the off diagonal matrix elements: 

(+|V|-) = {-\V\+) ~ exp(-khd) = 5 (57) 



18 




X 



Figure 2: A Schroedinger potential V(x) anal- 
ogous to the Higgs potential. 



+ 



are different from zero due to the non vanishing amplitude for a tunnel effect between the 
two vacua given in eq.(57), proportional to the exponential of minus the product of the 
distance d between the vacua and the height h of the barrier with k a constant (see Fig. 
2). After diagonalization the eigenvectors are (|+) + |— ))/v / 2 and (|+) — |— ))/v / 2, with 
different energies (the difference being proportional to 5). Suppose now that you have a sum 
of n equal terms in the potential, V = Yli^i x i)- Then the transition amplitude would 
be proportional to 5 n and would vanish for infinite n: the probability that all degrees 
of freedom together jump over the barrier vanishes. In this example there is a discrete 
number of minimum points. The case of a continuum of minima is obtained, always in the 
Schroedinger context, if we take 



with r = (x, y, z). Also in this case the ground state is unique: it is given by a state with 
total orbital angular momentum zero, an s-wave state, whose wave function only depends 
on |r|, independent of all angles. This is a superposition of all directions with the same 
weight, analogous to what happened in the discrete case. But again, if we replace a single 
vector r, with a vector field M(x), that is a different vector at each point in space, the 
amplitude to go from a minimum state in one direction to another in a different direction 
goes to zero in the limit of infinite volume. In simple words, the vectors at all points in 
space have a vanishing small amplitude to make a common rotation, all together at the 
same time. In the infinite volume limit all vacua along each direction have the same energy 
and spontaneous symmetry breaking can occur. 

A massless Goldstone boson correspond to a long range force. Unless the massless 
particles are confined, as for the gluons in QCD, these long range forces would be easily 
detectable. Thus, in the construction of the EW theory we cannot accept massless physical 
scalar particles. Fortunately, when spontaneous symmetry breaking takes place in a gauge 
theory, the massless Goldstone modes exist, but they are unphysical and disappear from 
the spectrum. Each of them becomes, in fact, the third helicity state of a gauge boson 
that takes mass. This is the Higgs mechanism [8-11] (it should be called Englert-Brout- 
Higgs mechanism, because of the simultaneous paper by Englert and Brout). Consider, for 




(58) 
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example, the simplest Higgs model described by the lagrangian [10,11] 



^ - fiV ' ' p-v/ri ' r-r r 2 



C = -- F 2 + \{d, + ze^Q)0| 2 + /x 2 0*0 - ttW) 2 • (59) 



Note the 'wrong' sign in front of the mass term for the scalar field 0, which is necessary 
for the spontaneous symmetry breaking to take place. The above lagrangian is invariant 
under the U(l) gauge symmetry 

-> A; = A ll - d^ix), 0^0' = exp[ieQ9(x)} 0. (60) 

For the U(l) charge Q we take Qcf) = —0, like in QED, where the particle is e~. Let 
0° = v 7^ 0, with v real, be the ground state that minimizes the potential and induces the 
spontaneous symmetry breaking. In our case v is given by v 2 = fi 2 /X. Making use of gauge 
invariance, we can do the change of variables 

<f>(x) -)• [v + -7^] exp -z- ■ 



V2 S " L Vy/2 1 ' 

A,(x)^A,-d„^W. (61) 
ev v2 

Then the position of the minimum at 0° = v corresponds to h = 0, and the lagrangian 
becomes 

C = --Fl + e 2 v 2 Al + \e 2 h 2 Al + v^/uMj + C(h) . (62) 

The field ((x) is the would-be Goldstone boson, as can be seen by considering only the 
terms in the lagrangian, i.e. setting A^ = in eq.(59). In fact in this limit the kinetic 
term d^d^C, remains but with no ( 2 mass term. Instead, in the gauge case of eq.(59), 
after changing variables in the lagrangian, the field ((x) completely disappears (not even 
the kinetic term remains), whilst the mass term e 2 v 2 A 2 l for A^ is now present: the gauge 
boson mass is M = y/2ev. The field h describes the massive Higgs particle. Leaving a 
constant term aside, the last term in eq.(62) is given by: 

C{h) = -dphPh - h 2 /i 2 + .... (63) 

where the dots stand for cubic and quartic terms in h. We see that the h mass term has 
the "right" sign, due to the combination of the quadratic terns in h that, after the shift, 
arise from the quadratic and quartic terms in 0. The h mass is given by m 2 = 2/i 2 . 

The Higgs mechanism is realized in well-known physical situations. It was actually 
discovered in condensed matter physics by Anderson [25]. For a superconductor in the 
Landau-Ginzburg approximation the free energy can be written as 

F = F + ^B 2 + I ( V — 2«eA)0| 2 /4m - a|0| 2 + /3|0| 4 . (64) 
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Here B is the magnetic field, |0| 2 is the Cooper pair (e~e~) density, 2e and 2m are the 
charge and mass of the Cooper pair. The 'wrong' sign of a leads to 7^ at the minimum. 
This is precisely the non-relativistic analogue of the Higgs model of the previous example. 
The Higgs mechanism implies the absence of propagation of massless phonons (states with 
dispersion relation u = kv with constant v). Also the mass term for A is manifested by the 
exponential decrease of B inside the superconductor (Meissner effect). But in condensed 
matter examples the Higgs field is not elementary but, rather, a condensate of elementary 
fields (like for the Cooper pairs). 

1.8 Quantization of Spontaneously Broken Gauge Theories: 



We have discussed in Sect. (1.6) the problems arising in the quantization of a gauge 
theory and in the formulation of the correct Feynman rules (gauge fixing terms, ghosts 
etc). Here we give a concise account of the corresponding results for spontaneously broken 
gauge theories. In particular we describe the gauge formalism [21], [24], [26]: in this 
formalism the interplay of transverse and longitudinal gauge boson degrees of freedom is 
made explicit and their combination leads to the cancellation from physical quantities of 
the gauge parameter £. We work out in detail an abelian example that later will be easy 
to generalize to the non abelian case. 

We restart from the abelian model of Eq.(59) (with Q = — 1). In the treatment pre- 
sented there the would-be Goldstone boson ((x) was completely eliminated from the la- 
grangian by a non linear field transformation formally identical to a gauge transformation 
corresponding to the U(l) symmetry of the lagrangian. In that description, in the new 
variables we eventually obtain a theory with only physical fields: a massive gauge boson 

with mass M = \/2ev and a Higgs particle h with mass = y/2fx. This is called a 
"unitary" gauge, because only physical fields appear. But if we work out the propagator 
of the massive gauge boson : 



we find that it has a bad ultraviolet behaviour due to the second term in the numerator. 
This choice does not prove to be the most convenient for a discussion of the ultraviolet 
behaviour of the theory. Alternatively one can go to an alternative formulation where 
the would-be Goldstone boson remains in the lagrangian but the complication of keeping 
spurious degrees of freedom is compensated by having all propagators with good ultraviolet 
behaviour ("renormalizable" gauges). To this end we replace the non linear transformation 
for in Eq.(61) with its linear equivalent (after all perturbation theory deals with the small 
oscillations around the minimum): 



Gauges 



.g^u - k^k v /M 2 
1 k 2 -M 2 + ie 



(65) 




(66) 
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Here we have only applied a shift by the amount v and separated the real and imaginary 
components of the resulting field with vanishing vacuum expectation value. If we leave 
as it is and simply replace the linearized expression for 0, we obtain the following quadratic 
terms (those important for propagators): 



£ q uad = - 1 -J2 F ^ FA " U+1 M2A ^ + 



4 

A 

+ l -(d»C? + MA^C + \{d»hY-h^ 2 (67) 

The mixing term between and <9 M £ does not allow to directly write diagonal mass 
matrices. But this mixing term can be eliminated by an appropriate modification of the 
covariant gauge fixing term given in Eq.(38) for the unbroken theory. We now take: 

AC GF = -^A^-iMQ 2 • (68) 

By adding ACqf to the quadratic terms in Eq.(67) the mixing term cancels (apart from a 
total derivative that can be omitted) and we have: 

£ qU ad = -\^f*f a ^ + \m 2 a,a»-^a,) 2 + 



+ 2(W-^C 2 +2(W-^V (69) 
We see that the ( field appears with a mass y/£M and its propagator is: 



> D < = v-w+u- (70) 

The propagators of the Higgs field h and of gauge field A^ are: 

= f-v+fc • (71) 

As anticipated, all propagators have a good behaviour at large k 2 . This class of gauges 
are called " gauges" [26]. Note that for £ = 1 we have a sort of generalization of the 
Feynman gauge with a Goldstone boson of mass M and a gauge propagator: 

Also for £ — > oo the unitary gauge description is recovered in that the would-be Goldstone 
propagator vanishes and the gauge propagator reproduces that of the unitary gauge in 
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Eq.(65). All £ dependence, including the unphysical singularities of the ( and propa- 
gators at k 2 = £M 2 , present in individual Feynman diagrams, must cancel in the sum of 
all contributions to any physical quantity. 

An additional complication is that a Faddeev-Popov ghost is also present in gauges 
(while it is absent in an unbroken abelian gauge theory). In fact under an infinitesimal 
gauge transformation with parameter 9(x) we have the transformations: 

Afj, -> A^-d^O 
(1-ie^ + ^-i^], (74) 

so that: 

8 Ay, = -dy.6, 5h = -e(6, 5( = e9V2(v + -^=) . (75) 
The gauge fixing condition d^A^ — £,M( = undergoes the variation: 

d^A" - iMC -> d^A" - £MC - [d 2 + £M 2 (1 + \)]6 , (76) 

where we used M = \/2ev. From this, recalling the discussion in Sect. (1.6), we see that 
the ghost is not coupled to the gauge boson (as usual for an abelian gauge theory) but has 
a coupling to the Higgs field h. The ghost lagrangian is: 

AC Ghost = v[d 2 + £M 2 (1 + A=)fo . (77) 

v v2 

The ghost mass is seen to be m gh = a/£M and its propagator is: 

w * = e-w + u - (78) 

The detailed Feynman rules follow for all the basic vertices involving the gauge boson, the 
Higgs, the would-be Goldstone boson and the ghost and can be easily derived, with some 
algebra, from the total lagrangian including the gauge fixing and ghost additions. The 
generalization to the non abelian case is in principle straightforward, with some formal 
complications involving the projectors over the space of the would-be Goldstone bosons 
and over the orthogonal space of the Higgs particles. But for each gauge boson that takes 
mass M a we still have a corresponding would-be Goldstone boson and a ghost with mass 
-y/£M a . The Feynman diagrams, both for the abelian and the non abelian listed 
explicitly, for example, in the Cheng and Li textbook in ref. [24]. 

We conclude that the renormalizability of non abelian gauge theories, also in presence 
of spontaneous symmetry breaking, was proven in the fundamental works of t'Hooft and 
Veltman [27] and discussed in detail in [28]. 
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2 QCD: The Theory of Strong Interactions 



2.1 Introduction 

This Chapter is devoted to a concise introduction to Quantum Chromo-Dynamics 
(QCD), the theory of strong interactions [4-6] (for a number of dedicated books on QCD, 
see [29], see also [30]). The main emphasis will be on ideas without too many technicalities. 
As an introduction we present here a broad overview of the strong interactions (for reviews 
of the subject, see, for example, [31,32]). Then some methods of non perturbative QCD 
will be briefly described including both analytic approaches and simulations of the theory 
on a discrete space-time lattice. Then we will proceed to the main focus of the Chapter 
which is on the principles and the applications of perturbative QCD that will be discussed 
in detail. 

As discussed in Chapter 1 the QCD theory of strong interactions is an unbroken gauge 
theory based on the group SU(3) of colour. The eight massless gauge bosons are the 
gluons g£ and the matter fields are colour triplets of quarks q? (in different flavours i). 
Quarks and gluons are the only fundamental fields of the Standard Model (SM) with strong 
interactions (hadrons). The QCD Lagrangian was introduced in Sect. 1.4, eqs. 28- 31. For 
quantization the classical Lagrangian in eq. (28) must be enlarged to contain gauge fixing 
and ghost terms, as described in Chapter 1. The Feynman rules of QCD are listed in Fig. 3. 
The physical vertices in QCD include the gluon-quark-antiquark vertex, analogous to the 
QED photon-fermion-antifermion coupling, but also the 3-gluon and 4-gluon vertices, of 
order e s and e 2 s respectively, which have no analogue in an abelian theory like QED. 

Why SU(N C = 3) co i our l The selection of SU(3) as colour gauge group is unique 
in view of a number of constraints, (a) The group must admit complex representations 
because it must be able to distinguish a quark from an antiquark [33]. In fact there are 
meson states made up of qq but not analogous qq bound states. Among simple groups 
this restricts the choice to SU(N) with N > 3, SO{AN + 2) with N > 2 (taking into 
account that SO (6) has the same algebra as SU(A)) and E(6). (b) The group must admit 
a completely antisymmetric colour singlet baryon made up of 3 quarks: qqq. In fact, 
from the study of hadron spectroscopy we know that the low lying baryons, completing 
an octet and a decuplet of (flavour) SU(3) (the approximate symmetry that rotates the 
3 light quarks u, d and s), are made up of three quarks and are colour singlets. The qqq 
wave function must be completely antisymmetric in colour in order to agree with Fermi 
statistics. Indeed if we consider, for example, a N* ++ with spin z-component +3/2, this 
is made up of (u ff - u ff - u f[~) in an s-state. Thus its wave function is totally symmetric 
in space, spin and flavour so that complete antisymmetry in colour is required by Fermi 
statistics. In QCD this requirement is very simply satisfied by e a b c q a q b q c where a, b, c are 
SU (3) co i our indices, (c) The choice of SU (N c = 3) co ; onr is confirmed by many processes 
that directly measure No- Some examples are listed here. The total rate for hadronic 
production in e + e~ annihilation is linear in Nc- Precisely if we consider R = R e + e - = 
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Figure 3: Feynman rules for QCD. Solid lines represent the quarks, curly lines the gluons, and dotted lines 
the ghosts (see Chapter 1). The gauge parameter is denoted by A. The 3-gluon vertex is written as if all 
gluon lines are outgoing. 
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Figure 4: Comparison of the data on R = a(e + e~ — > hadrons) / a po i n t{e. + e~ — > /i + /i~) with the QCD 
prediction [34] . Nc — 3 is indicated by the data points above 10 GeV (the bb threshold) and 40 GeV 
where the rise due to the Z resonance become appreciable. 



<r(e + e - — > hadrons) /a point (e + e~ — > above the bb threshold and below m z and we 

neglect small computable radiative corrections (that will be discussed later in Sect. 2.7) we 
have a sum of individual contributions (proportional to Q 2 , where Q is the electric charge 
in units of the proton charge) from qq final states with q = u, c, d, s, b: 

R « N c [2-^ + 3-i] « Nc 1 -^ (79) 

The data neatly indicate N c = 3 as seen from Fig. 4 [34] . The slight excess of the data with 
respect to the value 11/3 is due to the QCD radiative corrections ( Sect. 2.7). Similarly 
we can consider the branching ratio B(W~ — > e~u), again in Born approximation. The 
possible fermion-antifermion (//) final states are for / = e~ , r~,d,s (there is no 
f = b because the top quark is too heavy for bi to occur). Each channel gives the same 
contribution, except that for quarks we have Nc colours: 

w-^*-*) - (80) 

For Nc = 3 we obtain B = 11% and the experimental number is B = 10.7%. Another 
analogous example is the branching ratio B(t~ — > e~v e v T ). From the final state channels 
with f = e~ , d we find 

B(t~ e-v e v T ) « (81) 

For Nc = 3 we obtain B = 20% and the experimental number is B = 18% (the less accuracy 
in this case is explained by the larger radiative and phase-space corrections because the 
mass of t~ is much smaller than my/)- An important process that is quadratic in N c is 
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the rate r(7r° — > 27). This rate can be reliably calculated from a solid theorem in field 
theory which has to do with the chiral anomaly: 

IV -> 2 7 ) « (^) 2 ^| = (7.73 ± 0.04)(^) 2 eV (82) 

where the prediction is obtained for /„. = (130.7 ± 0.37) MeV. The experimental result 
is T = (7.7 ± 0.5) eV in remarkable agreement with Nc = 3. There are many more 
experimental confirmations that Nc = 3: for example the rate for Drell-Yan processes (see 
Sect. 2.9 ) is inversely proportional to Nc- 

2.2 Non Perturbative QCD 

The QCD lagrangian in eq. (28) has a simple structure but a very rich dynamical con- 
tent. It gives rise to a complex spectrum of hadrons, implies the striking properties of 
confinement and asymptotic freedom, is endowed with an approximate chiral symmetry 
which is spontaneously broken, has a highly non trivial topological vacuum structure (in- 
stantons, U(1)a symmetry breaking, strong CP violation (which is a problematic item in 
QCD possibly connected with new physics, like axions), ...), an intriguing phase transition 
diagram (colour deconfinement, quark-gluon plasma, chiral symmetry restoration, colour 
superconductivity, ...). 

How do we get testable predictions from QCD? On the one hand there are non pertur- 
bative methods. The most important at present is the technique of lattice simulations (for 
a recent review, see ref. [39]): it is based on first principles, it has produced very valuable 
results on confinement, phase transitions, bound states, hadronic matrix elements and so 
on, and it is by now an established basic tool. The main limitation is from computing 
power and therefore there is continuous progress and a lot of good perspectives for the 
future. Another class of approaches is based on effective lagrangians which provide simpler 
approximations than the full theory, valid in some definite domain of physical conditions. 
Typically at energies below a given scale L particles with mass larger than L cannot be 
produced and thus only contribute short distance effects as virtual states in loops. Under 
suitable conditions one can write down a simplified effective lagrangian where the heavy 
fields have been eliminated (one says "integrated out"). Virtual heavy particle short dis- 
tance effects are absorbed into the coefficients of the various operators in the effective 
Lagrangian. These coefficients are determined in a matching procedure, by requiring that 
the effective theory reproduces the matrix elements of the full theory up to power correc- 
tions. Chiral lagrangians are based on soft pion theorems [40] and are valid for suitable 
processes at energies below 1 GeV (for a recent, concise review see ref. [41] and references 
therein). Heavy quark effective theories [42] are obtained from expanding in inverse pow- 
ers of the heavy quark mass and are mainly important for the study of b and, to less 
accuracy, c decays (for reviews, see, for example, ref. [43]). Soft-collinear effective theories 
(SCET) [44], are valid for processes where quarks with energy much larger than their mass 
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appear. Light energetic quarks not only emit soft gluons, but also collinear gluons (a gluon 
in the same direction as the original quark), without changing their virtuality. In SCET 
the logs associated with these soft and collinear gluons are resummed. The approach of 
QCD sum rules [45,46] has led to interesting results but now appears not to offer much 
potential for further development. On the other hand, the perturbative approach, based 
on asymptotic freedom, still remains the main quantitative connection to experiment, due 
to its wide range of applicability to all sorts of "hard" processes. 

2.2.1 Progress in Lattice QCD 

A main approach to non perturbative problems in QCD is by simulations of the theory 
on the lattice, a technique started by K. Wilson in 1974 [55] which has shown continuous 
progress over the last decades. In this approach the QCD theory is reformulated on a 
discrete space time, an hypercubic lattice of sites (in the simplest realizations) with spacing 
a and 4- volume L 4 ; on each side there are N sites with L = Na. Over the years one 
has learned how to efficiently describe a field theory on a discrete space time and how 
to implement gauge symmetry, chiral symmetry and so on (for a recent review see, for 
example, ref. [39]). Gauge and matter fields are specified on the lattice sites and the 
path integral is computed numerically as a sum over the field configurations. Much more 
powerful computers than in the past now allow for a number of essential improvements. 
As eventually one is interested in the continuum limit, a — > 0, it is important to work 
with as fine lattice spacing a as possible. Methods have been developed for "improving" 
the lagrangian in such a way that the discretization errors vanish faster than linearly in a. 
Larger lattice volume (i. e. large L or N) is also useful as the dimensions of the lattice 
should be as large as possible in comparison with the dimensions of the hadrons to be 
studied. In many cases the volume corrections are exponentially damped, but this is not 
always the case. Lattice simulation is limited to large enough masses of light quarks: in fact, 
heavier quarks have shorter wavelenghts and can be accommodated in a smaller volume. 
In general computations are done for quark and pion masses heavier than in reality and 
then extrapolated to the physical values, but at present one can work with smaller quark 
masses than in the past. One can also take advantage of the chiral effective theory in order 
to control the chiral logs: log(m g /47r/ 7r ) and guide the extrapolation. A recent big step, 
made possible by the availability of more powerful dedicated computers, is the evolution 
from quenched (i.e. with no dynamical fermions) to unquenched calculations. In doing so 
an evident improvement in the agreement of predictions with the data is obtained. For 
example [39], modern unquenched simulations reproduce the hadron spectrum quite well. 
Calculations with dynamical fermions (which take into account the effects of virtual quark 
loops) imply the evaluation of the quark determinant which is a difficult task. How difficult 
depends on the particular calculation method. There are several approaches (Wilson, 
twisted mass, Kogut-Susskind staggered, Ginsparg- Wilson fermions), each with its own 
advantages and disadvantages (including the time it takes to run the simulation on a 
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computer): a compromise between efficiency and theoretical purity is needed. The most 
reliable lattice calculations are today for 2 + 1 light quarks (degenerate up and down 
quarks and a heavier strange quark s). The first calculations for 2 + 1 +1 including charm 
quarks are starting to appear. 

Lattice QCD is becoming increasingly predictive and plays a crucial role in different 
domains. For example, in flavour physics it is essential for computing the relevant hadronic 
matrix elements. In high temperature QCD the most illuminating studies of the phase 
diagram, the critical temperature and the nature of the phase transitions are obtained by 
lattice QCD: as we now discuss the best arguments to prove that QCD implies confinement 
come from the lattice. 

2.2.2 Confinement 

Confinement is the property that no isolated coloured charge can exist but only colour 
singlet particles. Our understanding of the confinement mechanism has much improved 
thanks to lattice simulations of QCD at finite temperatures and densities (for reviews, see, 
for example ref. [35,36]). For example, the potential between a quark and an antiquark has 
been studied on the lattice [37]. It has a Coulomb part at short distances and a linearly 
rising term at long distances: 



with N c the number of colours (N c = 3 in QCD). The scale dependence of a s (the distance 
r is Fourier-conjugate to momentum transfer) will be explained in detail later. The slope 
decreases with increasing temperature until it vanishes at a critical temperature Tc\ then 
above Tq the slope remains zero, as shown in Fig. 5. The value of the critical temperature 
is estimated around Tc ~ 175 MeV. 

The linearly rising term in the potential makes it energetically impossible to separate a 
q — q pair. If the pair is created at one space-time point, for example in e + e~ annihilation, 
and then the quark and the antiquark start moving away from each other in the center of 
mass frame, it soon becomes energetically favourable to create additional pairs, smoothly 
distributed in rapidity between the two leading charges, which neutralize colour and allow 
the final state to be reorganized into two jets of colourless hadrons, that communicate in 
the central region by a number of "wee" hadrons with small energy. It is just like the 
familiar example of the broken magnet: if you try to isolate a magnetic pole by stretching 
a dipole, the magnet breaks down and two new poles appear at the breaking point. 

Confinement is essential to explain why nuclear forces have very short range while 
massless gluon exchange would be long range. Nucleons are colour singlets and they cannot 





where 
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Figure 5: The potential between a quark and an antiquark computed on the lattice in the quenched 
approximation [37]. The upper panel shows that the slope of the linearly rising term decreases with 
temperature and, at the critical temperature Tc, it vanishes. At T > Tc the slope remains put at zero 
(lower panel). 

exchange colour octet gluons but only colourless states. The lightest colour singlet hadronic 
particles are pions. So the range of nuclear forces is fixed by the pion mass r ~ m^ 1 ~ 
10 -13 cm : V ~ exp(— m n r)/r. 

The phase transitions of colour deconfinement and of chiral restauration appear to 
happen together on the lattice [35,36,39] (see Fig. 6). A rapid transition is observed in 
lattice simulations where the energy density e(T) is seen to sharply increase near the critical 
temperature for deconfinement and chiral restauration (see Fig. 7). The critical parameters 
and the nature of the phase transition depend on the number of quark flavours nf and on 
their masses (see Fig.8). For example, for rif = 2 or 2+1 (i.e. 2 light u and d quarks and 1 
heavier s quark), T c ~ 175 MeV and e(T c ) ~ 0.5— 1.0 GeV/fm 3 . For realistic values of the 
masses m s and m u ^ the two phases are connected by a smooth crossover, while the phase 
transition becomes first order for very small or very large m u ^, s - Accordingly the hadronic 
phase and the deconfined phase are separated by a crossover region at small densities and 
by a critical line at high densities that ends with a critical point. Determining the exact 
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location of the critical point in T and hb is an important challenge for theory and is also 
important for the interpretation of heavy ion collision experiments. At high densities the 
colour superconducting phase is also present with bosonic diquarks acting as Cooper pairs. 




Figure 6: Order parameters for deconfinement (bottom) and chiral symmetry restoration (top), as a 
function of temperature [36,39]. On a finite lattice the singularities associated to phase transitions are 
not present but their development is indicated by a rapid rate of change. The vacuum expectation value 
of the quark-antiquark condensate, with increasing temperature, goes from the finite value that breaks 
chiral symmetry down to zero where chiral symmetry is restaured. In a comparable temperature range, 
the Wilson plaquette, the order parameter for deconfinement, goes from zero to a finite value. 

A large investment is being done in experiments of heavy ion collisions with the aim of 
finding some evidence of the quark gluon plasma phase. Many exciting results have been 
found at the CERN SPS in the past years, more recently at RHIC and now at the LHC in 
dedicated heavy ion runs [38] (the ALICE detector is especially designed for the study of 
heavy ion collisions). 

2.2.3 Chiral Symmetry in QCD and the Strong CP Problem 

In the QCD lagrangian, eq. (28), the quark mass terms are of the general form [mipLipR 
+h.c] (recall the definition of ipL,R in Sect. 1.5 and the related discussion). These terms 
are the only ones that show a chirality flip. In the absence of these terms, i.e. for m = 0, 
the QCD lagrangian would be invariant under independent unitary tranformations sep- 
arately on ipL and on ipR. Thus, if the masses of the Nf lightest quarks are neglected 
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Figure 7: The energy density divided by the 4rth power of the temperature, computed on the lattice with 
different number of sea flavours, shows a marked rise near the critical temperature [36,39]. The arrows 
on top show the limit for a perfect Bose gas (while the hot dense hadronic fluid is not expected to be a 
perfect gas). 

the QCD lagrangian is invariant under a global U(Nf)i <S> U{Nf)n chiral group. Consider 
Nf = 2: SU(2)v corresponds to the observed approximate isospin symmetry and U(l)y 
to the portion of baryon number associated with u and d quarks. Since no approximate 
parity doubling of light quark bound states is observed the U(2)a symmetry must be spon- 
taneously broken (for example, no opposite parity analogues of protons and neutrons exist 
with a few tens of MeV separation in mass from the ordinary nucleons). The breaking 
of chiral symmetry is induced by the VEV of a quark condensate: for Nf = 2 this is 
[ulUr + c?iC?K+h.c.]. A recent lattice calculation [49] has given for this condensate the 
value [234 ± 18 MeVf (in MS, N f = 2 + 1, with the physical m s value, at the scale of 
2 GeV). This scalar operator is an isospin singlet, so that it preserves U(2)v but breaks 
U (2) a (it transforms like (1/2,1/2) under U(2)l<^U (2)r but is a singlet under the diagonal 
group U (2)y). The pseudoscalar mesons are obvious candidates for the would-be Goldstone 
bosons associated with the breakdown of the axial group in that they have the quantum 
number of the broken generators: the three pions are the approximately massless Goldstone 
bosons (exactly massless in the limit of vanishing u and d quark masses) associated with 
the breaking of three generators of U(2) L (g)U(2) R down to SU(2) V (g) U(l) v (g) U(1) A - 
The couplings of Goldstone bosons are very special: in particular only derivative couplings 
are allowed. The pions as pseudo-Goldstone bosons have couplings that satisfy strong con- 
straints. An effective chiral lagrangian formalism [40] allows to systematically reproduce 
the low energy theorems implied by the approximate status of Goldstone particles for the 
pion and successfully describes QCD at energies at scales below ~ 1 GeV. 

The breaking mechanism for the remaining U(1)a arises from an even subtler mecha- 
nism. A state in the 77 — 77' space cannot be the associated Goldstone particle because the 
masses are too large [48] and the 77' mass does not vanish in the chiral limit [50] . Rather the 
conservation of the singlet axial current j§ = 2^7^75^ is broken by the Adler-Bell-Jackiw 
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Figure 8: Left: A schematic view of the QCD phase diagram. Right: On the lattice the nature of the 
phase transition depends on the number of quark flavours and their masses as indicated [39] 



anomaly [47]: 

d,f, = I(x) =N f °^J2 F ^ pA " U = N f^ Tr ( F ^) (85) 

A 

(recall that = T,F A t A and the normalization is Tr(t A t B ) = 1/25 AB ) with F A given in 
eq.31 and j£ being the u+d singlet axial current (the factor of Nf, in this case Nf = 2, in 
front of the right hand side takes into account that Nf flavours are involved) and 

K = \^ pApa ( 86 ) 

An important point is that the pseudoscalar quantity I(x) is a four divergence. Precisely 
one can check that: 

Tr(F^F^) = d% (87) 

with 

K = e, vXa Tr[A»{¥ x ° - -ie s A x A°)] (88) 

As a consequence the modified current j£ and its associated charge Q 5 appear to be still 
conserved: 

d^ = 9M-N f ^) = (89) 

and could act as modified chiral current and charge with an additional gluonic component. 
But actually this charge is not conserved due to the topological structure of the QCD 
vacuum (instantons) as discussed in the following (for an introduction, see ref. [51]). 

The configuration where all gauge fields are zero A A = can be called "the vacuum". 
However all configurations that are connected to A A = by a gauge transformation do 
also correspond to the same physical vacuum. For example, in a abelian theory all gauge 
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fields that can be written as the gradient of a scalar, = d^x( x ) are equivalent to 
= 0. In non abelian gauge theories there are some "large" gauge transformations 
that are topologically non trivial and correspond to non vanishing integer values of a 
topological charge, the "winding number". Taking SU(2) for simplicity (in QCD it could 
be any such subgroup of colour SU(3)) we can consider the following, time independent 
gauge transformation: 

nx(5) = f2 z Jrp (90) 

where d is a positive constant. Note that fl^ 1 = fi*. Starting from = (A , AA = (0, 0) 
(i=l,2,3) with A M = SA^r a /2 the gauge transformed potential by Q\ is (recall the general 
expression of a gauge transformation in eq. 15): 



A| 1} = --[y^i^)]^ 1 ^) (91) 

e 



For the vector potential A^ 1 ), which, being a pure gauge, is part of the "vacuum", the 
winding number n, defined in general by 

,3 



n 
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J d 3 xTr[Ai(x)Aj(x)A k (x)]e^ (92) 

is equal to 1: n = 1. Similarly, for obtained from Vt m = [f2i] m one has n = m. Given 

eq. 87 we could expect that the integrated four-divergence would vanish but instead one 
finds: 

J d*x Tr(F^F^) = ^Jd 4 x djf = d 3 x k ]1l™ = n + - n_ (93) 

for a configuration of gauge fields that vanish fast enough on the space sphere at infinity 
and the winding numbers are at time t = — , +oo ("instantons"). 

From the above discussion it follows that in QCD all gauge fields can be classified in 
sectors with different n: there is a vacuum for each n, \n), and Qi\n) = \n + 1) (not gauge 
invariant!). The true vacuum must be gauge invariant (up to a phase) and is obtained as 
a superposition of all \n): 

\9) = S+^e~ in V) (94) 

In fact: 

fl x |0) = ^e- ine \n + l) =e i9 \6) (95) 
If we compute the expectation value of any operator O in the 9 vacuum we find: 

(6\0\9) = S min e i(m - n)9<m| ° |n) (96) 

The path integral describing the O vacuum matrix element at 9 = must be modified to 
reproduce the extra phase, taking eq. 93 into account: 

(9\0\9) = J dAd^dip O exp[iS QC D + ^ J d 4 x Tr(F^F^)] (97) 
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This is equivalent to adding a 9 term to the QCD lagrangian: 



Cqcv = 6 ^j d ' x Tr(F^ u fn (98) 

The 9 term is parity P odd and charge conjugation C even, so that it introduces CP 
violation in the theory (and also time reversal T violation). A priori one would expect 9 
to be o(l). But it would contribute to the neutron electric dipole moment, according to 
d n (e'cm) ~ 3 10~ 16 #. The strong experimental bounds on d n ( d n (e'cm) < 3 10~ 26 [14]) 
imply that 9 must be very small: 9 < 1CT 10 . The so-called "Strong CP-problem" or "9- 
problem" is to find an explanation for such a small value [51,52]. An important point that 
is relevant for a possible solution is that a chiral transformation translates 9 by a fixed 
amount. By recalling eq. 89 we have: 

e iSQs\0} = \0-2N f 6) (99) 

To prove this relation we first observe that Q§ is not gauge invariant under Qi because it 
involves k^. 

Q^sQ- 1 = Q 5 - Q^Nf — l [ d 3 x k Q^ = Q 5 - 2N f (100) 

47T J 

It then follows that 

^l ie i& ^\9) = fixe^nr^ilfl) = e< e ~ 2N ^e iS ^\9) (101) 

which implies eq. 99. Thus in a chiral invariant theory one could dispose of 9. For this it 
would be sufficient that a single quark mass is zero and the obvious candidate would be 
m u = 0. But apparently this possibility has been excluded [52]. For non vanishing quark 
masses the transformation m — > \j\mU R needed to make the mass matrix hermitian (which 
implies 75 - free) and diagonal involves a chiral transformation that affects 9. Considering 
that U(N) = U(l) SU(N), that for hermitian m the argument of the determinant 
vanishes: ArgDetm = the transformation from a generic ml to a real and diagonal m 
gives: 

ArgDet m = = ArgDet U* L + ArgDet m + ArgDet U R = -2N f (6 L - 5 R ) + ArgDet m! 

(102) 

From this equation one derives the phase (Sr — Sl) of the chiral transformation and then, 
by eq. 99, the important result for the effective 9 value: 

9 eff = 9 + ArgDet m (103) 

As we have seen the small empirical value of 9 e ff poses a serious naturalness problem for 
the SM. Among the possible solutions perhaps the most interesting option is a mechanism 
proposed by Peccei-Quinn [53] . One assumes that the SM or an enlarged theory is invariant 
under an additional chiral symmetry U(1)pq acting on the fields of the theory. This 
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symmetry is spontaneously broken by the vacuum expectation value vpq of a scalar field. 
The associated Goldstone boson, the axion, is actually not massless because of the chiral 
anomaly. The parameter 9 is canceled by the vacuum expectation value of the axion field 
due to the properties of the associated potential, also determined by the anomaly. Axions 
could contribute to the Dark Matter in the Universe, if their mass falls in a suitable narrow 
range (for a recent review, see, for example, [54]. Alternative solutions to the ^-problem 
have also been suggested. Some of them can probably be discarded (for example, that 
the up quark is exactly massless), while other ones are still possible: for example, in 
supersymmetric theories, if the smallness of 9 could be guaranteed at the Planck scale by 
some feature of the more fundamental theory valid there, then the non rinormalization 
theorems of supersymmetry would preserve its small value throughout the running down 
to low energy. 

2.3 Massless QCD and Scale Invariance 

As discussed in Chapter 2, the QCD lagrangian in eq. (28) only specifies the theory at 
the classical level. The procedure for quantization of gauge theories involves a number of 
complications that arise from the fact that not all degrees of freedom of gauge fields are 
physical because of the constraints from gauge invariance which can be used to eliminate 
the dependent variables. This is already true for abelian theories and one is familiar with 
the QED case. One introduces a gauge fixing term (an additional term in the lagrangian 
density that acts as a Lagrange multiplier in the action extremization) . One can choose to 
preserve manifest Lorentz invariance. In this case, one adopts a covariant gauge, like the 
Lorentz gauge, and in QED one proceeds according to the formalism of Gupta-Bleuler [19] . 
Or one can give up explicit formal covariance and work in a non covariant gauge, like 
the Coulomb or the axial gauges, and only quantize the physical degrees of freedom (in 
QED the transverse components of the photon field). While this is all for an abelian 
gauge theory, in the non-abelian case some additional complications arise, in particular 
the necessity to introduce ghosts for the formulation of Feynman rules. As we have seen, 
there are in general as many ghost fields as gauge bosons and they appear in the form 
of a transformation Jacobian in the Feynman functional integral. Ghosts only propagate 
in closed loops and their vertices with gluons can be included as additional terms in the 
lagrangian density which are fixed once the gauge fixing terms and their infinitesimal gauge 
transformations are specified. Finally the complete Feynman rules in either the covariant 
or the axial gauges can be obtained and they appear in Fig. 3. 

Once the Feynman rules are derived we have a formal perturbative expansion but loop 
diagrams generate infinities. First a regularization must be introduced, compatible with 
gauge symmetry and Lorentz invariance. This is possible in QCD. In principle one can in- 
troduce a cut-off K (with dimensions of energy), for example, a' la Pauli-Villars [19]. But at 
present the universally adopted regularization procedure is dimensional regularization that 
we will briefly describe later on. After regularization the next step is renormalization. In 
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a renormalizable theory (which is the case for all gauge theories in 4 spacetime dimensions 
and for QCD in particular) the dependence on the cutoff can be completely reabsorbed in 
a redefinition of particle masses, of gauge coupling(s) and of wave function normalizations. 
After renormalization is achieved the perturbative definition of the quantum theory that 
corresponds to a classical lagrangian like in eq. (28) is completed. In the QCD Lagrangian 
of eq. (28) quark masses are the only parameters with physical dimensions (we work in the 
natural system of units h = c = 1). Naively we would expect that massless QCD is scale 
invariant. This is actually true at the classical level. Scale invariance implies that dimen- 
sionless observables should not depend on the absolute scale of energy but only on ratios 
of energy-dimensional variables. The massless limit should be relevant for the asymptotic 
large energy limit of processes which are non singular for m — > 0. 

The naive expectation that massless QCD should be scale invariant is false in the 
quantum theory. The scale symmetry of the classical theory is unavoidably destroyed by 
the regularization and renormalization procedure which introduce a dimensional parameter 
in the quantum version of the theory. When a symmetry of the classical theory is necessarily 
destroyed by quantization, regularization and renormalization one talks of an "anomaly". 
So, in this sense, scale invariance in massless QCD is anomalous. 

While massless QCD is finally not scale invariant, the departures from scaling are 
asymptotically small, logarithmic and computable. In massive QCD there are additional 
mass corrections suppressed by powers of m/E, where E is the energy scale (for processes 
that are non singular in the limit m — » 0). At the parton level (q and g) we can conceive 
to apply the asymptotic predictions of massless QCD to processes and observables (we use 
the word "processes" for both) with the following properties ("hard processes"), (a) All 
relevant energy variables must be large: 

Ei = ZiQ, Q » nij] zo scaling variables o(l) (104) 

(b) There should be no infrared singularities (one talks of "infrared safe" processes), (c) 
The processes concerned must be finite for m — > (no mass singularities). To possibly 
satisfy these criteria processes must be as "inclusive" as possible: one should include all 
final states with massless gluon emission and add all mass degenerate final states (given 
that quarks are massless also q — q pairs can be massless if "collinear", that is moving 
together in the same direction at the common speed of light). 

In perturbative QCD one computes inclusive rates for partons (the fields in the la- 
grangian, that is, in QCD, quarks and gluons) and takes them as equal to rates for hadrons. 
Partons and hadrons are considered as two equivalent sets of complete states. This is called 
"global duality" and it is rather safe in the rare instance of a totally inclusive final state. 
It is less so for distributions, like distributions in the invariant mass M ("local duality") 
where it can be reliable only if smeared over a sufficiently wide bin in M. 

Let us discuss more in detail infrared and collinear safety. Consider, for example, a 
quark virtual line that ends up into a real quark plus a real gluon (Fig. 9). 
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Figure 9: The splitting of a virtual quark into a quark and a gluon. 



For the propagator we have: 

propagator = — = — — — = — — • (105) 

p p 5 (p + k) 2 -m 2 2(p-k) 2E k E p l-f3 p cos9 v ' 

Since the gluon is massless, Ef, can vanish and this corresponds to an infrared singularity. 
Remember that we have to take the square of the amplitude and integrate it over the 
final state phase space, or, in this case, all together, dE^jE^. Indeed we get 1/E% from 
the squared amplitude and d 3 k/Ek ~ E^dE^ from the phase space. Also, for m — > 0, 



/3 P = Jl - m 2 1 'Ep 1 and (1 — /3 p cos8) vanishes at cos 8 = 1. This leads to a collinear 
mass singularity. 

There are two very important theorems on infrared and mass singularities. The first 
one is the Bloch-Nordsieck theorem [56]: infrared singularities cancel between real and 
virtual diagrams (see Fig. 10) when all resolution indistinguishable final states are added 
up. For example, for each real detector there is a minimum energy of gluon radiation that 
can be detected. For the cancellation of infrared divergences, one should add all possible 
gluon emission with a total energy below the detectable minimum. The second one is 
the Kinoshita-Lee, Nauenberg theorem [57]: mass singularities connected with an external 
particle of mass m are canceled if all degenerate states (that is with the same mass) are 
summed up. That is for a final state particle of mass m we should add all final states 
that in the limit m — Y have the same mass, also including gluons and massless pairs. 
If a completely inclusive final state is taken, only the mass singularities from the initial 
state particles remain (we shall see that they will be absorbed inside the non perturbative 
parton densities, which are probability densities of finding the given parton in the initial 
hadron) . 







Figure 10: The diagrams contributing to the total cross-section e + e — > hadrons at order a s . For simplicity, 
only the final state quarks and (virtual or real) gluons are drawn. 

Hard processes to which the massless QCD asymptotics can possibly apply must be 
infrared and collinear safe, that is they must satisfy the requirements from the Bloch- 
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Nordsieck and the Kinoshita-Lee-Nauenberg theorems. We give now some examples of 
important hard processes. One of the simplest hard processes is the totally inclusive cross 
section for hadron production in e + e~ annihilation, Fig. 11, parameterized in terms of 
the already mentioned dimensionless observable R = <r(e + e~ — > hadrons)/a point {e + e~ — >• 
H + H~). The pointlike cross section in the denominator is given by o v&mt = Ana 2 /3s, where 
s = Q 2 = 4E 2 is the squared total center of mass energy and Q is the mass of the exchanged 
virtual gauge boson. At parton level the final state is (qq + n g + n' q'q') and n and n' are 
limited at each order of perturbation theory. It is assumed that the conversion of partons 
into hadrons does not affect the rate (it happens with probability 1). We have already 
mentioned that in order for this to be true within a given accuracy an averaging over a 
sufficiently large bin of Q must be understood. The binning width is larger in the vicinity 
of thresholds: for example when one goes across the charm cc threshold the physical cross- 
section shows resonance bumps which are absent in the smooth partonic counterpart which 
however gives an average of the cross-section. 




Figure 11: The total cross-section e + e — > hadrons. 

A very important class of hard processes is Deep Inelastic Scattering (DIS) 

I + N->1' + X I = e ± ,/i ± ,z/,z/ (106) 

which has played and still plays a very important role for our understanding of QCD and 
of nucleon structure. For the processes in eq. (106), shown in Fig. 12, we have, in the lab 
system where the nucleon of mass m is at rest: 

O 2 

Q 2 = -q 2 = -(k-k 1 ) 2 = AEE' sin 2 9/2; mv = (p.q); x = ^ (107) 

In this case the virtual momentum q of the gauge boson is spacelike, x is the familiar 
Bjorken variable. The DIS processes in QCD will be extensively discussed in Sect. 2.8 

2.4 The Renormalization Group and Asymptotic Freedom 

In this section we aim at providing a reasonably detailed introduction to the renor- 
malization group formalism and the concept of running coupling which leads to the result 
that QCD has the property of asymptotic freedom. We start with a summary on how 
renormalization works. 
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Figure 12: Deep inelastic lepto-production. 



In the simplest conceptual situation imagine that we implement regularization of di- 
vergent integrals by introducing a dimensional cut-off K that respects gauge and Lorentz 
invariance. The dependence of renormalized quantities on K is eliminated by absorbing it 
into a redefinition of m (the quark mass: for simplicity we assume a single flavour here), 
the gauge coupling e (can be e in QED or e s in QCD) and the wave function renormal- 

1 /2 

ization factors Z q ; g for q and g, using suitable renormalization conditions (that is precise 
definitions of m, g and Z that can be implemented order by order in perturbation theory). 
For example we can define the renormalized mass m as the position of the pole in the quark 
propagator and, similarly, the normalization Z q as the residue at the pole: 

2 

Propagator = — — - — - + no — pole terms (108) 
p 2 — m 2 

The renormalized coupling e can be defined in terms of a renormalized 3-point vertex at 
some specified values of the external momenta. Precisely, we consider a one particle irre- 
ducible vertex (1PI). We recall that a connected Green function is the sum of all connected 
diagrams, while 1PI Green functions are the sum of all diagrams that cannot be separated 
into two disconnected parts by cutting only one line. 

We now become more specific by concentrating on the case of massless QCD. If we 
start from a vanishing mass at the classical (or "bare") level, tuq = 0, the mass is not 
renormalized because it is protected by a symmetry, chiral symmetry. The conserved 
currents of chiral symmetry are axial currents: (p/^sq. The divergence of the axial current 
gives, by using the Dirac equation, d^lqj^q) = Imq^^q. So the axial current and the 
corresponding axial charge are conserved in the massless limit. Actually the singlet axial 
current is not conserved due to the anomaly, but, since QCD is a vector theory we have 
not to worry about chiral anomalies in the present context. As there are no 75 around 
the chosen regularization preserves chiral symmetry besides gauge and Lorentz symmetry 
and the renormalized mass remains zero. The renormalized propagator has the form in 
eq. (108) with m — 0. 

The renormalized coupling e s can be defined from the renormalized 1PI 3-gluon vertex 
at a scale — /i 2 (Fig. 13): 

V bare (p 2 ,q 2 ,r 2 ) = ZV ren (p 2 ,q 2 ,r 2 ), Z = Z; 3 / 2 , V ren ( V, -/i 2 7 -fi 2 ) -+ e s (109) 
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We could as well use the quark-gluon vertex or any other vertex which coincides with 
e s o in lowest order (even the ghost-gluon vertex, if we want). With a regularization and 
renormalization that preserves gauge invariance we are guaranteed that all these different 
definitions are equivalent. 



+ ... 




Figure 13: Diagrams contributing to the 1PI 3-gluon vertex at the one-loop approximation level. 

Here Vb are is what is obtained from computing the Feynman diagrams including, for 
example, the 1-loop corrections at the lowest non trivial order. Vb are is defined as the scalar 
function multiplying the 3-gluon vertex tensor (given in Fig. 3), normalized in such a way 
that it coincides with e s o in lowest order. Hare contains the cut-off K but does not know 
about /i. Z is a factor that depends both on the cut-off and on \x but not on momenta. 
Because of infrared singularities the defining scale \x cannot vanish. The negative value 
— /i 2 < is chosen to stay away from physical cuts (a gluon with negative virtual mass 
cannot decay). Similarly, in the massless theory, we can define Z~ l as the inverse gluon 
propagator (the 1PI 2-point function) at the same scale — fi 2 (the vanishing mass of the 
gluon is guaranteed by gauge invariance). 

After computing all 1-loop diagrams indicated in Fig. 13 we have: 

K 2 

V bar e(p 2 ,p 2 ,p 2 ) = e s0 [1 + C« s0 " log — + ...] = 

K 2 -/i 2 

= [l + ca s -log + ...]e s0 [l + ca s0 ■ log— r-] 

— fi 2 p 1 

= Z v 1 e s0 [l + ca s -log— H 

K 2 -a 2 

= [l + da s -log + ...]e s [l + ca s ■ log— r-] 

—fi z p z 

= z/i 2 v ren (110) 

Note the replacement of a S Q with a s in the second step, as we work at 1-loop accuracy. 
Then we change e^o into e s given by eo = Z^ 3 ^ 2 Zye and this implies changing c into 
d in the first bracket. The definition of e s demands that one precisely specifies what 
is included in Z. For this, in a given renormalization scheme, a prescription is fixed 
to specify the finite terms that go into Z (i.e. the terms of order a s that accompany 
logfT 2 ). Then V ren is specified and the renormalized coupling is defined from it according 
to eq. (109). For example, in the momentum subtraction scheme we define V ren (p 2 , p 2 , p 2 ) = 
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e s + Vbare(p 2 , P 2 , P 2 ) — Vbare(—H 2 , —f^ 2 , —fJ 2 ), which is equivalent to say, at 1-loop, that all 
finite terms that do not vanish at p 2 = —fi 2 are included in Z. 

A crucial observation is that Hare depends on K but not on /i, which is only introduced 
when Z, V ren and hence a s are defined. (From here on, for shorthand, we write a to indicate 
either the QED coupling or the QCD coupling a s ). More in general for a generic Green 
function G, we similarly have: 



so that we have: 



or 



G bare (K ,a ,Pi) = Z G G ren (n ,a,p t ) (111) 
dGb are d 



d log /i 2 d log fx 2 



[Z G G ren ] = (112) 



Zcl^ + + ^^h]G ren = (113) 

a log fx 2 a log n 2 da Z G a log fi 2 



Finally the renormalization group equation (RGE) can be written as: 

r 9 . d 

+ /3q- 

alog/i 2 oa 

where 



[ aT~TT2 + ^(«)lr + ^(a)]G ren = (114) 



and 



a log n 2 
d\ogZ G 

lG{a) = ~d\^ (U6) 
Note that (3(a) does not depend on which Green function G we are considering; actually 
it is a property of the theory and of the renormalization scheme adopted, while 7g(q;) also 
depends on G. Strictly speaking the RGE as written above is only valid in the Landau 
gauge (A = 0). In other gauges an additional term that takes the variation of the gauge 
fixing parameter A should also be included. We omit this term, for simplicity, as it is not 
relevant at the 1-loop level. 

Assume that we want to apply the RGE to some hard process at a large scale Q, 
related to a Green function G that we can always take as dimensionless (by multiplication 
by a suitable power of Q). Since the interesting dependence on Q will be logarithmic we 
introduce the variable t as : 

O 2 

t = log^ (117) 
[i 2 

Then we can write G ren = F(t, a, Xi) where Xi are scaling variables (we often omit to write 
them in the following). In the naive scaling limit F should be independent of t, according to 
the classical intuition that massless QCD is scale invariant. To find the actual dependence 
on t, we want to solve the RGE 
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with a given boundary condition at t — (or Q 2 = /i 2 ): F(0, a). 

We first solve the RGE in the simplest case that 7g(«) = 0. This is not an unphysical 
case: for example, it applies to R = R e + e - = a(e + e~ —> hadrons) / a po i nt (e + e~ — > 
where the vanishing of 7 is related to the non renormalization of the electric charge in 
QCD (otherwise the proton and the electron charge would not exactly compensate each 
other: this will explained in Sect. 2.7). So we consider the equation: 



The solution is simply 



F(t,a) = F[0,a{t)) (120) 
where the "running coupling" a(t) is defined by: 

1 = L wf (121) 

Note that from this definition it follows that a(0) = a, so that the boundary condition is 
also satisfied. To prove that F[0, ait)} is indeed the solution, we first take derivatives with 
respect of t and a (the two independent variables) of both sides of eq. (121). By taking 
d/dt we obtain 

1 = (122) 

We then take d/da and obtain 

1 1 da(t) . , 

(3{a) (3{a{t) da 1 } 

These two relations make explicit the dependence of the running coupling on t and a: 

da(t) 



dt 

da(t) P(ct(t)) 
da (3(a) 



P(a(t)) (124) 

(125) 



Using these two equations one immediately checks that F[0, a(t)} is indeed the solution. 

Similarly, one finds that the solution of the more general equation with 7 7^ 0, eq. (118), 
is given by: 

r a{t) 'y(a') 

F{t, a) = MO, <*(*)] exp / -wr^da' (126) 

In fact the sum of the two derivatives acting on the factor F[0, a(t)] vanishes and the 
exponential is by itself a solution of the complete equation. Note that the boundary 
condition is also satisfied. 
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The important point is the appearance of the running coupling that determines the 
asymptotic departures from scaling. The next step is to study the functional form of the 
running coupling. From eq. (124) we see that the rate of change with t of the running 
coupling is determined by the (3 function. In turn (3(a) is determined by the \x dependence 
of the renormalized coupling through eq. (115). Clearly there is no dependence on \i of 
the basic 3-gluon vertex in lowest order (order e). The dependence starts at 1-loop, that 
is at order e 3 (one extra gluon has to be emitted and reabsorbed). Thus we obtain that in 
perturbation theory: 

<h oc e 3 (127) 



d log /i 2 

Recalling that a = e 2 /47i~, we have: 

da de A , . 

ex 2e- oc e 4 oca 2 128 

a log n z a log n z 

Thus the behaviour of /3(a) in perturbation theory is as follows: 

(3(a) = ±ba 2 [l + b'a + ...] (129) 

Since the sign of the leading term is crucial in the following discussion, we stipulate that 
always b > and we make the sign explicit in front. 

Let us make the procedure more precise for computing the 1-loop beta function in QCD 
(or, similarly, in QED). The result of the Hoop 1PI diagrams for V ren can be written down 
as: 

u 2 

V ren = e[l + aBz g \ogf-^\ (130) 



Vren satisfies the RGE: 



' + WniHI ~ -iMWren = (131) 



d log /i 2 da de 2 

With respect to eq. (114) the beta function term has been rewritten taking into account 

— 1/2 

that V ren starts with e and the anomalous dimension term arises from a factor Z g for 
each gluon leg. In general for a n-leg 1PI Green function V n fiare = Z g n ^ 2 V n , reni if all 
external legs are gluons. Note that in the particular case of V = V 3 that is used to define 
e other Z factors are absorbed in the replacement Z^Zg^e^ = e. At 1-loop accuracy we 
replace (3(a) = —ba 2 and r ) g (a) = jg ce. All together one obtains: 

b = 2(B 3g - | 7 «) (132) 

Similarly we can write the diagrammatic expression and the RGE for the 1PI 2-gluon 
Green function which is the inverse gluon propagator II (a scalar function after removing 
the gauge invariant tensor): 

u 2 

U ren = [l + aB 2g \ogf^ + ...} (133) 
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and 

+ p( a )°- - 7j(a) ]n„„ = (134) 

Notice that the normalization and the phase of II are specified by the lowest order term 
being 1. In this case the /3 function term is negligible being of order a 2 (because IT is a 
function of e only through a) and we obtain: 

Tf° = B 2g (135) 

Thus, finally: 

b = 2(B 3g - *-B 2g ) (136) 
By direct calculation at 1-loop one finds: 

N O 2 

QED : (3(a) ~ + ha 2 + b = ^ (137) 

i 

where Nc = 3 for quarks and Nc = 1 for leptons and the sum runs over all fermions of 
charge QiC that are coupled. Also, one finds: 

QCD : (3(a) ~ -ha 2 + b = 11Nc ~ 2llf ( 138 ) 

Yin 

where, as usual, n/ is the number of coupled (see below) flavours of quarks (we assume 
here that rif < 16 so that b > in QCD). If a(t) is small we can compute (3(a(t)) in 
perturbation theory. The sign in front of b then decides the slope of the coupling: a(t) 
increases with t (or Q 2 ) if (3 is positive at small a (QED), or a(t) decreases with t (or Q 2 ) 
if /3 is negative at small a (QCD). A theory like QCD where the running coupling vanishes 
asymptotically at large Q 2 is called (ultraviolet) "asymptotically free". An important result 
that has been proven [58] is that in 4 spacetime dimensions all and only non-abelian gauge 
theories are asymptotically free. 

Going back to eq. (121) we replace /3(a) ~ ± ba 2 , do the integral and perform a 
simple algebra. We find 

« ED : «M ~ T^bl (139) 



and 



QCD : a(t) ~ (140) 



A slightly different form is often used in QCD. Defining 1/a = 61og/i /A.q CD we can 
write: 

°(*) ~ 1 r = 2 F^J = r>2 — ( 141 ) 

5 + bt blogjf- + Mogg blogj^- 

QCD r JV QCD 

The parameter \i has been traded for the parameter Aqcd- We see that a(t) decreases 
logarithmically with Q 2 and that one can introduce a dimensional parameter Aqcd that 
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replaces fi. Often in the following we will simply write A for Aqqd- Note that it is clear 
that A depends on the particular definition of a, not only on the defining scale \i but also on 
the renormalization scheme (see, for example, the discussion in the next section). Through 
the parameter b, and in general through the (3 function, it also depends on the number rif 
of coupled flavours. It is very important to note that QED and QCD are theories with 
"decoupling" : up to the scale Q only quarks with masses m « Q contribute to the running 
of a. This is clearly very important, given that all applications of perturbative QCD so 
far apply to energies below the top quark mass m t . For the validity of the decoupling 
theorem [59] it is necessary that the theory where all the heavy particle internal lines 
are eliminated is still renormalizable and that the coupling constants do not vary with 
the mass. These requirements are true for the mass of heavy quarks in QED and QCD, 
but are not true in the electroweak theory where the elimination of the top would violate 
SU (2) symmetry (because the t and b left-handed quarks are in a doublet) and the quark 
couplings to the Higgs multiplet (hence to the longitudinal gauge bosons) are proportional 
to the mass. In conclusion, in QED and QCD, quarks with m » Q do not contribute 
to rif in the coefficients of the relevant (3 function. The effects of heavy quarks are power 
suppressed and can be taken separately into account. For example, in e + e~ annihilation 
for 2m c < Q < 2rrib the relevant asymptotics is for rif = 4, while for 2rrib < Q < 2m t 
rif = 5. Going accross the b threshold the j3 function coefficients change, so the a(t) slope 
changes. But a(t) is continuous, so that A changes so as to keep a(t) constant at the 
matching point at Q ~ o(2m&). The effect on A is large: approximately A 5 ~ 0.65A 4 
where A 4>5 are for rif = 4, 5. 

Note the presence of a pole in eqs. (139,140) at ±bat = 1, called the Landau pole, who 
realised its existence in QED already in the '50's. For \x ~ m e (in QED) the pole occurs 
beyond the Planck mass. In QCD the Landau pole is located for negative t or at Q < n 
in the region of light hadron masses. Clearly the issue of the definition and the behaviour 
of the physical coupling (which is always finite, when defined in terms of some physical 
process) in the region around the perturbative Landau pole is a problem that lies outside 
the domain of perturbative QCD. 

The non leading terms in the asymptotic behaviour of the running coupling can in 
principle be evaluated going back to eq. (129) and computing b' at 2-loops and so on. But 
in general the perturbative coefficients of (3(a) depend on the definition of the renormal- 
ized coupling a (the renormalization scheme), so one wonders whether it is worthwhile 
to do a complicated calculation to get b' if then it must be repeated for a different def- 
inition or scheme. In this respect it is interesting to remark that actually both b and 
b' are independent of the definition of a, while higher order coefficients do depend on 
that. Here is the simple proof. Two different perturbative definitions of a are related by 
a' ~ a(l + C\Ol + ...). Then we have: 



f3(a>) 



da' 



da 



(1 + 2ci« + ...) 



d log /i 2 



d log /i 2 
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= p(a)(l + 2 Cl a + ...) 

= ±ba 2 {l + b'a + ...)(1 + 2ci« + ...) 

= ±ba' 2 {l + b'a' + ...) (142) 

which shows that, up to the first subleading order, /3(a') has the same form as (3(a). 
In QCD (Nc = 3) one has calculated: 

b' = 153 ~ 19 "/ (143 ) 
2vr(33 - 2n f ) v ; 

By taking b' into account one can write the expression of the running coupling at next to 
the leading order (NLO): 

a(Q 2 ) = a LO (Q 2 )[l - b'a LO (Q 2 )log\og^ + ...] (144) 

where a^Q = blogQ 2 / A 2 is the LO result (actually at NLO the definition of A is modified 
according to b log fi 2 /A 2 = 1/a + b' log ba). 

Summarizing, we started from massless classical QCD which is scale invariant. But we 
have seen that the procedure of quantization, regularization and renormalization necessarily 
breaks scale invariance. In the quantum QCD theory there is a scale of energy, A, which 
from experiment is of the order of a few hundred MeV, its precise value depending on 
the definition, as we shall see in detail. Dimensionless quantities depend on the energy 
scale through the running coupling which is a logarithmic function of Q 2 /A 2 . In QCD 
the running coupling decreases logarithmically at large Q 2 (asymptotic freedom), while in 
QED the coupling has the opposite behaviour. 



2.5 More on the Running Coupling 

In the previous section we have introduced the renormalized coupling a in terms of the 
3-gluon vertex at p 2 = —fi 2 (momentum subtraction). The Ward identities of QCD then 
ensure that the coupling defined from other vertices like the qqg vertex are renormalized in 
the same way and the finite radiative corrections are related. But at present the universally 
adopted definition of a s is in terms of dimensional regularization [60], because of computa- 
tional simplicity, which is essential given the great complexity of present day calculations. 
So we now briefly review the principles of dimensional regularization and the definition of 
Minimal Subtraction (MS) [61] and Modified Minimal Subtraction (MS) [62]. The M~S 
definition of a s is the one most commonly adopted in the literature and a value quoted for 
it is normally referring to this definition. 

Dimensional Regularization (DR) is a gauge and Lorentz invariant regularization that 
consists in formulating the theory in D < 4 spacetime dimensions in order to make loop 
integrals ultraviolet finite. In DR one rewrites the theory in D dimensions (D is integer 
at the beginning, but then one realizes that the calculated expression of diagrams makes 
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sense at all D except for isolated singularities) . The metric tensor is extended into a D x D 
matrix g^ u = diag(l, — 1, —1, —1) and 4-vectors are given by = (k°, k 1 , k ^ 1 ). 
The Dirac 7 M are f(D) x f(D) matrices and the precise form of the function f(D) is 
not important. It is sufficient to extend the usual algebra in a straightforward way like 
{it^lv} = 29fM,uI, with I the D-dimensional identity matrix, 7 M 7^7 A1 = — (D — 2)^ v or 
Tr(YY) = f(D)g^. 

The physical dimensions of fields change in D dimensions and, as a consequence, the 
gauge couplings become dimensional cd = /i e e, where e is dimensionless, Z) = 4 — 2e and 
fi is a scale of mass (this is how a scale of mass is introduced in the DR of massless QCD!). 
In fact, the dimension of fields is determined by requiring that the action S = J dPxL is 
dimensionless. By inserting for L terms like m\i/\I/ or m 2 <ft(f) or e^^^A^ the dimensions 
of the fields and coupling are determined as: m, ^ ',<p, A^e = l,(D — l)/2,(D — 2)/2,(D — 
2)/2, (4 — D)/2, respectively. The formal expression of loop integrals can be written for 
any D. For example: 

d D k 1 _ T(2- J D/2)(-m 2 ) D / 2 - 2 ^ 



(2n) D (k 2 - m 2 ) 2 (4tt)^/ 2 
For D = 4 — 2e one can expand using: 

r(e) = - - 7£ + o(e), 1e = 0.5772 (146) 

For some Green function G, normalized to 1 in lowest order, (like V/e with V the 3-g 
vertex function at the symmetric point p 2 = q 2 = r 2 , considered in the previous section) 
we typically find at 1-loop: 

-u 2 1 

G bare = 1 + a (-^-) e [S(-+log47r-7E) + A + o(e)] (147) 



In MS one rewrites this as (diagram by diagram: this is a virtue of the method): 

Gfyare ZG ren 

Z = 1 + a [5(- + log47r-7 S )] 

-a 2 

G ren = 1 + a [5 log -^- + A] (148) 

Here Z stands for the relevant product of renormalization factors. In the original MS 
prescription only 1/e was subtracted (that clearly plays the role of a cutoff) and not also 
log47r and jg. Later, since these constants always appear from the expansion of T functions 
it was decided to modify MS into MS. Note that the MS definition of a is different than 
that in the momentum subtraction scheme because the finite terms (those beyond logs) are 
different. In particular here the order a correction to G ren does not vanish at p 2 = — fi 2 . 

The third [63] and fourth [64] coefficients of the QCD f} function are also known in 
the MS prescription (recall that only the first two coefficients are scheme independent). 
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The calculation of the last term involved the evaluation of some 50,000 4-loop diagrams. 
Translated in numbers, for rif = 5 one obtains : 

f3(a) = -0.610a 2 [l + 1.261...- + 1.475. ..(-) 2 + 9.836.. ,(-) 3 ...] (149) 

7T 7T 7T 

It is interesting to remark that the expansion coefficients are of order 1 or 10 (only for the 
last one), so that the MS expansion looks reasonably well behaved. 

2.6 On the Non-convergence of Perturbative Expansions 

It is important to keep in mind that the QED and QCD perturbative series, after 
renormalization, have all their coefficients finite, but the expansion does not converge. 
Actually the perturbative series is not even Borel summable (for reviews, see, for example 
refs. [65]). After the Borel resummation, for a given process one is left with a result which is 
ambiguous by terms typically down by exp —n/(ba), with n an integer and b is the absolute 
value of the first /3 function coefficient. In QED these corrective terms are extremely small 
and not very important in practice. On the contrary in QCD a = a s (Q 2 ) ~ 1/(6 log Q 2 /A 2 ) 
and the ambiguous terms are of order (1/Q 2 ) n , that is are power suppressed. It is interesting 
that, through this mechanism, the perturbative version of the theory is able to somehow 
take into account the power suppressed corrections. A sequence of diagrams with factorial 
growth at large order n is made up by dressing gluon propagators by any number of 
quark bubbles together with their gauge completions (renormalons).The problem of the 
precise relation between the ambiguities of the perturbative expansion and the higher 
twist corrections has been discussed in recent years [65]. 
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2.7 e + e Annihilation and Related Processes 
2.7.1 R e+e - 

The simplest hard process is R = R e + e - = a(e + e~ — > hadrons)/a po i nt (e + e~ — > /i + yU~) 
that we have already introduced. R is dimensionless and in perturbation theory is given 
by R = Nc^iQiFtyjQZs), where F = 1 + o(a s ). 2 We have already mentioned 
that for this process the "anomalous dimension" function vanishes: 7(a s ) = because of 
electric charge non renormalization by strong interactions. Let us recall how this happens 
in detail. The diagrams that are relevant for charge renormalization in QED at 1-loop are 
shown in Fig. 14. The Ward identity that follows from gauge invariance in QED imposes 
that the vertex (Zy) and the self-energy (Zf) renormalization factors cancel and the only 
divergence remains in Z 1 , the vacuum polarization of the photon. So the charge is only 
renormalized by the photon vacuum polarization blob, hence it is universal (the same factor 
for all fermions, independent of their charge) and is not affected by QCD at 1-loop. It is 
true that at higher orders the photon vacuum polarization diagram is affected by QCD 
(for example, at 2-loops we can exchange a gluon between the quarks in the loop) but 
the renormalization induced by the divergent logs from the vacuum polarization diagram 
remain independent of the nature of the fermion to which the photon line is attached. The 
gluon contributions to the vertex (Zy) and to the self-energy (Zf) cancel because they have 
exactly the same structure as in QED, and there is no gluon contribution to the photon 
blob at 1-loop, so that j(a s ) = 0. 





Figure 14: Diagrams for charge renormalization in QED at 1-loop (the blob, in each diagram, represents 
the loop). 



At 1-loop the diagrams relevant for the computation of R are shown in Fig. 15. There 
are virtual diagrams and also real diagrams with one additional gluon in the final state. 
Infrared divergences cancel between the interference term of the virtual diagrams and the 
absolute square of the real diagrams, according to the Bloch-Nordsieck theorem. Similarly 

2 Actually starting from the order a 2 there are some "singlet" terms proportional to [^,- Qi] 2 . These 
small terms are included in F by dividing and multiplying by J2 i Q 2 . 
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there are no mass singularities, in agreement with the Kinoshita-Lee-Nauenberg theorem, 
because the initial state is purely leptonic and all degenerate states that can appear at the 
given order are included in the final state. Given that j(cx s ) = the RGE prediction is 
simply given, as we have already seen, by F(t, a s ) = F[0, a s (t)}. This means that if we do, 
for example, a 2-loop calculation, we must obtain a result of the form: 

F(t,a s ) = 1 + cia s (l-ba s t) + c 2 a 2 s + o(a 3 s ) (150) 

In fact, taking into account the expression of the running coupling in eq. (140): 

(X 

a «(*) ~ rrr 7 ~ a *Q ~ bast + ••••) ( 151 ) 

1 + ba s t 

eq. 150 can be rewritten as 

F(t,a s ) = 1 + Cl a 9 (t) + c 2 a 2 (t) + o(a 3 s (t)) = F[0,a s (t)} (152) 

The content of the RGE prediction is, at this order, that there are no a s t and (a s t) 2 terms 
(the leading log sequence must be absent) and the term of order a 2 s t has the appropriate 
coefficient to be reabsorbed in the transformation of a s into a s (t). 




2 




Figure 15: Real and virtual diagrams relevant for the computation of R at 1-loop accuracy (the initial 
e + e~ has been omitted to make the drawing simpler). 

At present the first 4 coefficients c\, ...,04 have been computed in the MS scheme (the 
references are: for c 2 [66], for c 3 [67] and for c 4 [68]). Clearly c\ = l/n does not depend 
on the definition of ot s but c n with n > 2 do. The subleading coefficients also depend 
on the scale choice: if instead of expanding in a s (Q) we decide to choose a s (Q/2) the 
coefficients c n n > 2 change. In the MS scheme, for 7-exchange and n/ = 5, which are 
good approximations for 2m b « Q « m z , one has: 

F[0,a s {t)\ = 1 + ^1 + i.409...(^) 2 - 12.8.... (^) 3 - 80.0... (^) 4 + ... (153) 

7T TT TT TT 

Similar perturbative results at 3-loop accuracy also exist for Rz = T(Z — > hadrons) /Y(Z — > 
leptons), R T = T(r — > v T + hadrons) /Y{r — > Ur+leptons), etc. We will discuss these results 
in Sect. 2.10 where we deal with measurements of a s . 
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The perturbative expansion in powers of a s (t) takes into account all contributions that 
are suppressed by powers of logarithms of the large scale Q 2 ("leading twist" terms). In 
addition there are corrections suppressed by powers of the large scale Q 2 ("higher twist" 
terms). The pattern of power corrections is controlled by the light-cone Operator Product 
Expansion (OPE) [69, 70] which (schematically) leads to: 

m 2 < OlTrfF^F^llO > < 0|O 6 |0 > /1K „, 
F = pert. + r 2 — + r 4 ' L ^ Jl + ... + r 6 ^ + ... (154) 

Here m 2 generically indicates mass corrections, for example from b quarks, beyond the 
b threshold, while top quark mass corrections only arise from loops, vanish in the limit 
m t — > oo and are included in the coefficients as those in eq. (153) and the analogous ones 
for higher twist terms; F^ u = Xm-^uI^j ^6 is typically a 4-fermion operator, etc. For 
each possible gauge invariant operator the corresponding negative power of Q 2 is fixed by 
dimensions. 

We now consider the light-cone OPE in some more detail. R e + e ~ ~ Tl(Q 2 ) where Tl(Q 2 ) 
is the scalar spectral function related to the hadronic contribution to the imaginary part 
of the photon vacuum polarization T^ v \ 

T„, = {-g^Q 2 + q^MQ 2 ) = J d 4 xex P i(q-x) < 0| jt(x)J„(0)|0 > = 

= 5^<0|jt(0)|n><n|J,(0)|0>(27r) 4 5 4 (g-p n ) (155) 



For Q 2 — > oo the x 2 — > region is dominant. The light cone OPE is valid to all orders in 
perturbation theory. Schematically, dropping Lorentz indices, for simplicity, near x 2 ~ 
we have: 

oo 

j\x)J(0) = I(x 2 ) + E(x 2 ) ^c n (x 2 )x^...x^ ■ . Mn (0) +less sing, terms (156) 

n=0 

Here I(x 2 ), E(x 2 ),..., c n (x 2 ) are c-number singular functions, O n is a string of local op- 
erators. E(x 2 ) is the singularity of free field theory, I(x 2 ) and c n (x 2 ) in the interacting 
theory contain powers of log (/i 2 x 2 ). Some O n are already present in free field theory, other 
ones appear when interactions are switched on. Given that Tl(Q 2 ) is related to the Fourier 
transform of the vacuum expectation value of the product of currents, less singular terms 
in x 2 lead to power suppressed terms in l/Q 2 . The perturbative terms, like those in eq. 
152, come from I(x 2 ) which is the leading twist termand the dominant logarithmic scaling 
violations induced by the running coupling are the logs in I(x 2 ). 

2.7.2 The Final State in e + e Annihilation 

Experiments on e + e~ annihilation at high energy provide a remarkable possibility of 
systematically testing the distinct signatures predicted by QCD for the structure of the 
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final state averaged over a large number of events. Typical of asymptotic freedom is the 
hierarchy of configurations emerging as a consequence of the smallness of a s (Q 2 ). When all 
corrections of order a s (Q 2 ) are neglected one recovers the naive parton model prediction for 
the final state: almost collinear events with two back-to-back jets with limited transverse 
momentum and an angular distribution as (1+cos 2 9) with respect to the beam axis (typical 
of spin 1/2 parton quarks: scalar quarks would lead to a sin 2 # distribution). At order 
a s (Q 2 ) a tail of events is predicted to appear with large transverse momentum px ~ <5/2 
with respect to a suitably defined jet axis (for example the thrust axis, see below). This 
small fraction of events with large pt mostly consists of three-jet events with an almost 
planar topology. The skeleton of a three-jet event, at leading order in a s (Q 2 ), is formed 
by three hard partons qqg, the third being a gluon emitted by a quark or antiquark line. 
At order a 2 (Q 2 ) a hard perturbative non planar component starts to build up and a small 
fraction of four-jet events qqgg or qqqq appear, and so on. 

Event shape variables defined from the set of 4-momenta of final state particles are 
introduced to quantitatively describe the topological structure of the final state energy 
flow [71]. The most well known event shape variable is thrust (T) [72] defined as: 

T = max ^iJfVJT 1 ( 157 ) 

Li \Pi\ 

where the maximization is in terms of the axis defined by the unit vector n^: the thrust 
axis is the axis that maximizes the sum of the absolute values of the longitudinal momenta 
of the final state particles. The thrust T varies between 1/2, for a spherical event, to 
1 for a collinear (2-jet) event. Event shape variables are important for QCD tests and 
measurements of a s and also for more practical purposes like a laboratory for assessing 
the reliability of event simulation programmes and a tool for the separation of signals and 
background. 

A quantitatively specified definition of jets and of the number of jets in one event (jet 
counting) must be introduced for precise QCD tests and for measuring a s , which must be 
infrared safe (i.e. not altered by soft particle emission or collinear splittings of massless 
particles) in order to be computable at parton level and as much as possible insensitive to 
the transformation of partons into hadrons (see, for example, ref. [73]). For e + e~ physics 
one has used a jet algorithm based on a resolution parameter y cut and a suitable pair 
variable; for example [74]: 

2mm(E 2 ,Ej)(l-cos%) 

Vij = (158) 

s 

Note that 1 — cos 8 ij ~ 0f,-/2 so that the relative transverse momentum k\ is involved (hence 
the name kr algorithm). The particles i,j belong to different jets for y^ > y cut . Clearly the 
number of jets becomes a function of y cu t- there are more jets for smaller y cut . 

Recently, motivated by the LHC experiments there has been a flurry of improved jet 
algorithm studies: it is essential that a correct jet finding is implemented by LHC ex- 
periments for an optimal matching of theory and experiment [75,76]. In particular the 
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existing sequential recombination algorithms like [74], [77] and Cambridge/Aachen [78] 
have been generalized. In this recursive definitions one introduces distances dij between 
particles or clusters of particles i and j and diB between i and the beam (B). The inclu- 
sive clustering proceeds by identifying the smallest of the distances and, if it is a d^, by 
recombining particles i and j, while, if it is diB, calling i a jet and removing it from the 
list. The distances are recalculated and the procedure repeated until no i and j are left. 
The extension relative to the kx [77] and Cambridge/ Aachen [78] algorithms lies in the 
definition of the distance measures: 

2 2 A r 

= minik^k^)^-, (159) 

where A?- = (yi — yj) 2 + ((pi — (fij) 2 and kn, Hi and 4>i are respectively the transverse 
momentum, rapidity and azimuth of particle i. R is the radius of the jet, i.e. the radius 
of a cone that, by definition, contains the jet. The exponent p fixes the relative power of 
the energy versus geometrical (Ay) scales. For p = 1 one has the inclusive kx algorithm. 
It can be shown in general that for p > the behaviour of the jet algorithm with respect 
to soft radiation is rather similar to that observed for the kx algorithm. The case p = is 
special and it corresponds to the inclusive Cambridge/ Aachen algorithm [78]. At first sight 
surprisingly, taking p to be negative also yields an algorithm that is infrared and collinear 
safe and has sensible phenomenological behaviour. For p = — 1 one obtains the recently 
introduced anti-A^ jet-clustering algorithm [79] which has particularly stable jet boundaries 
with respect to soft radiation and is suitable for practical use by the experiments. 

2.8 Deep Inelastic Scattering 

Deep Inelastic Scattering (DIS) processes have played and still play a very important 
role for our understanding of QCD and of nucleon structure. This set of processes actually 
provides us with a rich laboratory for theory and experiment. There are several structure 
functions that can be studied, Fi(x,Q 2 ), each a function of two variables. This is true 
separately for different beams and targets and different polarizations. Depending on the 
charges of I and V (see eq. (106)) we can have neutral currents (7,Z) or charged currents in 
the l-V channel (Fig. 12). In the past DIS processes were crucial for establishing QCD as 
the theory of strong interactions and quarks and gluons as the QCD partons. At present 
DIS remains very important for quantitative studies and tests of QCD. The theory of 
scaling violations for totally inclusive DIS structure functions, based on operator expansion 
or diagrammatic techniques and renormalization group methods, is crystal clear and the 
predicted Q 2 dependence can be tested at each value of x. The measurement of quark 
and gluon densities in the nucleon, as functions of x at some reference value of Q 2 , which 
is an essential starting point for the calculation of all relevant hadronic hard processes, is 
performed in DIS processes. At the same time one measures a s (Q 2 ) and the DIS values 
of the running coupling can be compared with those obtained from other processes. At all 
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times new theoretical challenges arise from the study of DIS processes. Recent examples 
(see the following) are the so-called "spin crisis" in polarized DIS and the behaviour of 
singlet structure functions at small x as revealed by HERA data. In the following we will 
review the past successes and the present open problems in the physics of DIS. 

The cross-section a ~ L^W^ is given in terms of the product of a leptonic (L^) and 
a hadronic [W^ u ) tensor. While L^ u is simple and easily obtained from the lowest order 
electroweak (EW) vertex plus QED radiative corrections, the complicated strong interac- 
tion dynamics is contained in W^ u . The latter is proportional to the Fourier transform of 
the forward matrix element between the nucleon target states of the product of two EW 
currents: 

= J d A xexpi(q-x) <p\Jl(x)J u (0)\p> (160) 

Structure functions are defined starting from the general form of W^ v given Lorentz in- 
variance and current conservation. For example, for EW currents between unpolarized 
nucleons we have: 

= (- 9flv + -^-) W 1 (is,Q ) + (p M - -j-qjipv - -j-qv) ~ 2 

e^x P p X q p W 3 (v,Q 2 ) 



2m 2 



W 3 arises from VA interference and is absent for pure vector currents. In the limit Q 2 » 
m 2 , x fixed, the structure functions obey approximate Bjorken scaling which in reality is 
broken by logarithmic corrections that can be computed in QCD: 

mW^Q 2 ) ->■ F 1 (x) 

vW 2t3 (is,Q 2 ) F 2)3 (x) (161) 
The 7 — N cross-section is given by (Wi = Wi(Q 2 , u)): 

[2sm 2 -W x + cos 2 -W 2 ] (162) 



dQ 2 dv Q A E 1 2 2 

while for the h> — N ox v — N cross-section one has: 

d<r u,v GpE' mly 2 . o 2®-<it 1 E + E' . J , 

-— — - = — - — (—. — ! ^-) • 2 sm 2 -W 1 + cos 2 -W 2 ± sm 2 -W 3 \ (163) 

dQ 2 dv 2nE K Q 2 + m 2 w J 1 2 1 2 m 2 Sl K ' 

(Wi for photons, v and v are all different, as we shall see in a moment). 

In the scaling limit the longitudinal and transverse cross sections are given by: 

a L ~ -[— - Fl (x)} 
&rh,lh ~ -[Fi(x) ± F 3 (x)] 

°T = &RH + &LH (164) 
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where L, RH, LH refer to the helicity 0, 1, -1, respectively, of the exchanged gauge vector 
boson. For the photon case F 3 = and a RH = a LH . 

In the '60's the demise of hadrons from the status of fundamental particles to that of 
bound states of constituent quarks was the breakthrough that made possible the construc- 
tion of a renormalizable field theory for strong interactions. The presence of an unlimited 
number of hadrons species, many of them with large spin values, presented an obvious dead- 
end for a manageable field theory. The evidence for constituent quarks emerged clearly 
from the systematics of hadron spectroscopy. The complications of the hadron spectrum 
could be explained in terms of the quantum numbers of spin 1/2, fractionally charged, u, d 
and s quarks. The notion of colour was introduced to reconcile the observed spectrum with 
Fermi statistics. But confinement that forbids the observation of free quarks was a clear 
obstacle towards the acceptance of quarks as real constituents and not just as fictitious 
entities describing some mathematical pattern (a doubt expressed even by Gell-Mann at 
the time). The early measurements at SLAC of DIS dissipated all doubts: the observation 
of Bjorken scaling and the success of the "naive" (not so much after all) parton model of 
Feynman imposed quarks as the basic fields for describing the nucleon structure (parton 
quarks) . 

In the language of Bjorken and Feynman the virtual 7 (or, in general, any gauge 
boson) sees the quark partons inside the nucleon target as quasi-free, because their (Lorentz 
dilated) QCD interaction time is much longer than r 7 ~ 1/Q, the duration of the virtual 
photon interaction. Since the virtual photon 4-momentum is spacelike, we can go to a 
Lorentz frame where _E 7 = (Breit frame). In this frame q = {E 1 = 0; 0, 0, Q) and the 
nucleon momentum, neglecting the mass m « Q, is p = (Q/2x; 0,0, —Q/2x) (note that 
this correctly gives q 2 = —Q 2 and x = Q 2 /2(p ■ q)). Consider (Fig. 16) the interaction of 
the photon with a quark carrying a fraction y of the nucleon 4-momentum: p q = yp (we are 
neglecting the transverse components of p q which are of order m). The incoming parton 
with p q = yp absorbs the photon and the final parton has 4-momentum p' q . Since in the 
Breit frame the photon carries no energy but only a longitudinal momentum Q, the photon 
can only be absorbed by those partons with y = x: then the longitudinal component of 
p q = yp is —yQ/2x = —Q/2 and can be flipped into +Q/2 by the photon. As a result, the 
photon longitudinal momentum +Q disappears, the parton quark momentum changes of 
sign from —Q/2 into +Q/2 and the energy is not changed. So the structure functions are 
proportional to the density of partons with fraction x of the nucleon momentum, weighted 
with the squared charge. Also, recall that the helicity of a massless quark is conserved in a 
vector (or axial vector) interaction (recall Sect. 1.5). So when the momentum is reversed 
also the spin must flip. Since the process is collinear there is no orbital contribution and 
only a photon with helicity ±1 (transverse photon) can be absorbed. Alternatively, if 
partons were spin zero only longitudinal photons would instead contribute. 

Using these results, which are maintained in QCD at leading order, the quantum num- 
bers of the quarks were confirmed by early experiments. The observation that R = olIo t ^- 
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Figure 16: Schematic diagram for the interaction of the virtual photon with a parton quark in the Breit 
frame. 



implies that the charged partons have spin 1/2. The quark charges were derived from 
the data on the electron and neutrino structure functions: 

4 1 4 1 

F ep = -u(x) + -d(x) + ; F en = -d(x) + -u(x) + .... 

y y y y 

F vp = F Pn = 2d(x) + ; F un = F Pp = 2u(x) + (165) 

where F ~ 2F\ ~ F2/X and u(x), d(x) are the parton number densities in the proton 
(with fraction x of the proton longitudinal momentum), which, in the scaling limit, do not 
depend on Q 2 . The normalization of the structure functions and the parton densities are 
such that the charge relations hold: 

/ [u(x) — u(x)]dx = 2, / [d(x) — d(x)]dx = 1, / [s(x) — s(x)]dx = (166) 
Jo Jo Jo 

Also it was proven by experiment that at values of Q 2 of a few GeV 2 , in the scaling region, 
about half of the nucleon momentum, given by the momentum sum rule: 



(qi(x) + qi(x)) + g(x)]xdx = 1 (167) 



is carried by neutral partons (gluons). 

In QCD there are calculable log scaling violations induced by a s (t). The parton rules 
in eq.165 can be summarized in the schematic formula: 

F(x,t) = [ dy^a point (-,a s (t)) + o(-U (168) 

J x y y ^ 

Before QCD corrections a po i nt = e 2 5(x/y — 1) and F = e 2 qo(x) (here we denote by e the 
charge of the quark in units of the positron charge, i.e. e = 2/3 for the u quark). QCD 
modifies a point at order a s via the diagrams of Fig. 17. From a direct computation of the 
diagrams one obtains a result of the following form: 

a point (z,a s (t)) ~ e 2 [5(z-l) + ^{t ■ P(z) + f{z))\ (169) 

Note that the y integral in eq. 168 is from x to 1, because the energy can only be lost 
by radiation before interacting with the photon (which eventually wants to find a fraction 
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g 2 



(a) (b) (c) (d) 

Figure 17: First order QCD corrections to the virtual photon-quark cross-section. 



x, as we have explained). For y > x the correction arises from diagrams with real gluon 
emission. Only the sum of the two real diagrams in Fig. 17 is gauge invariant, so that 
the contribution of one given diagram is gauge dependent. But in an axial gauge, which 
for this reason is some times also called the "physical gauge", the diagram of Fig. 17(c), 
among real diagrams, gives the whole i-proportional term at < x < 1. It is obviously not 
essential to go to this gauge, but this diagram has a direct physical interpretation: a quark 
in the proton has a fraction y > x of the parent 4-momentum; it then radiates a gluon and 
looses energy down to a fraction x before interacting with the photon. The log arises from 
the virtual quark propagator, according to the discussion of collinear mass singularities in 
eq. (105). In fact in the massless limit one has: 



propagator 



1 



1 



-1 



(k-h) 2 
1 



2E k E h 
-1 



1 — cos 9 



AE k E h sin 2 9/2 



oc 



Pt 



(170) 



where px is the transverse momentum of the virtual quark. So the square of the propagator 
goes like l/p?- But there is a p\ factor in the numerator, because in the collinear limit, 
when 9 = and the initial and final quarks and the emitted gluon are all aligned, the 
quark helicity cannot flip (vector interaction) so that the gluon should carry helicity zero 
while a real gluon can only have ±1 helicity. Thus the numerator vanishes as p\ in the 
forward direction and the cross-section behaves as: 



a 



\dp 2 T 
Pt 



logQ 2 



(171) 



Actually the log should be read as log Q 2 /m 2 because in the massless limit a genuine mass 
singularity appears. In fact the mass singularity connected with the initial quark line is 
not cancelled because we do not have the sum of all degenerate initial states [57], but 
only a single quark. But in correspondence to the initial quark we have the (bare) quark 
density qo(y) that appears in the convolution integral. This is a non perturbative quantity 
determined by the nucleon wave function. So we can factorize the mass singularity in a 
redefinition of the quark density: we replace go (y) q(y, t) = qo(y) + Aq(y,t) with: 



Aq(x,t) 



2^ Jx y 



p 



X , 



(172) 
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Here the factor of £ is a bit symbolic: it stands for \ogQ 2 /m 2 but what we exactly put 
under Q 2 depends on the definition of the renormalized quark density, which also fixes the 
exact form of the finite term f(z) in eq. (169). 

The effective parton density q(y,t) that we have defined is now scale dependent. In 
terms of this scale dependent density we have the following relations, where we have also 
replaced the fixed coupling with the running coupling according to the prescription derived 
from the RGE: 

F(x,t) = f 1 dy q -^e 2 [5( X - - 1) + °^-f(-))} = e 2 q(x,t) + o(a s (t)) 

Jx y y 2f 1/ 

d f a s (t) f 1 q(y,t) x 2 

-7T90M) = -s — / d y p " + o(a 8 (t) 173 

dt 2tt J x y y 



We see that in lowest order we reproduce the naive parton model formulae for the structure 
functions in terms of effective parton densities that are scale dependent. The evolution 
equations for the parton densities are written down in terms of kernels (the "splitting 
functions" [80]) that can be expanded in powers of the running coupling. At leading order, 
we can interpret the evolution equation by saying that the variation of the quark density 
at x is given by the convolution of the quark density at y times the probability of emitting 
a gluon with fraction x/y of the quark momentum. 

It is interesting that the integro-differential QCD evolution equation for densities can 
be transformed into an infinite set of ordinary differential equations for Mellin moments [5] . 
The Mellin moment f n of a density f(x) is defined as: 

/„ = I dxx n - x f{x) (174) 
Jo 

By taking moments of both sides of the second of eqs.(173) one finds, with a simple 
interchange of the integration order, the simpler equation for the n-th moment: 

j t q n {t) = °^--Pn-q n (t) (175) 

To solve this equation we observe that: 

, Quit) P n f* , s , p n f as{t) da' , , 

log^-( = f - / a s (tdt = f - / — 176 
q n (0) 2tt J 2tt J as -ba' 

where we used eq. (124) to change the integration variable from dt to da{t) (denoted as 
da') and (3(a) ~ —ba 2 + .... Finally the solution is: 

Init) = [-%]^-q n (0) (1T7) 
a s (t) 

The connection of these results with the RGE general formalism occurs via the light 
cone OPE (recall eq. (160) for W^ v and eq. (156) for the OPE of two currents). In the 
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case of DIS the c-number term I(x 2 ) does not contribute, because we are interested in the 
connected part of the matrix element < p|...|p > — < 0|...|0 >. The relevant terms are: 

oo 

jt( x )j(o) = Eix^^c^x^x^.-.x^ -O™ Mn (0) + less sing, terms (178) 

n=0 

A formally intricate but conceptually simple argument based on the analiticity properties 
of the forward virtual Compton amplitude shows that the Mellin moments M n of structure 
functions are related to the individual terms in the OPE, precisely to the Fourier transform 
c n (Q 2 ) (we will write it as c n (t, a)) of the coefficient c n (x 2 ) times a reduced matrix element 
h n from the operators O n : < p\0^ ... Mn (0)|p >= Kp^.-.p^: 



•1 



c n < p\O n \p >->■ M n = dxx n - L F(x) (179) 

Jo 

Since the matrix element of the products of currents satisfy the RGE so do the moments 
M n . Hence the general form of the Q 2 dependence is given by the RGE solution (see 
eq. (126)): 

ra(t) 7nK) 



M n (t,a) = c„[0,a(t)]exp / ^—^da'-h n {a) (180) 
In lowest order, identifying in the simplest case M n with q n , we have: 

ln {a) = + P{a) = - ba 2 + ... (181) 

and 

q n (t) = g n (0)exp / ?%-±dct = [— S -A^ ■ q n (0) (182) 

Ja Pipe) a 8 (t) 

which exactly coincides with eq. (177). 

Up to this point we have implicitly restricted our attention to non-singlet (under the 
flavour group) structure functions. The Q 2 evolution equations become non diagonal as 
soon as we take into account the presence of gluons in the target. In fact the quark which 
is seen by the photon can be generated by a gluon in the target (Fig. 18). 

The quark evolution equation becomes: 

j f QiM = °^-fa®P qq ] + ^l9®P g9 ] (183) 
where we introduced the shorthand notation: 

[q®p] = [P®q] = I dy^^--P{-) (184) 

j x y y 

(it is easy to check that the so-defined convolution, like an ordinary product, is commuta- 
tive). At leading order, the interpretation of eq. (183) is simply that the variation of the 
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Figure 18: Lowest order diagram for the interaction of the virtual photon with a parton gluon. 

quark density is due to the convolution of the quark density at a higher energy times the 
probability of finding a quark in a quark (with the right energy fraction) plus the gluon 
density at a higher energy times the probability of finding a quark (of the given flavour i) 
in a gluon. The evolution equation for the gluon density, needed to close the system 3 , can 
be obtained by suitably extending the same line of reasoning to a gedanken probe sensitive 
to colour charges, for example a virtual gluon. The resulting equation is of the form: 

j t 9&t) = °^-[J2iQi + Qi)®P 9q ] + °^[9®P 99 \ (185) 

i 

The explicit form of the splitting functions in lowest order [80,81,83] can be directly 
derived from the QCD vertices [80] . They are a property of the theory and do not depend 
on the particular process the parton density is taking part in. The results are : 

4 1 + x 2 3 

Pqq = 3 [ (l-x) + + 2 6{1 ~ X)] + ° {as) 

41 + (l-x) 2 
P gq = K — '- + o(a s ) 

Pgg = \[X 2 + {1~X) 2 } + 0(a s ) 



3 Thc evolution equations are now often called DGLAP equations (Dokshitzer Gribov Lipatov Altarelli 
Parisi). The first article by Gribov and Lipatov was published in 72 [81] (even before the works by Gross 
and Wilczek and by Politzer!) and was followed in 74 by a paper by Lipatov [82] (these dates correspond 
to the publication in russian) . All these articles refer to an abelian vector theory (treated in parallel with a 
pseudoscalar theory). Seen from the point of view of the evolution equations, these papers, in the context 
of the abelian theory, ask the right question and extract the relevant logarithmic terms from the dominant 
class of diagrams. But from their formal presentation the relation to real physics is somewhat hidden (in 
this respect the 74 paper by Lipatov makes some progress and explicitly refers to the parton model). The 
article by Dokshitser [83] was exactly contemporary to that by Altarelli Parisi [80]. It now refers to the 
non abelian theory (with running coupling) and the discussion is more complete and explicit than in the 
Gribov-Lipatov articles. But, for example, the connection to the parton model, the notion of the evolution 
as a branching process and the independence of the kernels from the process are not emphasized. 
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Figure 19: A representative selection ol data on the proton electromagnetic structure function from 
Collider (HERA) and fixed target experiments [14] that clearly shows the pattern of scaling violations. 
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For a generic non singular weight function f(x), the "+" distribution is defined as: 



The 5(1 — x) terms arise from the virtual corrections to the lowest order tree diagrams. 
Their coefficient can be simply obtained by imposing the validity of charge and momentum 
sum rules. In fact, from the request that the charge sum rules in eq. (166) are not affected 
by the Q 2 dependence one derives that 



which can be used to fix the coefficient of the 5(1 —x) terms of P qq . Similarly, by taking the 
t-derivative of the momentum sum rule in eq. (167) and imposing its vanishing for generic 
qi and g, one obtains: 



At higher orders the evolution equations are easily generalized but the calculation of 
the splitting functions rapidly becomes very complicated. For many years the splitting 
functions were only completely known at NLO accuracy [84]: a s P ~ a s Pi + °%p2 + •••• 
Then in recent years the NNLO results P3 have been first derived in analytic form for the 
first few moments and, then the full NNLO analytic calculation, a really monumental work, 
was completed in 2004 by Moch, Vermaseren and Vogt [85]. Beyond leading order a precise 
definition of parton densities should be specified. One can take a physical definition: for 
example, quark densities can be defined as to keep the LO expression for the structure 
function F 2 valid at all orders, the so called DIS definition [86], and the gluon density 
could be defined starting from Fl, the longitudinal structure function. Alternatively one 
can adopt a more abstract specification as, for example, in terms of the MS prescription. 
Once the definition of parton densities is fixed, the coefficients that relate the different 
structure functions to the parton densities at each fixed order can be computed. Similarly 
the higher order splitting functions also depend, to some extent, from the definition of 
parton densities, and a consistent set of coefficients and splitting functions must be used 
at each order. 

The scaling violations are clearly observed by experiment (Fig. 19) and their pattern 
is very well reproduced by QCD fits at NLO Fig. 20 [87]. These fits provide an impressive 
confirmation of a quantitative QCD prediction, a measurement of qi(x, Qi) and g(x ) Qf j ) 1 
at some reference value Ql of Q 2 , and a precise measurement of a s (Q 2 ). 
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;ure 20: A NLO QCD fit to the combined HERA data with Q 2 > 3.5 GeV 2 ): x 2 /dof = 574/582 [87]. 
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Figure 21: A more detailed view of the NLO QCD fit to a selection of the HERA data [87]. 
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Figure 22: The longitudinal structure function Fl measured by HI at HERA, as function of Q 2 for different 
values of x. The theoretical curves are obtained from different sets of parton densities as indicated. 

2.8.1 The Longitudinal Structure Function 

After SLAC established the dominance of the transverse cross section it took 40 years 
to get meaningful data on the longitudinal structure function F L (see eq.(164)! These data 
represent an experimental highlight of recent years. They have been obtained by HI at 
HERA [88]. The data are shown in Fig. (22). For spin 1/2 charged partons Fl vanishes 
asymptotically. In QCD Fl starts at order a s {Q 2 ). At LO the simple, 30 years old, formula 
is valid (for N f = 4) [89]: 



y 3 



\F 2 {y,Q*) + A \g{y,Q*){l- X -] 
3 9 v 



(190) 



The o(a 2 s ) [90] and o(af) [91] corrections are at present also known. One would not have 
expected that it would take such a long time to have a meaningful test of this simple 
prediction! And in fact better data would be highly desirable. But how and when they 
will be obtained is at present not clear at all. 

2.8.2 Large and Small x Resummations for Structure Functions 

At values of x either near or near 1 (with Q 2 large) those terms of higher order in a s 
in both the coefficients or the splitting functions which are multiplied by powers of log 1/x 
or log (1 — x) eventually become important and should be taken into account. Fortunately 
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the sequences of leading and subleading logs can be evaluated at all orders by special 
techniques and resummed to all orders. 

For x ~ 1 resummation [92] I refer to the recent papers [93], [94] (the latter also 
involving higher twist corrections, which are important at large x) where a list of references 
to previous work can be found. 

More important is the small x resummation because the singlet structure functions are 
large in this domain of x (while all structure functions vanish near x = 1). Here we will 
briefly summarize the small-x case for the singlet structure function which is the dominant 
channel at HERA, dominated by the sharp rise of the gluon and sea parton densities at 
small x. The small x data collected by HERA can be fitted reasonably well even at the 
smallest measured values of x by the NLO QCD evolution equations, so that there is no 
dramatic evidence in the data for departures. This is surprising also in view of the fact that 
the NNLO effects in the evolution have recently become available and are quite large [85] . 
Resummation effects have been shown to resolve this apparent paradox. For the singlet 
splitting function the coefficients of all LO and NLO corrections of order [a s (Q 2 ) logl/x] n 
and a s {Q 2 )[a s {Q 2 ) log l/x] n , respectively, are explicitly known from the Balitski, Fadin, 
Kuraev, Lipatov (BFKL) analysis of virtual gluon- virtual gluon scattering [95], [96]. But 
the simple addition of these higher order terms to the perturbative result (with subtraction 
of all double counting) does not lead to a converging expansion (the NLO logs completely 
overrule the LO logs in the relevant domain of x and Q 2 ). A sensible expansion is only 
obtained by a proper treatment of momentum conservation constraints, also using the 
underlying symmetry of the BFKL kernel under exchange of the two external gluons, 
and especially, of the running coupling effects (see the analysis in [97, 98] and references 
therein). In Fig. 23 we present the results for the dominant singlet splitting function 
xP gg (x, a s {Q 2 )) for a s (Q 2 ) ~ 0.2. We see that while the NNLO perturbative splitting 
function sharply deviates from the NLO approximation at small x, the resummed result 
only shows a moderate dip with respect to the NLO perturbative splitting function in the 
region of HERA data, and the full effect of the true small x asymptotics is only felt at 
much smaller values of x. The related effects are not very important for most processes at 
the LHC but could become relevant for next generation of hadron colliders. 

2.8.3 Polarized Deep Inelastic Scattering 

Polarized DIS is a subject where our knowledge is still far from satisfactory in spite of a 
great experimental effort (for recent reviews, see, for example, [99]). One main question is 
how the proton helicity is distributed among quarks, gluons and orbital angular momentum: 

l -AY, + Ag + L z = l - (191) 

Experiments with polarized leptons on polarized nucleons are sensitive to the polarized 
parton densities Aq = q + — q_, the difference of quark densities with helicity plus and 
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Figure 23: The dominant singlet splitting function xP gg (x 7 a s (Q 2 )) for a s (Q 2 ) ~ 0.2. The resummed 
results from ref. [97] (labeled ABF) and from ref. [98] (CCSS), which are in good agreement among them, 
are compared with the LO, NLO and NNLO perturbative results. 



minus in a proton with helicity plus. These differences are related to the quark matrix 
elements of the axial current. The polarized densities satisfy evolution equations analogous 
to eqs. 183,185 but with modified splitting functions that were derived in ref. [80] (the 
corresponding anomalous dimensions were obtained in ref. [100]). The measurements have 
shown that the quark moment AS is small (the "spin crisis" started by ref.( [101]): values 
from recent fits [102-106] are in the range AS ~ 0.2 — 0.3 (see Table 1): in any case, a 
less pronounced crisis than it used to be in the past. From the spin sum rule one obtains 
that either Ag + L z is relatively large or there are contributions to AS at very small x 
outside of the measured region. Denoting, for short hand, by Ag the first moment of the 
net helicity carried by the sum q + q we have the relations [103, 104]: 

a 3 = Au - Ad = (F + D)(l + e 2 ) = 1.269 ± 0.003 (192) 

a 8 = Aw + Ad - 2As = (3F - D)(l + e 3 ) = 0.586 ± 0.031 (193) 

where the F and D couplings are defined in the SU(3) flavour symmetry limit and e 2 and 
e 3 describe the SU(2) and SU(3) breakings, respectively. From the measured first moment 
of the structure function g x one obtains the value of a = AS: 



I', = / dxg x {x) = — 



a>3 + + 4a ) 



(194) 
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Table 1: The first moments of the 


singlet and gluon polarized PDFs at the scale Q 2 = 


AGeV 2 [102]. 


with the result, at Q 2 ~ 4GeV 2 








a = AS : 


= Au + Ad + As = a 8 + 3As 


~ 0.25 


(195) 


In turn, in the SU(3) limit e 2 = 


e 3 = 0, one then obtains: 






Aw ~ 0. 


82, Ad 0.45, As ~ - 


-0.11 


(196) 



This is an important result! Given F, D and Ti we know Au, Ad, As and AS in the 
SU(3) limit which should be reasonably accurate. The x distribution of g\ is known down 
to x ~ 10~ 4 on proton and deuterium and the 1st moment of g\ does not seem to get much 
from the unmeasured range at small x (also theoretically g\ should be smooth at small 
x [107]). The value of As ~ —0.11 from totally inclusive data and SU(3) appears to be 
at variance with the value extracted from single particle inclusive DIS (SIDIS) where one 
obtains a nearly vanishing result for As in a fit to all data [103] that leads to puzzling 
results. There is, in fact, an apparent tension between the 1st moments as determined by 
using the approximate SU (3) symmetry and from fitting the data on SIDIS (x > 0.001) (in 
particular for the strange density). But the adequacy of the SIDIS data is questionable (in 
particular of the kaon data which fix As) and of their theoretical treatment (for example, 
the application of parton results at too low an energy and the ambiguities on the kaon 
fragmentation function). 

AS is conserved in perturbation theory at LO (i.e. it does not evolve in Q 2 ). For 
conserved quantities we would expect that they are the same for constituent and for parton 
quarks. But actually the conservation of AS is broken by the axial anomaly and, in fact, in 
perturbation theory beyond LO the conserved density is actually AS' = AS + nf/2na s Ag 
[108]. Note that also a s Ag is conserved in LO, that is Ag ~ logQ 2 . This behaviour is 
not controversial but it will take long before the log growth of Ag will be confirmed by 
experiment! But by establishing this behaviour one would show that the extraction of 
Ag from the data is correct and that the QCD evolution works as expected. If Ag was 
large enough it could account for the difference between partons (AS) and constituents 
(AS'). From the spin sum rule it is clear that the log increase should cancel between 
Ag and L z . This cancelation is automatic as a consequence of helicity conservation in 
the basic QCD vertices. Ag can be measured indirectly by scaling violations and directly 
from asymmetries, e.g. in SIDIS. Existing measurements by HERMES, COMPASS, and 
at RHIC are still crude but show no hint of a large Ag at accessible values of x and Q 2 . 
Present data, affected by large errors (see Table 1) are consistent with a sizable contribution 
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of Ag to the spin sum rule in eq.(191) but there is no indication that a s Ag effects can 
explain the difference between constituents and parton quarks. 

2.9 Hadron Collider Processes and Factorization 

There are three classes of hard processes: those with no hadronic particles in the initial 
state, like e + e~ annihilation, those initiated by a lepton and a hadron, like DIS, and 
those with two incoming hadrons. The parton densities, defined and measured in DIS, 
are instrumental to compute hard processes initiated by collisions of two hadrons, like pp 
(Tevatron) or pp (LHC). Suppose you have a hadronic process of the form hi+h 2 — > X + all 
where hi are hadrons and X is some triggering particle or pair of particles or one or more 
jets which specify the large scale Q 2 relevant for the process, in general somewhat, but 
not much, smaller than s, the total c.o.m. squared mass. For example, X can be a W or 




Figure 24: Diagram for the Factorization Theorem 

a Z or a virtual photon with large Q 2 (Drell-Yan processes), or a jet at large transverse 
momentum p?, or a pair of heavy ( of mass M) quark- ant iquark. By "all" we mean a totally 
inclusive collection of hadronic particles. The Factorization Theorem (FT) states that for 
the total cross-section or some other sufficiently inclusive distribution we can write, apart 
from power suppressed corrections, the expression (see also Fig. 24): 

cr(s,r) = / dx 1 dx 2 piA(xi,Q 2 )p2B(x2,Q 2 )o- A B(xix 2 s,T) (197) 

AB 

Here r = Q 2 /s is a scaling variable, p iA are the densities for a parton of type A inside the 
hadron hi, gab is the partonic cross-section for parton-A + parton-B— > X + all'. Here 
alV is the partonic version of all: a totally inclusive collection of quarks, antiquarks and 
gluons. This result is based on the fact that the mass singularities that are associated with 
the initial legs are of universal nature, so that one can reproduce the same modified parton 
densities, by absorbing these singularities into the bare parton densities, as in DIS. Once 
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the parton densities and a s are known from other measurements, the prediction of the rate 
for a given hard process is obtained with not much ambiguity (e.g from scale dependence 
or hadronization effects). At least a NLO calculation of the reduced partonic cross-section 
gab is needed in order to correctly specify the scale and in general the definition of the 
parton densities and of the running coupling in the leading term. The residual scale and 
scheme dependence is often the most important source of theoretical error. An important 
question is: to what extension is the FT proven? In perturbation theory up to NNLO 
it has been explicitly checked to hold for many processes: if corrections exist we already 
know that they must be small (we stress that we are only considering totally inclusive 
processes). At all orders the most in depth discussions have been carried out in refs. [109], 
in particular for Drell-Yan processes. The LHC experiments offer a wonderful opportunity 
for testing the FT by comparing precise theoretical predictions with accurate data on a 
large variety of processes (for a recent review, see, for example, ref. [110]). 

A great effort has been and is being devoted to the theoretical preparation and inter- 
pretation of the LHC experiments. For this purpose very difficult calculations are needed 
at NLO and beyond because the strong coupling, even at the large Q 2 values involved, is 
not that small. New powerful techniques for amplitude calculations have been developed. 
An interesting development at the interface between string theory and QCD is twistor 
calculus. A precursor work was the Parke- Taylor result in 1986 [111] on the amplitudes 
for n incoming gluons with given helicities [112]. Inspired by dual models, they derived 
a compact formula for the maximum non vanishing helicity violating amplitude (with n- 
2 plus and 2 minus helicities) in terms of spinor products. Using the relation between 
strings and gauge theories in twistor space Witten developed in '03 [113] a formalism in 
terms of effective vertices and propagators that allows to compute all helicity amplitudes. 
The method, alternative to other modern techniques for the evaluation of Feynman dia- 
grams [115], leads to very compact results. Since then rapid progress followed (for reviews, 
see [114]): the method was extended to include massless external fermions [116] and also 
external EW vector bosons [117] and Higgs particles [118]. The level already attained is 
already important for multijet events at the LHC. The study of loop diagrams came next. 
The basic idea is that loops can be fully reconstructed from their unitarity cuts. First 
proposed by Bern, Dixon and Kosower [119] the technique was revived by Britto, Cachazo 
and Feng [120] and then perfected by Ossola, Papadopoulos and Pittau [121] and further 
extended to massive particles in ref. [122]. For a recent review of these new methods see 
ref. [123]. In parallel also the activity on event simulation has received a big boost from 
the LHC preparation (see, for example, the review [124]). Powerful techniques for the 
generation of numerical results at NLO for processes with complicated final states have 
been developed: the matching of matrix element calculation together with the modeling 
of parton showers has been realised in packages like Black Hat [125] (on-shell methods for 
loops) used in association with Sherpa [126] (for real emission), or POWHEG BOX [127], 
or aMCONLO [129], the automated version of the general framework MC@NLO [128]. 
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In a complete simulation the matrix element calculation, improved by resummation of 
large logs, provides the hard skeleton (with large branchings) while the parton shower 
is constructed by a sequence of factorized collinear emissions fixed by the QCD splitting 
functions. In addition, at low scales a model of hadronization completes the simulation. 
The importance of all the components, matrix element, parton shower and hadronization 
can be appreciated in simulations of hard events compared with Tevatron and LHC data. 

A partial list of examples of recent NLO calculations in pp collisions, obtained with 
these techniques is: W + 3 jets [130], Z, 7* + 3 jets [131], W, Z + 4 jets [132], W + 
5 jets [133], ttbb [134], it + 2 jets [135], it W [136], WW+ 2 jets [137], WWbb [138], 
bbbb [139] etc. Here in the following we present in more detail a number of important, 
simplest, examples without any pretension to completeness. 

2.9.1 Vector Boson Production 

Drell-Yan processes which include lepton pair production via virtual 7, W or Z ex- 
change, offer a particularly good opportunity to test QCD. This process, among those 
quadratic in parton densities with a totally inclusive final state, is perhaps the simplest 
one from a theoretical point of view. The large scale is specified and measured by the in- 
variant mass squared Q 2 of the lepton pair which itself is not strongly interacting (so there 
are no dangerous hadronization effects). The QCD improved parton model leads directly 
to a prediction for the total rate as a function of Q 2 . The value of the LO cross-section 
is inversely proportional to the number of colours Nq because a quark of given colour can 
only annihilate with an antiquark of the same colour to produce a colourless lepton pair. 
The order a s (Q 2 ) corrections to the total rate were computed long ago and found to be 
particularly large, when the quark densities are defined from the structure function F2 
measured in DIS at q 2 = —Q 2 [86,140]. The ratio <J corr jo^o of the corrected and the Born 
cross-sections, was called K-factor [141], because it is almost a constant in rapidity. In 
recent years also the NLO full calculation of the K-factor was completed, a very remark- 
able calculation [142]. Over the years the QCD predictions for W and Z production, a 
better testing ground than the older fixed target Drell-Yan experiments, have been tested 
at CERN SppS and Tevatron energies and now at the LHC. Q ~ mw,z is large enough 
to make the prediction reliable (with a not too large K-factor) and the ratio \fr = Qj \fs 
is not too small. Recall that in lowest order Xxx^s = Q 2 so that the parton densities are 
probed at x values around y/r. We have \fr = 0.13 — 0.15 (for W and Z production, 
respectively) at x/s = 630 GeV (CERN SppS Collider) and y/r = 0.04 - 0.05 at the Teva- 
tron. At the LHC at 8 TeV or at 14 TeV one has \Jt ~ 10 or ~ 6 -10 -3 , respectively. A 
comparison of the experimental total rates for W, Z with the QCD predictions at hadron 
colliders [143] is shown in Fig. 25. It is also important to mention that the cross-sections 
for di-boson production (i.e. WW, WZ, ZZ, Wj, Zj) have been measured at the Tevatron 
and the LHC and are in fair agreement with the SM prediction [144] . The typical precision 
is comparable to or better than the size of NLO corrections. 
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Figure 25: Data vs. theory for W and Z production at hadron colliders [143]. 



The calculation of the W/Z px distribution is a classic challenge in QCD. For large p?, 
for example p T ~ o(m w ), the px distribution can be reliably computed in perturbation 
theory, which was done up to NLO in the late '70's and early '80's [145]. A problem arises in 
the intermediate range Aqcd « Pt « %, where the bulk of the data is concentrated, 
because terms of order a s (p T ) log m^ v /p T become of order 1 and should included to all 
orders [146]. At order a s we have: 



lg. (1 + ^ (p | ) + £ log !| + (198) 



where A, B, C, D are coefficients of order a s . The "+" distribution is defined in complete 
analogy with eq. (187): 



2 2 

Ptmax ("Ptmax 



g(z)f(z) + dz = / [g(z)-g(0)]f(z)dz (199) 



o 



The content of this, at first sight mysterious, definition is that the singular "+" terms do 
not contribute to the total cross-section. In fact for the cross-section the weight function 
is g(z) = 1 and we obtain: 



Ptmax 



a = a [(l + A) + / D(z)dz] (200) 



The singular terms, of infrared origin, are present at the non completely inclusive level but 
disappear in the total cross-section. Solid arguments have been given that these singu- 
larities exponentiate. Explicit calculations in low order support the exponentiation which 
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leads to the following expression: 




1 da 



I 



<fb 
Arc 



exp (— ib ■ px) (1 + A) exp S(b) 



(201) 



with: 
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At large pt the perturbative expansion is recovered. At intermediate px the infrared pt 
singularities are resummed (the Sudakov log terms, which are typical of vector gluons, are 
related to the fact that for a charged particle in acceleration it is impossible not to radiate, 
so that the amplitude for no soft gluon emission is exponentially suppressed). A delicate 
procedure for matching perturbative and resummed terms is needed [147]. However this 
formula has problems at small pt, for example, because of the presence of a s under the 
integral for S(b): presumably the relevant scale is of order k T . So it must be completed 
by some non perturbative ansatz or an extrapolation into the soft region [146]. All the 
formalism has been extended to NLO accuracy [148], where one starts from the perturbative 
expansion at order a 2 sl and generalises the resummation to also include NLO terms of order 
a s (p T ) 2 log m^y/p^. The comparison with the data is very impressive. In Fig. 26 we see the 
Pt distribution as predicted in QCD (with a number of variants that mainly differ in the 
approach to the soft region) compared with some recent data at the Tevatron [149]. The 
W and Z pt distributions have also been measured at the LHC and are in fair agreement 
with the theoretical expectation [150]. 

The rapidity distributions of the produced W and Z have also been measured with fair 
accuracy at the Tevatron and at the LHC and predicted at NLO [151]. As a representative 
example of great significance we show in Fig. 27 the combined LHC results for the W 
charge asymmetry (defined as A ~ (W + — W~)/ (W + + W~) ) as a function of the pseudo 
rapidity rj [152]. These data combine the ATLAS and CMS results at smaller values of rj 
with those of the LHCb experiments at larger r\ (in the forward direction). This is a very 
important input for the disentangling of the different quark parton densities. 

2.9.2 Jets at Large Transverse Momentum 

Another simple and important process at hadron colliders is the inclusive production 
of jets at large energy a/s and transverse momentum p?- A comparison of the data with 
the QCD NLO predictions in pp or pp collisions is shown in Fig. 28 [153]. This is a 
particularly significant test because the rates at different c.o.m. energies and, for each 
energy, at different values of pt, span over many orders of magnitude. This steep behaviour 
is determined by the sharp falling of the parton densities with increasing x. Also, the 
corresponding values of y/s and px are large enough to be well inside the perturbative 
region. The overall agreement of the data from ISR, UA1,2, STAR (at RHIC), CDF/DO, 
and now ATLAS/CMS is indeed spectacular. In fact, the uncertainties on the resulting 
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Figure 26: QCD predictions for the W Pt distribution compared with recent DO data at the Tevatron (^/s 
= 1.8 TeV) [149] [148]. 



experiment /theory ratio, due to systematics and to ambiguities on parton densities, value 
of a s , scale choice and so on, which can reach a factor of 2-3, are much smaller than the 
spread of the cross-section values over many orders of magnitude. 

Similar results also hold for the production of photons at large px- The ATLAS data 
[154], shown in Fig. 29, are in fair agreement with the theoretical predictions. For the 
same process less clear a situation was found with fixed target data. Here, first of all, the 
experimental results show some internal discrepancies. Also, the pt accessible values being 
smaller, the theoretical uncertainties are larger. 

2.9.3 Heavy Quark Production 

We now discuss heavy quark production at colliders. The totally inclusive cross sections 
have been known at NLO since a long time [156]. The resummation of leading and next 
to leading logarithmically enhanced effects in the vicinity of the threshold region have also 
been studied [157]. The bottom production at the Tevatron has for some time represented 
a problem: the total rate and the pt distribution of b quarks observed at CDF and DO 
appeared in excess of the prediction, up to the largest measured values of px [158]. But 
this is a complicated problem, with different scales being present at the same time: \/s, 
Pt, frib. Finally the discrepancy has been solved by better taking into account a number 
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Figure 27: The combined LHC results for the W charge asymmetry (defined as A ~ (W + — W )/ (W + 
W~ ) ) as a function of the pseudo rapidity r\ [152]. 



of small effects from resummation of large logarithms, the difference between b hadrons 
and b partons, the inclusion of better fragmentation functions etc. [159]. At present the 
LHC data on b production are in satisfactory agreement with the theoretical predictions 
(Fig. 30 [160]). 

The top quark is really special: its mass is of the order of the Higgs VEV or its Yukawa 
coupling is of order 1 (in this sense it is the only "normal" case among all quarks and 
charged leptons). Due to its heavy mass it decays so fast that it has no time to be bound 
in a hadron: thus it can be studied as a quark. It is very important to determine its 
mass and couplings for different precision predictions of the SM. Perhaps the top quark 
could be particularly sensitive to new heavy states or have a connection to the Higgs 
sector in beyond the SM theories. Thus top physics has attracted much attention both 
from the experimental side, at hadron colliders, and from the theoretical point of view. 
In particular, the top-antitop inclusive cross-section has been measured in pp collisions at 
the Tevatron [161] and now in pp collisions at the LHC [162]. The QCD prediction is 
completely known at NLO [156] and almost completed at NNLO [163]. The agreement 
of theory and experiment is good for the best available parton density functions together 
with the values of a s and of m t measured separately (the top mass is measured from the 
invariant mass of the decay products), as can be seen from Fig. 31 [163]. The mass of 
the top (and the value of a s ) can be determined from the cross section, assuming that 
QCD is correct, and compared with the more precise value from the decay final state. The 
value of the pole top mass derived in ref. [164] from the cross-section data, using the best 
available parton densities with the correlated value of a s , is: mP t ole = 173.3 ± 2.8 GeV to 
be compared with the value measured at the Tevatron by the CDF and DO collaborations 
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Figure 28: Jet production cross-section at pp or pp colliders, as function of px [153]. The theory predictions 
are from NLO perturbative calculations with state-of-the-art parton densities with the corresponding value 
of a s plus a non perturbative correction factor due to hadronization and the underlying event, obtained 
using Monte Carlo event generators 



m exp _ 2 -l. O Q GeV. This quoted error is clearly too optimistic especially if one would 
identify this value with the pole mass which it resembles to. This error is only adequate 
within the specific procedure used by the experimental collaborations to define their mass 
(including which Montecarlo, which assumptions on higher order terms, non perturbative 
effects etc). The problem is how to export this value in other processes. Leaving aside the 
thorny issue of the precise relation of m e t xp with mP° le it is clear that there is a good overall 
consistency. 

The inclusive forward-backward asymmetry, Afb, m the ti rest frame has been mea- 
sured by both the CDF [165] and DO [166] collaborations and found to be in excess of the 
SM prediction, by about 2 a [167]. For CDF the discrepancy increases at large ti invariant 
mass and reaches about 2.5 a for M t i > 450 GeV. A related observable at the LHC is the 
charge asymmetry in tt production, Ac- In contrast to Afb, the combined value of Ac 
reported by ATLAS [168] and CMS [169] agrees with the SM, within the still limited pre- 
cision of the data. Is this a real discrepancy? It is too early to say, but it is an interesting 
effect that indeed is receiving much attention [170]. 
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Figure 29: Single photon production in pp colliders as function of px [154] 



2.9.4 Higgs Boson Production 

We now turn to the discussion of the SM Higgs inclusive production cross-section [171]. 
The most important Higgs production modes are gluon fusion, vector boson fusion, Higgs 
strahlung, and associated production with top quark pairs. Some typical Feynman dia- 
grams for those different modes are depicted in Fig. 32. The predicted rates are shown in 
Fig. 33 [172]. 

The most important channel at the LHC is Higgs production via g + g — )■ H. The 
amplitude is dominated by the top quark loop [173]. The NLO corrections turn out to 
be particularly large [174], as seen in Fig. 34. Higher order corrections can be computed 
either in the effective lagrangian approach, where the heavy top is integrated away and 
the loop is shrunk down to a point [175] (the coefficient of the effective vertex is known to 
a s accuracy [176]), or in the full theory. At the NLO the two approaches agree very well 
for the rate as a function of m# [177]. The NNLO corrections have been computed in the 
effective vertex approximation [178] (see Fig. 34). Beyond fixed order resummation of large 
logs were carried out [179]. Also the NLO EW contributions have been computed [180]. 
Rapidity (at NNLO) [181] and pt distributions (at NLO) [183] have also been evaluated. 
At smaller pt the large logarithms [Log(pr /it^h)}"' have been resummed in analogy with 
what was done long ago for W and Z production [184]. For additional recent works on 
Higgs physics at colliders see, for example, [185]. 
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Figure 30: The b production pt distribution at the LHC [160]. 




At different places in the previous pages we have seen examples of resummation of large 
logs. This is a very important chapter of modern QCD. The resummation of soft gluon logs 
enter in different problems and the related theory is subtle. I refer the reader here to some 
recent papers where additional references can be found [186]. A particularly interesting 
related development has to do with the so called non global logs (see, for example, [187]). 
If in the measurement of an observable some experimental cuts are introduced, which is a 
very frequent case, then a number of large logs can arise from the corresponding breaking 
of inclusiveness. It is also important to mention the development of software for the 
automated implementation of resummation (see, for example, [188]). 

2.10 Measurements of a s 

Very precise and reliable measurements of a s (mz) are obtained from e + e~ colliders (in 
particular LEP), from deep inelastic scattering and from the hadron Colliders (Tevatron 
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Figure 32: Representative Feynman diagrams for the Higgs production cross-section mechanisms: a) gluon 
fusion; b) Vector boson fusion (V=W, Z); c) Higgsstrahlung from a Z boson (an analogue diagram can be 
drawn for the W boson); d) tt associated production 



and LHC). The "official" compilation due to Bethke [189] and included in the 2012 edition 
of the PDG [14] is reproduced here in figs. 35. The agreement among so many different 
ways of measuring a s is a strong quantitative test of QCD. However for some entries the 
stated error is taken directly from the original works and is not transparent enough as 
seen from outside (e.g. the lattice determination). In my opinion one should select few 
theoretically cleanest processes for measuring a s and consider all other ways as tests of 
the theory. Note that in QED a is measured from one single very precise and theoretically 
clean observable (one possible calibration process is at present the electron g-2 [270]). The 
cleanest processes for measuring a s are the totally inclusive ones (no hadronic corrections) 
with light cone dominance, like Z decay, scaling violations in DIS and perhaps r decay 
(but, for r, the energy scale is dangerously low). We will review these cleanest methods 
for measuring a s in the following. 



2.10.1 a s from e + e colliders 

The totally inclusive processes for measuring a s at e + e _ colliders are hadronic Z decays 
(Ri, ah, <Ji, Yz) and hadronic r decays. As we have seen in Sect. 2.7.1, for a quantity like 
Ri we can write a general expression of the form: 

T(Z,t -» hadrons) EW 

R\ - -^t? — n + — r H I 1 + "qcd + o NP ) (206) 

1 (Z, r leptons) 

where R EW is the electroweak-corrected Born approximation, 8qcd, 8np are the pertur- 
bative (logarithmic) and non perturbative (power suppressed) QCD corrections. For a 
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Figure 33: The production cross sections at the LHC for a Higgs with mass Mjj ~ 125 GeV different cm. 
energies [172]. 
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Figure 34: The Higgs gluon fusion cross section in LO, NLO and NLLO [182]. 



measurement of a s (in the following we always refer to the MS definition of a s ) at the 
Z resonance peak one can use all the information from Ri, V z = 31^ + Th + Ti nv and 
ap = 12ttT iT p / (m 2 z T 2 z ) (F=h or 1). In the past the measurement from Ri was preferred 
(by itself it leads to a s (mz) = 0.1226 ± 0.0038, a bit on the large side) but after LEP 
there is no reason for this preference. In all these quantities a s enters through I\, but 
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Figure 35: Left: Summary of measurements of cx s {mz), used as input for the world average value of 
refs. [14,189]. The yellow band is the proposed average: a s (m z ) — 0.1184 ± 0.0007. Right: Summary of 
measurements of a s as a function of the respective energy scale Q 

the measurements of, say, Tz, Ri and o\ are really independent as they are affected by an 
entirely different systematics: Y z is extracted from the line shape, Ri and ai are measured 
at the peak but R\ does not depend on the absolute luminosity while <j\ does. The most 
sensitive single quantity is 07. It gives a s {m z ) = 0.1183 ± 0.0030. The combined value 
from the measurements at the Z (assuming the validity of the SM and the observed Higgs 
mass) is [190]: 

a s (m z ) = 0.1187 ± 0.0027 (204) 

By adding all other electroweak precision electroweak tests (in particular mw) one similarly 
finds [213]: 

a s (m z ) = 0.1186 ± 0.0026 (205) 

These results have been obtained from the 8qcd expansion up to and including the C3 term 
of order a 3 . But by now the c 4 term (NNNLO!) has also been computed [68] for inclusive 
hadronic Z and r decay. This remarkable calculation of about 20.000 diagrams, for the 
inclusive hadronic Z width, led to the result, for rif = 5 and a s = a s {mz)l^'- 

S QCD = [1 + a s + 0.76264 a 2 s - 15.49 a 3 s - 68.2 a 4 s + . . .] (206) 

This result can be used to improve the value of a s (mz) from the EW fit given in eq. (205) 
that becomes: 

a s (m z ) = 0.1190 ± 0.0026 (207) 
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Note that the error shown is dominated by the experimental errors. Ambiguities from 
higher perturbative orders [191], from power corrections and also from uncertainties on the 
Bhabha luminometer (which affect 07^) [192] are very small. In particular, having now 
fixed tuh does not decrease the error significantly [193, 194]. The main source of error is 
the assumption of no new physics, for example in the Zbb vertex that could affect the 
prediction. 

We now consider the measurement of a s (mz) from r decay. R T has a number of 
advantages that, at least in part, tend to compensate for the smallness of m T = 1.777 GeV. 
First, R T is maximally inclusive, more than R e + e -(s), because one also integrates over all 
values of the invariant hadronic squared mass: 

R T = 1 r (1 - ^fl m U T (s) (208) 
ii J m l T m l T 

As we have seen, the perturbative contribution is now known at NNNLO [68]. Analyticity 
can be used to transform the integral into one on the circle at \s\ = m 2 .: 

If ds 



R r = 7^f ^(1 - — ) 2 n T (s) (209) 

2m J\ s \ =m 2 mi m 2 . 

Also, the factor (1 — ^fj) 2 is important to kill the sensitivity the region Re[s] = m 2 . 
where the physical cut and the associated thresholds are located. Still the sensitivity to 
hadronic effects in the vicinity of the cut is a non negligible source of theoretical error that 
the formulation of duality violation models try to decrease. But the main feature that 
has attracted attention to r decays for the measurement of a s {mz) is that even a rough 
determination of Aqcd at a low scale Q ~ m T leads to a very precise prediction of a s at the 
scale mz, just because in log Q/Aqcd the value of Aqcd counts less and less as Q increases. 
The absolute error on a s shrinks by a factor of about one order of magnitude going from 
a s {m T ) to a s (mz)- Still I find a little suspicious that to obtain a better measurement of 
ots{mz) you have to go down to lower and lower energy scales. And in fact, in general, 
one finds that the decreased control of higher order perturbative and of non perturbative 
corrections makes the apparent advantage totally illusory. For a s from R T the quoted 
amazing precision is obtained by taking for granted that corrections suppressed by 1/m 2 
are negligible. The argument is that in the massless theory, the light cone expansion is 
given by: 

ZERO < Q 4 > < Q 6 > , om , 
+ c 4 — + c 6 — +••• (210) 



>NP 



ml mi ' m 6 



In fact there are no dim-2 Lorentz and gauge invariant operators. For example, g^g 11 
is not gauge invariant. In the massive theory, the ZERO is replaced by the light quark 
mass-squared m? . This is still negligible if m is taken as a lagrangian mass of a few MeV. 
If on the other hand the mass were taken to be the constituent mass of order Aqcd-, 
this term would not be negligible at all and would substantially affect the result (note 
that a s (m T )/n ~ 0.1 ~ (0.6 GeV/m T ) 2 and that Aq C d for 3 flavours is large). The 
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principle that coefficients in the operator expansion can be computed from the perturbative 
theory in terms of parton masses has never been really tested (due to ambiguities on the 
determination of condensates) and this particular case with a ZERO there is unique in 
making the issue crucial. Most distinguished people believe the optimistic version. I am 
not convinced that the gap is not filled up by ambiguities of Q(k 2 Q CD /m 2 ) from S pert [195]. 

There is a vast and sophisticated literature on a s from r decay. Unbelievably small 
errors are obtained in one or the other of several different procedures and assumptions that 
have been adopted to end up with a specified result. With time there has been an increasing 
awareness on the problem of controlling higher orders and non perturbative effects. In 
particular fixed order perturbation theory (FOPT) has been compared to resummation of 
leading beta function effects in the so called contour improved perturbation theory (CIPT). 
The results are sizably different in the two cases and there have been many arguments in the 
literature on which method is best. One important progress comes from the experimental 
measurement of moments of the r decay mass distributions, defined by modifying the 
weight function in the integral in eq.(208). In principle one can measure a s from the sum 
rules obtained from different weight functions that emphasize different mass intervals and 
different operator dimensions in the light cone operator expansion. A thorough study of 
the dependence of the measured value of a s on the choice of the weight function and in 
general of higher order and non perturbative corrections has appeared in ref. [196] and I 
advise the interested reader to look at that paper and the references therein. 

We consider here the recent evaluations of a s from r decay based on the NNNLO pertur- 
bative calculations [68] and different procedures for the estimate of all sorts of corrections. 
From the papers given in refs. [197] we obtain an average value and error that agrees with 
the Erler and Langacker values given in the PDG'12 [14]: 

a s (m T ) = 0.3285 ±0.018 (211) 

or 

a s (m z ) = 0.1194 ±0.0021 (212) 

In any case, one can discuss the error, but what is true and remarkable, is that the central 
value of a s from r decay, obtained at very small Q 2 , is in good agreement with all other 
precise determinations of a s at more typical LEP values of Q 2 . 

2.10.2 a s from Deep Inelastic Scattering 

In principle DIS is expected to be an ideal laboratory for the determination of a s but 
in practice the outcome is still to some extent open. QCD predicts the Q 2 dependence 
of F(x, Q 2 ) at each fixed x, not the x shape. But the Q 2 dependence is related to the x 
shape by the QCD evolution equations. For each x-bin the data allow to extract the slope 
of an approximately straight line in dlogF(x, Q 2 ) / dlogQ 2 : the log slope. The Q 2 span and 
the precision of the data are not much sensitive to the curvature, for most x values. A 
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single value of Aqqd must be fitted to reproduce the collection of the log slopes. For the 
determination of a s the scaling violations of non-singlet structure functions would be ideal, 
because of the minimal impact of the choice of input parton densities. We can write the 
non-singlet evolution equations in the form: 

d j r,, A a s{t) f 1 dyF(y,t) ,x 

where P qq is the splitting function. At present NLO and NNLO corrections are known. It is 
clear from this form that, for example, the normalization error on the input density drops 
away, and the dependence on the input is reduced to a minimum (indeed, only a single 
density appears here, while in general there are quark and gluon densities). Unfortunately 
the data on non-singlet structure functions are not very accurate. If we take the difference 
of data on protons and neutrons, F p — F n , experimental errors add up in the difference and 
finally are large. The F 3u ^ data are directly non-singlet but are not very precise. Another 
possibility is to neglect sea and glue in F 2 at sufficiently large x. But by only taking data at 
x > x one decreases the sample, introduces a dependence on x and an error from residual 
singlet terms. A recent fit to non singlet structure functions in electro- or muon-pro duct ion 
extracted from proton and deuterium data, neglecting sea and gluons at x > 0.3 (error to 
be evaluated) has led to the results [198]: 

a s (m z ) = 0.1148 ±0.0019(exp) + ? (NLO) (214) 
a s (m z ) = 0.1134 ±0.0020(ea;p) + ? (NNLO) (215) 

The central values are rather low and there is not much difference between NLO and 
NNLO. The question marks refer to the uncertainties from the residual singlet component 
at x > 0.3 and also to the fact that the old BCDMS data, whose systematics has been 
questioned, are very important at x > 0.3 and push the fit towards small values of a s . 

When one measures a s from scaling violations in F 2 measured with e or /x beams, the 
data are abundant, the statistical errors are small, the ambiguities from the treatment of 
heavy quarks and the effects of the longitudinal structure function F L can be controlled, 
but there is an increased dependence on input parton densities and especially a strong 
correlation between the result on a s and the adopted parametrization of the gluon density. 
In the following we restrict our attention to recent determinations of a s from scaling vio- 
lations at NNLO accuracy, as, for example, those in refs. [199,200] that report the results, 
in the order: 

a s (m z ) = 0.1134 ±0.0011(exp) + ? (216) 
a s (m z ) = 0.1158 ±0.0035 (217) 

In the first line my question mark refers to the issue of the o^-gluon correlation. In fact 
a s tends to slide towards low values (a s ~ 0.113 — 0.116) if the gluon input problem is not 
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fixed. Indeed, in the second line from ref. [200], the large error also includes an estimate 
of the ambiguity from the gluon density parametrization. One way to restrict the gluon 
density is to use the Tevatron high p? jet data to fix the gluon parton density at large x 
that, via the momentum conservation sum rule, also constrains the small x values of the 
same density. Of course in this way one has to go outside the pure domain of DIS. Also, 
the jet rates have been computed at NLO only. In a simultaneous fit of a s and the parton 
densities from a set of data that, although dominated by DIS data, also contains Tevatron 
jets and Drell- Yan production, the result was [201]: 

a s (m z ) =0.1171 ±0.0014 + ? (218) 

The authors of ref. [201] attribute their larger value of a s to a more flexible parametrization 
of the gluon and the inclusion of Tevatron jet data that are important to fix the gluon at 
large x. An alternative way to cope with the gluon problem is to drastically suppress 
the gluon parametrization rigidity by adopting the neural network approach. With this 
method, in ref. [202], from DIS data only, treated at NNLO accuracy, the following value 
was obtained: 

a s (m z ) = 0.1166 ± 0.0008(exp) ± 0.0009(U) + ? (219) 

where the stated theoretical error is that quoted by the authors within their framework, 
while the question mark has to do with possible additional systematics from the method 
adopted. Interestingly, in the same approach, by also including the Tevatron jets and the 
Drell- Yan data not much difference is found: 

a s (m z ) = 0.1173 ± 0.0007(ea;p) ± 0.0009(t/i) + ? (220) 

We see that when the gluon input problem is suitably addressed the fitted value of a s is 
increased. 

As we have seen there is some spread of results, even among the most recent determi- 
nations based on NNLO splitting functions. We tend to favour determinations from the 
whole DIS set of data (i.e. beyond the pure non singlet case) and with attention paid 
to the gluon ambiguity problem (even if some non DIS data from Tevatron jets at NLO 
have to be included). A conservative proposal for the resulting value of a s from DIS, that 
emerges from the above discussion is something like: 

a s (m z ) = 0.1165 ±0.0020 (221) 

The central value is below those obtained from Z and r decays but perfectly compatible 
with those results. 

2.10.3 Recommended value of a s (m z ) 

According to my proposal to calibrate a s (m z ) from the theoretically cleanest and most 
transparent methods, identified as the totally inclusive, light cone operator expansion dom- 
inated processes, I collect here my understanding of the results: from Z decays and EW 
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precision tests, eq.(205): 



a s {m z ) = 0.1190 ±0.0026; 



(222) 



from scaling violations in DIS, eq.(221): 



a s {m z ) = 0.1165 ±0.0020; 



(223) 



from R T , eq.(212): 



a s {m z ) = 0.1194 ±0.0021. 



(224) 



If one wants to be on the safest side one can take the average of Z decay and DIS: 



This is my recommended value. If one adds to the average the rather conservative R T value 
and error given above in eq.212, that takes into account the dangerous low energy scale of 
the process, one obtains: 



Note that this is essentially coincident with the "official" average with a moderate increase 
of the error. 

2.10.4 Other a s (m z ) Measurements as QCD Tests 

There are a number of other determinations of a s that are important because they arise 
from qualitatively different observables and methods. Here I will give a few examples of 
the most interesting measurements. 

A classic set of measurements is from a number of infrared safe observables related 
to event rates and jet shapes in e + e~ annihilation. One important feature of these mea- 
surements is that they can be repeated at different energies in the same detector, like the 
JADE detector in the energy range of PETRA (most of the intermediate energy points 
in Fig. 35-right panel are from this class of measurements) or the LEP detectors from 
LEP1 to LEP2 energies. As a result one obtains a striking direct confirmation of the 
running of the coupling according to the renormalization group prediction. The pertur- 
bative part is known at NNLO [203] and resummations of leading logs arising from the 
vicinity of cuts and/or boundaries have been performed in many cases using effective field 
theory methods. The main problem of these measurements is the possible large impact 
of non perturbative hadronization effects on the result and therefore on the theoretical 
error. According to ref. [189] a summary result that takes into account the central val- 
ues and the spread from the JADE measurements, in the range 14 to 46 GeV, at PE- 
TRA is given by: a s (m z ) = 0.1172 ± 0.0051, while from the ALEPH data at LEP, in 
the range 90 to 206 GeV, the reported value [204] is a s {m z ) = 0.1224 ± 0.0039. It is 
amazing to note that among the related works there are a couple of papers by Abbate 
et al [205, 206] where an extremely sophisticated formalism is developed for the thrust 



a s {m z ) = 0.1174 ±0.0016. 



(225) 



a s {m z ) = 0.1184 ±0.0011. 



(226) 
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distribution, based on NNLO perturbation theory with resummations at NNNLL plus 
a data/theory based estimate of non perturbative corrections. The final quoted results 
are unbelievably precise: a s (mz) = 0.1135 ± 0.0011 from the tail of the Thrust distri- 
bution [205] and a s (mz) = 0.1140 ± 0.0015 from the first moment of the Thrust distri- 
bution [206]. I think that this is a good example of an underestimated error which is 
obtained within a given machinery without considering the limits of the method itself. 
Another allegedly very precise determination of a s (mz) is obtained from lattice QCD by 
several groups [207] with different methods and compatible results. A value that sum- 
marizes these different results is [14] a s {mz) = 0.1185 ± 0.0007. With all due respect to 
lattice people I think this small error is totally unplausible. But we have shown that a 
sufficiently precise measure of ot s {m z ) can be obtained, eqs. (225,226), by only using the 
simplest processes where the control of theoretical errors is maximal. One is left free to 
judge whether a further restriction of theoretical errors is really on solid ground. 
The value of A (for nj = 5) which corresponds to eq. (225) is: 

A 5 = 202 ± 18 MeV (227) 

while the value from eq. (226) is: 

A 5 = 213 ± 13 MeV (228) 

A is the scale of mass that finally appears in massless QCD. It is the scale where a s (A) 
is of order 1. Hadron masses are determined by A. Actually the p mass or the nucleon 
mass receive little contribution from the quark masses (the case of pseudoscalar mesons 
is special, as they are the pseudo Goldstone bosons of broken chiral invariance). Hadron 
masses would be almost the same in massless QCD. 

2.11 Conclusion 

We have seen that perturbative QCD based on asymptotic freedom offers a rich variety 
of tests and we have described some examples in detail. QCD tests are not as precise as 
for the electroweak sector. But the number and diversity of such tests has established a 
very firm experimental foundation for QCD as a theory of strong interactions. The physics 
content of QCD is very large and our knowledge, especially in the non perturbative domain, 

is still very limited but progress both from experiment (Tevatron, RHIC, LHC ) and 

from theory is continuing at a healthy rate. And all the QCD predictions that we were 
able to formulate and to test appear to be in very good agreement with experiment. 

The field of QCD appears as one of great maturity but also of robust vitality with many 
rich branches and plenty of new blossoms. I may mention the very exciting explorations of 
Supersymmetric extensions of QCD and the connections with string theory (for a recent 
review and a list of references see ref. [208]). In particular N=4 SUSY QCD (that is, with 
4 spinor charge generators) has a vanishing beta function and is loop finite. In the limit 
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N c — > oo with A = e 2 s Nc fixed planar diagrams are dominant. There is progress towards a 
solution of planar N=4 SUSY QCD. The large A limit corresponds by the AdS/CFT duality 
(Anti de Sitter/ Conformal Field Theory), a string theory concept, to the weakly coupled 
string (gravity) theory on AdS*, ' S§ (the 10 dimensions are compactified in a 5-dimensional 
Anti de Sitter space times a 5-dimensional sphere. By moving along this very tentative 
route one can transfer some results (assumed to be of sufficiently universal nature) from 
the computable weak limit of the associated string theory to the non perturbative ordinary 
QCD domain. Further away on this line there are studies on the N = 8 Supergravity, 
related to N = 4 SUSY Yang-Mills, which has been proven finite up to 4 loops. It could 
possibly lead to a finite field theory of gravity in 4 dimensions. 
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3 The Theory of Electroweak Interactions 



3.1 Introduction 

In this Chapter, we summarize the structure of the standard EW theory 4 and specify 
the couplings of the intermediate vector bosons W ± , Z and of the Higgs particle with the 
fermions and among themselves, as dictated by the gauge symmetry plus the observed 
matter content and the requirement of renormalizability. We discuss the realization of 
spontaneous symmetry breaking and of the Higgs mechanism. We then review the phe- 
nomenological implications of the EW theory for collider physics (that is we leave aside the 
classic low energy processes that are well described by the "old" weak interaction theory 
(see, for example, [212])). For this discussion we split the lagrangian into two parts by 
separating the terms with the Higgs field: 

C -£gauge ~\~ ^Higgs • (229) 

Both terms are written down as prescribed by the SU(2) {7(1) gauge symmetry and 
renormalizability, but the Higgs vacuum expectation value (VEV) induces the spontaneous 
symmetry breaking responsible for the non vanishing vector boson and fermion masses. 



3.2 The Gauge Sector 

We start by specifying £ ga uge, which involves only gauge bosons and fermions, according 
to the general formalism of gauge theories discussed in Chapter 1: 

1 3 1 - 

A=l 

This is the Yang-Mills lagrangian for the gauge group 577(2) ® {7(1) with fermion matter 
fields. Here 

B^ v = d^B u -d u B^ and = d^W? - dvWf - geABC W*W? (231) 

are the gauge antisymmetric tensors constructed out of the gauge field associated with 
i7(l), and corresponding to the three SU{2) generators; €abc are the group structure 
constants [see Eqs. (233,234)] which, for 577(2), coincide with the totally antisymmetric 
Levi-Civita tensor (with £123 = 1; recall the familiar angular momentum commutators). 
The normalization of the SU(2) gauge coupling g is therefore specified by Eqs. (231). As 
discussed in Sect. 1.5 the standard EW theory is a chiral theory, in the sense that ipL 
and behave differently under the gauge group (so that parity and charge conjugation 
non conservation are made possible in principle). Thus, mass terms for fermions (of the 
form ipLipR + h.c.) are forbidden in the symmetric limit. In the following by ipL,n we 

4 Some recent textbooks are listed in ref. [209]; see also refs. [210], [211] 
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mean column vectors, including all fermion types in the theory that span generic reducible 
representations of SU(2) ® U(l). In the absence of mass terms, there are only vector and 
axial vector interactions in the lagrangian and those have the property of not mixing i/jl 
and ipR. Fermion masses will be introduced, together with W and Z masses, by the 
mechanism of symmetry breaking. The covariant derivatives D^l,r are explicitly given 
by 



A=l 



i>L,R , (232) 



where t^ R and \j2Y LR are the SU(2) and U(l) generators, respectively, in the reducible 
representations ipL,R- The commutation relations of the SU(2) generators are given by 

[tttf] = * tAB C t C L and [t%t B R ] = ie ABC t c R . (233) 

We use the normalization as in eq. (11) [in the fundamental representation of SU(2)]. The 
electric charge generator Q (in units of e, the positron charge) is given by 

Q = t 3 L + 1/2 Y L = t% + 1/2 Y R . (234) 

Note that the normalization of the U(l) gauge coupling g' in (232) is now specified as a 
consequence of (234). Note that t R ipR = 0, given that, for all known quark and leptons, 
ipR is a singlet. But in the following, we keep t R ipR for generality, in case one day a non 
singlet right-handed fermion is discovered. 

3.3 Couplings of Gauge Bosons to Fermions 

All fermion couplings of the gauge bosons can be derived directly from Eqs. (230) and 
(232). The charged fields are described by W^' 2 , while the photon and weak neutral 
gauge boson are obtained from combinations of and B^. The charged- current (CC) 
couplings are the simplest. One starts from the W^ 2 terms in Eqs. (230) and (232) which 
can be written as: 

g^Wl+fWl) = g{[^ + it 2 )/V2^(Wi-iW 2 )/V2]+h.c.} 

= g{[(t + W-)/V2]+h.c.} , (235) 

where = t 1 ± it 2 and W ± = (W 1 ± iW 2 )/V2. By applying this generic relation to L 
and R fermions separately, we obtain the vertex 

V^ w = #7/, [{ttl V2){1 - 7B )/2 + (tyV2)(l + 75 )/2] W~ + h.c. (236) 

Given that tR = for all fermions in the SM, the charged current is pure V — A. In 
the neutral-current (NC) sector, the photon A^ and the mediator Z^ of the weak NC are 
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orthogonal and normalized linear combinations of and W 3 : 

= cos OwB^ + sin 6 W W* , 
Zf, = -sin9 w B^ + cos 9 W W 3 . (237) 

and conversely: 

W 3 = sin 9 w A il + cos 9^Z ll , 

Bp = cos 9\vA^ — sin 9-^rZ^ . (238) 

Equations (237) define the weak mixing angle 9w- We can rewrite the W 3 and B^ terms 
in Eqs. (230) and (232) as follows: 

gt 3 W 3 + sfY/ZBn = [gt 3 sin9 w + g'(Q-t 3 )cos9 w )A fM + 

+ [gt 3 cos 9 W ~g'(Q- t 3 ) sin 9 w ]Z tl , (239) 

where Eq. (234) for the charge matrix Q was also used. The photon is characterized by 
equal couplings to left and right fermions with a strength equal to the electric charge. Thus 
we immediately obtain 

g sin 6^ = g' cos 9w = e , (240) 

so that: 

tg 9 W = g'/g (241) 

Once 9w has been fixed by the photon couplings, it is a simple matter of algebra to derive 
the Z couplings, with the result 

V^z = tt^-H-^lO- - 7s) + 4(1 + 75) - 2Q sin 2 OwtyZ* , (242) 

Z COS u\y 

where V^z is a notation for the vertex. Once again, recall that in the minimal SM, t 3 R = 
and t 3 L = ±1/2. 

In order to derive the effective four-fermion interactions that are equivalent, at low 
energies, to the CC and NC couplings given in Eqs. (236) and (242), we anticipate that 
large masses, as experimentally observed, are provided for W ± and Z by £ffi g gs- For left- 
left CC couplings, when the momentum transfer squared can be neglected (with respect 
to m^y) in the propagator of Born diagrams with single W exchange (see, for example, the 
diagram for /1 decay in Fig. 3.1), from Eq. (236) we can write 

4 C ff ° - s4-07m(1 - 7 5 )t^W(l - T6)tl1>] • (243) 

By specializing further in the case of doublet fields such as v e — e" or v^ — pr, we obtain 
the tree-level relation of g with the Fermi coupling constant Gp precisely measured from 
\x decay [see Chapter 1, Eqs. (2), (3)]: 

G F /V2 = g 2 /8m 2 w . (244) 
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\x~ W 



Figure 36: The Born diagram for fi decay. 



f f 




a b 

Figure 37: Diagrams for (a) the W and (b) the Z widths in Born approximation. 



By recalling that g sinOw = e, we can also cast this relation in the form 

m w = /iBorn/ sin 9 W , (245) 

with 

/iBom = (na/y/2G F ) 1/2 ~ 37.2802 GeV , (246) 

where a is the fine-structure constant of QED (a = e 2 /47r = 1/137.036). 

In the same way, for neutral currents we obtain in Born approximation from Eq. (242) 
the effective four-fermion interaction given by 

C c ~ y/2 G FPo ^[-]^7 M [-# , (247) 

where 

[...]= t£(l - 75) + 4(1 + 7s) - 2Q sin 2 9 W (248) 

and 

p0 = m *1t*o ■ (249) 

771 z COS t>w 

All couplings given in this section are valid at tree level and are modified in higher orders 
of perturbation theory. In particular, the relations between mw and sin^jy [Eqs. (245) 
and (246)] and the observed values of p (p = p at tree level) in different NC processes, are 
altered by computable EW radiative corrections, as discussed in Sect. (3.11). 

The partial width T(W — > ff) is given in Born approximation by the simplest diagram 
in Fig. 3.2 and one readily obtains from Eq.(236) with tu = 0, in the limit of neglecting 
the fermion masses and summing over all possible /' for a given /: 

T(W -> ff) = N c °f^ = N C -^-, (250) 
67TV2 12snr6V 
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where Nq = 3 or 1 is the number of colours for quarks or leptons, respectively, and 
the relations Eqs.(240, 244) have been used. Here and in the following expressions for 
the Z widths the one loop QCD corrections for the quark channels can be absorbed in a 
redefinition of Nq- Nq —> 3[l + a s (mz)/ir+ ■■■]■ Note that the widths are particularly large 
because the rate already occurs at order g 2 or Gp. The experimental values of the W total 
width and the leptonic branching ratio (the average of e, /i and r modes) are [14], [213] 
(see Sect. (3.11)): 



2.085 ± 0.042 GeV, B(W ->• lv x ) = 10.80 ± 0.09. (251) 



The branching ratio B is in very good agreement with the simple approximate formula, 
derived from Eq.(250): 

B(W lvx) ~ — — ) 2W , ~ 10.8%. (252) 

v ' 2-3(1 + a s (m|)/vr) + 3 v ; 

The denominator corresponds to the sum of the final states d'u, s'c, e~v e , ii~v^ t~v t (for 
the definition of d' and s' see Eq. (291)). 

For £r = the Z coupling to fermions in Eq.(242) can be cast into the form: 



9 J, ni j j _ J 

v ipfiPfZ 2 c 

with: 



v */+f* = v^br^f 1 ^ ~ Ai^i^ , (253) 

Cub u ly 



9a = , 9v/9a = 1 - MQA sm 2 9 W . (254) 
and t 3 / = ±1/2 for up-type or down-type fermions. In terms of g^,v given in Eqs. (254) (the 
widths are proportional to (g 2 , +g\)), the partial width T(Z — > ff) in Born approximation 
(see the diagram in Fig. 3.2), for negligible fermion masses, is given by: 

T(Z^ff) = N c am 2 Z [l + (l-4|g / |sin 2 M 2 ] 
12 sm 2v w 

= iV cPo ^|[l + (l-4|Q / |sin 2 M 2 ]- (255) 

where po = / m 2 z cos 2 9w is given in Eq. (280). The experimental values of the Z total 
width and of the partial rates into charged leptons (average of e, /i and r), into hadrons 
and into invisible channels are [14], [213]: 

T z = 2.4952 ± 0.0023 GeV, 

Ti+i- = 83.984 ± 0.086 MeV, 

T h = 1744.4 ± 2.0 MeV, 

T inv = 499.0 ± 1.5 MeV. (256) 

The measured value of the Z invisible width, taking radiative corrections into account, 
leads to the determination of the number of light active neutrinos [14], [213]: 

N v = 2.9840 ± 0.0082, (257) 
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well compatible with the 3 known neutrinos u e , and v T \ hence there exist only the three 
known sequential generations of fermions (with light neutrinos), a result with important 
consequences also in astrophysics and cosmology. 

At the Z peak, besides total cross sections, various types of asymmetries have been mea- 
sured. The results of all asymmetry measurements are quoted in terms of the asymmetry 
parameter Af, defined in terms of the effective coupling constants, gy and g A , as: 

A = 2 |£g£ = 2 9 * V, / A f , A f FB = -A e A f . (258) 
' 9v 2 + 9 f A 1 + (^M) 2 4 ' 

The measurements are: the forward-backward asymmetry {Ap B = (3/4)A e Af), the tau 
polarization (A T ) and its forward backward asymmetry (A e ) measured at LEP, as well as 
the left-right and left-right forward-backward asymmetry measured at SLC (A e and Af, 
respectively). Hence the set of partial width and asymmetry results allows the extraction 
of the effective coupling constants: widths measure (g v + g\) and asymmetries measure 
9v/9A- 

The top quark is heavy enough that it can decay into a real bW pair, which is by far 
its dominant decay channel. The next mode, t — > sW , is suppressed in rate by a factor 
\V ts \ 2 ~ 1.710 -3 , see Eqs. (296-298). The associated width, neglecting nib effects but 
including 1-loop QCD corrections in the limit mw = 0, is given by (we have omitted a 
factor \V t b\ 2 that we set equal to 1) [214]: 

Gpm? . m?,,,,, m?, Ar 2a s (mz) .27T 2 5 N n . 

r(( ^ bW) = 8^1 (1 ~ 1% » (1 + 2 ^ M1 - ^tr 1 ^ - 5> + -I- (259) 

The top quark lifetime is so short, about 0.510~ 24 s, that it decays before hadronizing or 
forming toponium bound states. 



3.4 Gauge Boson Self-interactions 

The gauge boson self- interactions can be derived from the F^ u term in £ ga uge, by using 
Eq. (237) and W ± = (W 1 ±iW 2 ) / y/2. Defining the three-gauge-boson vertex as in Fig. 3.3 
(with all incoming lines), we obtain (V = 7, Z) 

Vw-w+v = igw-w+v[g»v(p - q)x + 9^\{r-p) v + g vX (q - r)J , (260) 

with 

9w-w+i = 9 sin6>vy = e and gw~w+z = 9 cos 0w . (261) 

Note that the photon coupling to the W is fixed by the electric charge, as imposed by 
QED gauge invariance. The ZWW coupling is larger by a cot 6w factor. This form of the 
triple gauge vertex is very special: in general, there could be departures from the above 
SM expression, even restricting us to Lorentz invariant, em gauge symmetric and C and P 
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Figure 38: The 3- and 4-gauge boson vertices. The cubic coupling is of order g, while the quartic one is of 
order g 2 . 

conserving couplings. In fact some small corrections are already induced by the radiative 
corrections. But, in principle, more important could be the modifications induced by some 
new physics effect. The experimental testing of the triple gauge vertices has been done 
mainly at LEP2 and at the Tevatron [215]. 




Figure 39: The lowest order diagrams for e + e — > W + W . 

As a particularly important example, the cross-section and angular distributions for 
the process e + e~ — > W + W~ have been studied at LEP2. In Born approximation the 
Feynman diagrams for the LEP2 process are shown in Fig. 3.4 [216]. Besides neutrino 
exchange which only involves the well established charged current vertex, the triple weak 
gauge vertices Vw-w+v appear in the 7 and Z exchange diagrams. The Higgs exchange 
is negligible because the electron mass is very small. The analytic cross section formula 
in Born approximation can be found, for example, in ref. [14] (in the section "Cross- 
section formulae for specific processes"). The experimental data are compared with the 
SM prediction in Fig. 40. The agreement, within the present accuracy, is good. Note 
that the sum of all three exchange amplitudes has a better high energy behaviour than 
its individual components. This is due to cancellations among the amplitudes implied by 
gauge invariance, connected to the fact that the theory is renormalizable (the cross-section 
can be seen as a contribution to the imaginary part of the e + e~ — > e + e~ amplitude). 

The quartic gauge coupling is proportional to g 2 eABcW B W c eADEW D W E . Thus in the 
term with A=3 we have 4 charged W's. For A=l or 2 we have 2 charged W's and 2 W 3 's, 
each W3 being a combination of 7 and Z according to Eq. (238). With a little algebra the 
quartic vertex can be cast in the form: 

Vwwvv — i9wwvv[^9^v9\p ~ 9^x9u P — 9^p9ux} , (262) 
where, fi and v refer to W + W + in the 4W vertex and to VV in the WWW case and: 
9wwww — 9' \ 9ww-y-y — — ^ 2 , 9ww~{Z = —eg cos9 w , gwwzz = —g 2 cos 2 9 w . (263) 

In order to obtain these result for the vertex the reader must duly take into account the 
factor of -1/4 in front of F 2 in the lagrangian and the statistical factors which are equal 
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Figure 40: The measured production cross section for e + e — > W + W compared to the SM and to 
fictitious theories not including trilinear gauge couplings, as indicated 

to 2 for each pair of identical particles (like W + W + or 77, for example). The quartic 
coupling, being quadratic in g, hence small, could not be directly tested so far. 



We now turn to the Higgs sector of the EW lagrangian [10]. Until recently this simplest 
realization of the EW symmetry breaking was a pure conjecture. But on July '12 the AT- 
LAS and CMS Collaborations at the CERN LHC have announced [217,218] the discovery 
of a particle with mass m# ~ 12QGeV that very much looks like the long sought Higgs 
particle. A more precise measurement of its couplings and the proof that its spin is zero 
are necessary before this identification can be confirmed. But the following description of 
the Higgs sector of the SM can now be read with this striking development in mind. 

The Higgs lagrangian is specified by the gauge principle and the requirement of renor- 
malizability to be 



where is a column vector including all Higgs fields; in general it transforms as a reducible 
representation of the gauge group SU(2)l <g) U{1). In the Minimal SM it is just a complex 
doublet. The quantities V (which include all coupling constants) are matrices that make the 



3.5 The Higgs Sector 



(264) 
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Yukawa couplings invariant under the Lorentz and gauge groups. The potential V(0V)j 
symmetric under SU{2) L <g> £7(1), contains, at most, quartic terms in so that the theory- 
is renormalizable: 

Vitpfi = + 2 A (^V) 2 (265) 

As discussed in Chapter 1, spontaneous symmetry breaking is induced if the minimum 
of V, which is the classical analogue of the quantum mechanical vacuum state, is not a 
single point but a whole orbit obtained for non-vanishing values. Precisely, we denote 
the vacuum expectation value (VEV) of 0, i.e. the position of the minimum, by v (which 
is a doublet): 



(O|0(x)|O) = 




(266) 



The reader should be careful that, for economy of notation, the same symbol is used for 
the doublet and for the only non zero component of the same doublet. The fermion mass 
matrix is obtained from the Yukawa couplings by replacing 4>(x) by v: 

M = ^ L M^ R + ^PrMHl , (267) 

with 

M = T • v . (268) 

In the MSM, where all left fermions ipL are doublets and all right fermions ipR are singlets, 
only Higgs doublets can contribute to fermion masses. There are enough free couplings in 
T so that one single complex Higgs doublet is indeed sufficient to generate the most general 
fermion mass matrix. It is important to observe that by a suitable change of basis we can 
always make the matrix Ai Hermitian (so that the mass matrix is 75-free) and diagonal. 
In fact, we can make separate unitary transformations on ipi and ipR according to 

1/ L = Uif> L , i>' R = WiP R (269) 

and consequently 

M -t M' = U ] MW . (270) 

This transformation produces different effects on mass terms and on the structure of the 
fermion couplings in £ symm , because both the kinetic terms and the couplings to gauge 
bosons do not mix L and R spinors. The combined effect of these unitary rotations leads to 
the phenomenon of mixing and, generically, to flavour changing neutral currents (FCNC), 
as we shall see in Sect. 3.6. If only one Higgs doublet is present, the change of basis that 
makes M. diagonal will at the same time diagonalize the fermion-Higgs Yukawa couplings. 
Thus, in this case, no flavour-changing neutral Higgs vertices are present. This is not 
true, in general, when there are several Higgs doublets. But one Higgs doublet for each 
electric charge sector i.e. one doublet coupled only to u-type quarks, one doublet to d- 
type quarks, one doublet to charged leptons and possibly one for neutrino Dirac masses, 
would also be all right, because the mass matrices of fermions with different charges are 
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diagonalized separately. For several Higgs doublets in a given charge sector it is also possible 
to generate CP violation by complex phases in the Higgs couplings. In the presence of six 
quark flavours, this CP-violation mechanism is not necessary. In fact, at the moment, the 
simplest model with only one Higgs doublet could be adequate for describing all observed 
phenomena. 

We now consider the gauge-boson masses and their couplings to the Higgs. These effects 
are induced by the (D^Y (D^cfi) term in £m g gs [Eq. (264)], where 



d, + igJ2 tA K + i 9'( Y /2)B, 



A=l 



(271) 



Here t A and Y/2 are the SU(2) <S> U(l) generators in the reducible representation spanned 
by 0. Not only doublets but all non-singlet Higgs representations can contribute to gauge- 
boson masses. The condition that the photon remains massless is equivalent to the condi- 
tion that the vacuum is electrically neutral: 

Q\v) = (? + ±Y)\v) = 0. (272) 
We now explicitly consider the case of a single Higgs doublet: 

The charged W mass is given by the quadratic terms in the W field arising from £ffi gg s, 
when <f)(x) is replaced by v in Eq.(266). By recalling Eq.(235), we obtain 

m 2 w W+W~^ = g 2 \(t + v/V2)\ 2 W+W-v , (274) 

whilst for the Z mass we get [recalling Eqs. (237-239)] 

-m|Z M Z" = \[gcos9 w t 3 - g' sm9 w (Y/2)]v\ 2 Z^ , (275) 

where the factor of 1/2 on the left-hand side is the correct normalization for the definition 
of the mass of a neutral field. By using Eq. (272), relating the action of t 3 and Y/2 on the 
vacuum v, and Eqs. (241), we obtain 

i m | = (gcosOw + g' sin6 w ) 2 \t 3 v\ 2 = (g 2 / cos 2 d w )\t 3 v\ 2 . (276) 
For a Higgs doublet, as in eq.(273), we have 

\t+ v \ 2 = v 2 , \t 3 v\ 2 = l/Av 2 , (277) 



so that 



1 oV 

2 L 2 2 2 y u /c,no\ 

m w = -g v , m z = . (278) 

2 2 cos z ow 
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Note that by using Eq. (244) we obtain 

v = 2- 3 / 4 G~ 1/2 = 174.1 GeV . (279) 
It is also evident that for Higgs doublets 

Po = 2 m \ ft = 1 • (280) 
m z z cos^ u w 

This relation is typical of one or more Higgs doublets and would be spoiled by the existence 
of Higgs triplets etc. In general, 

"° - mm (281) 

for several Higgs bosons with VEVs Vi, weak isospin ti, and z-component tf . These results 
are valid at the tree level and are modified by calculable EW radiative corrections, as 
discussed in Sect. 3.11. 

The measured values of the W (combined from the LEP and Tevatron experiments) 
and Z masses (from LEP) are [14], [213]: 

mw = 80.385 ± 0.015 GeV, m z = 91.1876 ± 0.0021 GeV. (282) 

In the minimal version of the SM only one Higgs doublet is present. Then the fermion- 
Higgs couplings are in proportion to the fermion masses. In fact, from the fermion / 
Yukawa couplings gsffifL&fii + h.c), the mass mj is obtained by replacing <fi by v, so 
that 77i f — g<hffV. In the minimal SM three out of the four Hermitian fields are removed 
from the physical spectrum by the Higgs mechanism and become the longitudinal modes of 
W + , W~ , and Z. The fourth neutral Higgs is physical and should presumably be identified 
with the newly discovered particle at ~ 126 GeV. If more doublets are present, two more 
charged and two more neutral Higgs scalars should be around for each additional doublet. 

The couplings of the physical Higgs H can be simply obtained from £Hi ggs , by the 
replacement (the remaining three hermitian fields correspond to the would-be Goldstone 
bosons that become the longitudinal modes of W ± and Z): 

*m=(S8M* + wv?>) ■ (283) 

[so that (D M 0)tpM0) = l/2(d^H) 2 + ...], with the results 

2 

C[H, W, Z] = g 2 ^=W+W~^H + 9 —W+W-^H 2 + 
v 2 4 

+ g 2 ^^- Z,Z^H+ 9 ——Z,Z^H 2 . (284) 

y 2 v / 2cos 2 ^ 8cos 2 ^ 1 ; 
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Note that the trilinear couplings are nominally of order g 2 , but the adimensional coupling 
constant is actually of order g if we express the couplings in terms of the masses according 
to Eqs.(278): 

a 2 

C[H,W,Z] = gmwW+W-^H + jW+W-^H 2 + 

+ n^Ta- Z ^ H+ 9i 9 lfi Z ^H 2 . (285) 
2 cos 2 W 8 cos 2 v w 

Thus the trilinear couplings of the Higgs to the gauge bosons are also proportional to the 
masses (at fixed g: if instead Gp is kept fixed then, by eq. 244, g is proportional to m w , 
and the Higgs couplings are quadratic in m w ). The quadrilinear couplings are of order g 2 . 
Recall that to go from the lagrangian to the Feynman rules for the vertices the statistical 
factors must be taken into account: for example, the Feynman rule for the ZZHH vertex 
is ig^g 2 /2 cos 2 9w- 

The generic coupling of H to a fermion of type f is given by (after diagonalization): 

C[H,$,^] = ^WH, (286) 

with 

% = = 2^Gfm f . (287) 

The Higgs self couplings are obtained from the potential in Eq.(265) by the replacement 
in Eq.(283). Given that, from the minimum condition: 




,, 2 



(288) 



one obtains: 



V = -n\v +^L) 2 + ^-( v + ^Ly = -^+ f, 2 H 2 + -^Lh 3 + f-H 4 (289) 



v/2 2t; 2V v/2 7 2 ^ y/2v 8v 

The constant term can be omitted in our context. We see that the Higgs mass is positive 
(compare with Eq.(265)) and is given by: 

m 2 H = 2/i 2 = 2\v 2 (290) 

By recalling the value of v in Eq.(279), we see that for m# ~ 126 GeV A is small, A ~ 0.26, 
and the Higgs self interaction is in the perturbative domain. 

The difficulty of the Higgs search is due to the fact that it is heavy and coupled in 
proportion to mass: it is a heavy particle that must be radiated by another heavy particle. 
So a lot of phase space and of luminosity are needed. At LEP2 the main process for 
Higgs production was the Higgs-strahlung process e + e~ — > ZH shown in Fig. 3.5 [219]. 
The alternative process e + e~ — > HuP, via WW fusion, also shown in Fig. 3.5 [220], has 
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a smaller cross-section at LEP2 energies but would become important, even dominant at 
higher energy e + e~ colliders, like the ILC or CLIC (the corresponding ZZ fusion process 
has a much smaller cross-section). The analytic formulae for the cross-sections of both 
processes can be found, for example, in [216]. The direct experimental limit on mn from 
LEP2 was m H > 114 GeV at 95% c.L. The phenomenology of the SM Higgs particle and 
its production and detection at hadron colliders will be discussed in Sects. 3.13, 3.16. 



V 




a b 



Figure 41: Higgs production diagrams in Born approximation for e + e annihilation: (a) The Higgs- 
strahlung process e + e~ — > ZH, (b) the WW fusion process e + e~ — > Hvv. 



3.6 The CKM Matrix 

Weak charged current vertices are the only tree level interactions in the SM that change 
flavour: for example, by emission of a W + an up-type quark is turned into a down-type 
quark, or a V\ neutrino is turned into a l~ charged lepton (all fermions are letf-handed) . 
If we start from an up quark that is a mass eigenstate, emission of a W + turns it into 
a down-type quark state d' (the weak isospin partner of u) that in general is not a mass 
eigenstate. The mass eigenstates and the weak eigenstates do not coincide and a unitary 
transformation connects the two sets: 













M 




1 = VD 


\Vj 









D' = Is' J =V ys j =VD (291) 

where V is the Cabibbo-Kobayashi-Maskawa (CKM) matrix [221]. In analogy with D we 
can denote by U the column vector of the 3 up quark mass eigenstates. Thus in terms of 
mass eigenstates the charged weak current of quarks is of the form: 

J+cx U^(l- l5 )t + VD (292) 

where 

V = UlU d (293) 

Here U u and Ud are the unitary matrices that operate on left-handed doublets in the 
diagonalization of the u and d quarks, respectively (see Eq.(269)). Since V is unitary (i.e. 
yyt = y^y = 1) and commutes with T 2 , T 3 and Q (because all d-type quarks have the 
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same isospin and charge), the neutral current couplings are diagonal both in the primed 
and unprimed basis (if the down-type quark terms in the Z current are written in terms of 
weak isospin eigenvectors as DTD', then by changing basis we get DVTVD and V and 
r commute because, as seen from eq.(248), V is made of Dirac matrices and of T3 and Q 
generator matrices). It follows that DTD' = DTD. This is the GIM mechanism [16] that 
ensures natural flavour conservation of the neutral current couplings at the tree level. 

For N generations of quarks, V is a NxN unitary matrix that depends on N 2 real 
numbers (N 2 complex entries with N 2 unitarity constraints). However, the 2N phases 
of up- and down-type quarks are not observable. Note that an overall phase drops away 
from the expression of the current in eq.(292), so that only 2N — 1 phases can affect V. In 
total, V depends on N 2 — 2N + 1 = (N — l) 2 real physical parameters. A similar counting 
gives N(N — l)/2 as the number of independent parameters in an orthogonal NxN matrix. 
This implies that in V we have N(N — l)/2 mixing angles and (N — l) 2 — N(N — l)/2 = 
(N — 1)(N — 2)/2 phases: for N = 2 one mixing angle (the Cabibbo angle 9c) and no 
phases, for iV = 3 three angles (6*i 2 , #13 and #23) and one phase <p etc. 

Given the experimental near diagonal structure of V a convenient parametrization is the 
one proposed by Maiani [222]. It can be cast in the form of a product of three independent 



2x2 block matrices (stj and C{ j £1X6 shorthands for sin^ and cos 6^): 



V 




-■up 





1 






(294) 



The advantage of this parametrization is that the 3 mixing angles are of different orders of 
magnitude. In fact, from experiment we know that S12 = A, S23 ~ o(A 2 ) and si 3 ~ o(A 3 ), 
where A = sin^c is the sine of the Cabibbo angle, and, as order of magnitude, sy can be 
expressed in terms of small powers of A. More precisely, following Wolfenstein [223] one 
can set: 

s 12 = A, s 23 = AX 2 , si 3 e-^ = A\ 3 (p-t V ) (295) 



o(A 4 ). (296) 



As a result, by ne£ 


dectinj 


7 terms of higher order in 


A 


one 


can write down 




V u d v us 


v ub - 




r 1-- 
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AX 3 (p — it]) ' 


V = 


V c d V cs 


v cb 




-A 
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A 2 
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AX 2 




V td V ts 


v tb _ 




_AX 3 (1 - p-irj) 




AX 2 


1 



It has become customary to make the replacement p, rj — > p, 7] with: 

/'-"/= ^7== ~ (p-^)(l + A 2 /2 + ...). 



p — 17] 

VT^A 1 



(297) 



The best values of the CKM parameters as obtained from experiment are continuously 
updated in refs. [224,225] (a survey of the current status of the CKM parameters can also 
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be found in ref. [14]). A Summer 2012 fit [224] led to the values (compatible values, within 
stated errors, are given in ref. [225]): 

A = 0.2254 ± 0.0006 
A = 0.824 ±0.013 

p = 0.139 ±0.021; fj = 0.352 ± 0.016 (298) 
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Figure 42: The unitarity triangle corresponding to eq.(299). 

In the SM the non vanishing of the fj parameter (related to the phase <p in Eqs. 294 and 
295) is the only source of CP violation in the quark sector (we shall see that new sources 
of CP violation very likely arise from the neutrino sector). Unitarity of the CKM matrix 
V implies relations of the form Vb a V* a = 5& c . In most cases these relations do not imply 
particularly instructive constraints on the Wolfenstein parameters. But when the three 
terms in the sum are of comparable magnitude we get interesting information. The three 
numbers which must add to zero form a closed triangle in the complex plane (unitarity 
triangle), with sides of comparable length. This is the case for the t-u triangle shown in 
Fig. 42 (or, what is equivalent in first approximation, for the d-b triangle): 

V td V: d + V ts V; s + V tb V: b = (299) 

All terms are of order A 3 . For rj=0 the triangle would flatten down to vanishing area. 
In fact the area of the triangle, J of order J ~ rjA 2 X 6 , is the Jarlskog invariant [227] 
(its value is independent of the parametrization). In the SM, in the quark sector, all CP 
violating observables must be proportional to J, hence to the area of the triangle or to 
f]. Its experimental value is J ~ (3.2 ± 0.1) 10~ 5 [224]. A direct and by now very solid 
evidence for J non vanishing has been first obtained from the measurements of e and e' in 
K decay. Additional direct evidence has more recently been collected from experiments on 
B decays at beauty factories, at the Tevatron and at the LHC (in particular by the LHCb 
experiment). Very recently some evidence for CP violation in D decays has also been 
reported by the LHCb experiment [226]. The angles (3 (the most precisely measured), a 
and 7 have been determined with fair precision. The angle measurements and the available 
information on the magnitude of the sides, taken together, are in good agreement with the 
predictions from the SM unitary triangle (see Fig. 43) [224] [225]. Some alleged tensions 
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Figure 43: Constraints in the p, rj plane including the most recent data inputs (summer 2012) in the global 



are not convincing either because of their poor statistical significance or because of lack 
of confirmation from different potentially sensitive experiments or because the associated 
theoretical error estimates can be questioned. 

As we have discussed, due to the GIM mechanism, there are no flavour changing neutral 
current (FCNC) transitions at the tree level in the SM. Transitions with |AF| = 1,2 
are induced at one loop level. In particular, meson mixing, i.e. M — > M off diagonal 
\AF\ = 2 mass matrix elements (with M = K, D or B neutral mesons), are obtained 
from box diagrams. For example, in the case of K° — K° mixing the relevant transition 
is sd — > sd (see Fig. 44). In the internal quark lines all up-type quarks are exchanged. In 
the amplitude, two vertices and the connecting propagator (with virtual four momentum 
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Figure 44: Box diagrams describing K° — K° mixing at the quark level at 1-loop. 
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Pfj) at one side contribute a factor (w, = u, c,t): 

Fgim = Yl K.jT^— V ** ' ( 30 °) 
*-r< p-m ui 

which, in the limit of equal m u i, is clearly vanishing due to the unitarity of the CKM matrix 
V. Thus the result is proportional to mass differences. For K° — K° mixing the contribution 
of virtual u quarks is negligible due to the small value of m u and the contribution of the 
t quark is also small due to the mixing factors V t *V t d ~ o(A 2 X 5 ). The dominant c quark 
contribution to the real part of the box diagram quark-level amplitude is approximately of 
the form (see, for example, [228]): 

ReH box = ^ m 2ife(V£Krf)VO Ai=2 , (301) 

l07T Z 

where r\\ ~ 0.85 is a QCD correction factor and As=2 = (IlI^sl sij^di is the relevant 
4-quark, dimension-6, operator. The r)i factor arises from gluon exchanges among the 
quark legs of the 4-quark operator. Indeed the coefficients of the operator expansion 
that arises when the heavy particles exchanged are integrated away obey renormalization 
group equations and the associated logarithms can be resummed (the first calculation of 
resummed QCD corrections to weak non leptonic amplitudes was performed in refs. [229] ; 
for a pedagogical introduction see, for example, ref. [230]). To obtain the K° — K° mixing 
amplitude the matrix element between meson states must be taken of As=2 which is 
parametrized in terms of a "Br parameter" which is defined in such a way that B% = 1 
for vacuum state insertion between the two currents: 

(K°\0 As = 2 \K°) = ^f K m 2 K B K , (302) 

where B K ~ 0.75 (this is the renormalization group independent definition usually denoted 
as Bk) and fx ~ 113 MeV, the kaon pseudoscalar constant, are best evaluated by QCD 
lattice simulations [231]. Clearly to the charm parton contribution in eq.(301) additional 
non perturbative terms must be added, some of them of o(m 2 K /ml), because the smallness of 
m c makes a completely partonic dominance inadequate. In Eq.(301) the factor o(m 2 /m 2 v ) 
is the "GIM suppression" factor (1/m^ is hidden in Gp according to eq. (244)). 

For B mixing the dominant contribution is from the t quark. In this case, the partonic 
dominance is more realistic and the GIM factor o{m 2 /m 2 v ) is actually larger than 1. More 
recently also D mixing has been observed [232] . In the corresponding box diagrams down- 
type quarks are involved. But starting from D ~ cu the b-quark contribution is strongly 
suppressed by the CKM angles (given that V c bV* b ~ o(Xq). The masses of the d and s 
quarks are too small for a partonic evaluation of the box diagram and non perturbative 
terms cannot be neglected. This makes a theoretical evaluation of mixing and CP violation 
effects for D mesons problematic. 

All sorts of transitions with |AF| = 1 are also induced at loop level. For example, 
an effective vertex Z — >■ tc, which does not exist at tree level, is generated at 1-loop (see 
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a b c 

Figure 45: Examples of \AF\ = 1 transitions at the quark level at 1-loop: (a) Diagram for a Z — s- t c 
vertex, (b) b — > s 7, (c) a "penguin" diagram for b — > s e + e~ . 

Fig. 45). Similarly, transitions involving photons or gluons are also possible, like t — > c g or 
b — > s 7 (Fig. 45) or b — )■ s g. For light fermion exchange in the loop the GIM suppression 
is also effective in \AF\ = 1 amplitudes. For example, analogous leptonic transitions like 
fi->e7orr4|i7 also exist but in the SM are extremely small and out of reach for 
experiments, because the tiny neutrino masses enter in the GIM suppression factor. But 
new physics effects could well make these rare lepton flavour violating processes accessible 
to experiment. In fact, the present limits already pose stringent constraints on models of 
new physics. Of particular importance is the recent bound [233] obtained by the MEG 
Collaboration at SIN, near Zurich, Switzerland, on the branching ratio for fi — > e 7: 
B(fi — > e 7) < 5.7 1CT 13 at 90%. The external Z, photon or gluon can be attached to 
a pair of light fermions, giving rise to an effective four fermion operator, as in "penguin 
diagrams" like the one shown in Fig. 45 for b — > s l + l~ . The inclusive rate B — > X s 7 (here 
B stands for Bj) with X s a hadronic state containing a unit of strangeness corresponding 
to an s-quark, has been precisely measured. The world average result for the branching 
ratio with E 1 > 1.6 GeV is [234]: 

B(B -> X s -f) exp = (3.55 ± 0.26)-10" 4 . (303) 

The theoretical prediction for this inclusive process is to a large extent free of uncertainties 
from hadronization effects and is accessible to perturbation theory as the b-quark is heavy 
enough. The most complete result at order a 2 s is at present from ref. [235] and refs. therein: 

B(B -> X s 7 ) tfc = (2.98 ± 0.26)10" 4 . (304) 

Note that the theoretical value has recently become smaller than the experimental value. 
The fair agreement between theory and experiment imposes stringent constraints on pos- 
sible new physics effects. 

Related processes are B s ^ — > n + fi~. These decay are very rare in the SM, their pre- 
dicted branching ratio being B(B S -»■ /i + /i~) ~ (3.23 ± 0.27) 10~ 9 , B(B d ~ 
(1.07 ± 0.10) 10~ 10 [236]. These very small expected branching ratios result from being 
FCNC processes with helicity suppression in the purely leptonic final state (the decaying 
meson has spin zero and the muon pair is produced by vector exchange in the SM). Many 
models of new physics beyond the SM predict large deviations. Thus these processes pose 
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very stringent tests to the SM. Recently the LHCb experiment has reached the sensitiv- 
ity to observe the B s mode and has published the important result B(B S — > n + n~) = 
3.2^2' 5 10~ 9 [237], in very good agreement with the SM (in the same paper the bound 
B(Bd — > < 9.4 1CT 10 at 95% c.l. was quoted). Another very demanding test of the 

SM has been passed! 

In conclusion, the CKM theory of quark mixing and CP violation has been precisely 
tested in the last decade and turns out to be very successful. The expected deviations 
from new physics at the EW scale did not appear so far. Either the scale of new physics 
is very far, in the 10 3 — 10 5 TeV range, or its flavour structure must be very special. For 
example, the class of models with minimal flavour violation (MFV) [238], where the SM 
Yukawa couplings are the only flavour symmetry breaking terms also beyond the SM, have 
been much studied and represent a sort of limiting baseline. 

3.7 Neutrino Mass and Mixing 

In the minimal version of the SM the right handed neutrinos which have no gauge 
interactions, are not present at all. With no vr no Dirac mass is possible for neutrinos. If 
lepton number conservation is also imposed, then no Majorana mass is allowed either and, 
as a consequence, all neutrinos are massless. But, at present, from neutrino oscillation 
experiments, we know that at least 2 out of the 3 known neutrinos have non vanishing 
masses (for reviews, see, for example, refs. [239]: the two mass squared differences measured 
from solar (Am^) and atmospheric oscillations (Amy are given by AmJ 2 ~ 8 10~ 5 eV 2 
and Am^ ~ 2.5 1CT 3 eV 2 [240-242]. The absolute values of the masses are very small, 
with an upper limit of a fraction of eV, obtained from laboratory experiments (tritium 
j3 decay near the end point: m v < 2 eV [14], absence of visible neutrinoless double (3 
decay : \m ee \ < 0.3 — 0.7 eV (m ee is a combination of neutrino masses; for a review, see, 
for example [243]) and from cosmological observations: m u < 0.1 — 0.7 eV (depending 
on the cosmological model assumptions) [244]. If um are added to the minimal model 
and lepton number is imposed by hand, then neutrino masses would in general appear 
as Dirac masses, generated by the Higgs mechanism, like for any other fermion. But, for 
Dirac neutrinos, to explain the extreme smallness of neutrino masses, one should allow for 
very small Yukawa couplings. However, we stress that, in the SM, baryon B and lepton 
L number conservation, which are not guaranteed by gauge symmetries (as is the case 
for the electric charge Q), are understood as "accidental" symmetries. In fact the SM 
lagrangian should contain all terms allowed by gauge symmetry and renormalizability, but 
the most general renormalizable lagrangian (i.e. consisting of operator dimension d < 4), 
built from the SM fields, compatible with the SM gauge symmetry, in the absence of vm, is 
automatically B and L conserving (however, non perturbative instanton effects break the 
conservation of B+L while preserving B-L, as discussed in Sect. 3.8). In the presence of 
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m. 



(306) 



Urn, this is no more true and the right handed Majorana mass term is allowed: 

Mrr = v c lR M i3 Vj R = vfxCMijVjR , (305) 

where uf R = CvJ R is the charge conjugated neutrino field and C is the charge conjugation 
matrix in Dirac spinor space. The Majorana mass term is an operator of dimension d = 3 
with AL = 2. Since the vm are gauge singlets the Majorana mass Mrr is fully allowed by 
the gauge symmetry and a coupling with the Higgs is not needed to generate this type of 
mass. As a consequence, the entries of the mass matrix M^- do not need to be of the order 
of the EW symmetry breaking scale v and could be much larger. If one starts from the 
Dirac and RR Majorana mass terms for neutrinos, the resulting mass matrix, in the L, R 
space, has the form: 

mo 
m D M 

where mo and M are the Dirac and Majorana mass matrices (M is the matrix in 
eq.(305)). The corresponding eigenvalues are 3 very heavy neutrinos with masses of order 
M and 3 light neutrinos with masses 

m v = - m T D M- l m D , (307) 

which are possibly very small if M is large enough. This is the see-saw mechanism for 
neutrino masses [245]. Note that if no um exist a Majorana mass term could still be built 
out of VjL. But VjL have weak isospin 1/2, being part of the left handed lepton doublet I. 
Thus, the left handed Majorana mass term has total weak isospin equal to 1 and needs 2 
Higgs fields to make a gauge invariant term. The resulting mass term: 

5 = XlJXijljHH/M , (308) 

with M a large scale (apriori comparable to the scale of Mrr) and A a dimensionless 
coupling generically of o(l), is a non renormalizable operator of dimension 5, first pointed 
out by S. Weinberg [246]. The corresponding mass terms are of the order m v ~ Xv 2 /M, 
hence of the same generic order of the light neutrino masses from eq.(307). 

In conclusion, neutrino masses are believed to be small because neutrinos are Majorana 
particles with masses inversely proportional to the large scale M of energy where L non 
conservation is induced. This corresponds to an important enlargement of the original 
minimal SM where no vr was included and L conservation was imposed by hand (actually 
it is true "accidentally" only at the renormalisable level, but is violated by the non renor- 
malisable Weinberg operator and by instanton effects). It is interesting that the observed 
magnitudes of the mass squared splittings of neutrinos are well compatible with a scale 
M remarkably close to the Grand Unification scale, where indeed L non conservation is 
naturally expected. In fact, for m v ~ a/ AmJ tm ~ 0.05 eV (see Table(3)) and m u w m 2 D /M 
with mo ~ v ~ 200 GeV we find M ~ 10 15 GeV which indeed is an impressive indication 
for M GUT . 
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In the previous Section we have discussed flavour mixing for quarks. But, clearly, 
given that non vanishing neutrino masses have been established, a similar mixing matrix 
is also introduced in the leptonic sector. We assume in the following that there are only 
two distinct neutrino oscillation frequencies, the atmospheric and the solar frequencies 
(both of them now also confirmed by experiments where neutrinos are generated on the 
earth like K2K, KamLAND, MINOS). At present the bulk of neutrino oscillation data 
are well reproduced in terms of three light neutrino species. However, some (so far not 
compelling) evidence for additional "sterile" neutrino species (i.e. not coupled to the weak 
interactions, as demanded by the LEP limit on the number of "active" neutrinos) are 
present in some data. We discuss here 3-neutrino mixing, which is in any case a good 
approximate framework to discuss neutrino oscillations, while for possible sterile neutrinos 
we refer to the comprehensive review in ref. [247]. 

Neutrino oscillations are due to a misalignment between the flavour basis, v ' = [y e , z/ M , u T ), 
where v e is the partner of the mass and flavour eigenstate e~ in a left-handed (LH) weak 
isospin SU(2) doublet (similarly for and z/ r )) and the mass eigenstates v = (^1,^2,^3) 
[239,248,249]: 

v' = Uv , (309) 

where U is the unitary 3 by 3 mixing matrix. Given the definition of U and the transfor- 
mation properties of the effective light neutrino mass matrix m„: 

u' T m v u' = v T U T m v Uv (310) 
U T m v U = Diag (mi, m 2 , m 3 ) = m diag , 

we obtain the general form of m v (i.e. of the light v mass matrix in the basis where the 
charged lepton mass is a diagonal matrix): 

m v = U*m diag U ] . (311) 

The matrix U can be parameterized in terms of three mixing angles #12, #23 and #13 (0 < 
&ij < 7r /2) and one phase ip (0 < tp < 27r) [250], exactly as for the quark mixing matrix 
Vckm- The following definition of mixing angles can be adopted: 

/I \ / C13 si 3 e*n / Cia S12 0\ 
U = C 23 S 23 1 -8X2 C 12 (312) 

\0 -s 23 c 23 / \-s 13 e-^ cia / V l) 

where Sy = sin^-, Cy = cos 9ij. In addition, if v are Majorana particles, we have two 
more phases [251] given by the relative phases among the Majorana masses mi, m 2 and 
m 3 . If we choose m 3 real and positive, these phases are carried by m 12 = (wi^le 1 ^ 1 ' 2 . 
Thus, in general, 9 parameters are added to the SM when non-vanishing neutrino masses 
are included: 3 eigenvalues, 3 mixing angles and 3 CP violating phases. 



109 



In our notation the two frequencies, Am'j/AE (I = sun,atm), are parametrized in 
terms of the v mass eigenvalues by 

Am lun = \ Am u\, Am ltm = \ Am ls\ ■ ( 313 ) 

where Amf 2 = \rri2\ 2 — \mi\ 2 > and Am^ = m\ — \m2\ 2 . The numbering 1,2,3 corresponds 
to a definition of the frequencies and in principle may not coincide with the ordering from 
the lightest to the heaviest state. "Normal hierarchy" is the case where 777,3 is the largest 
mass in absolute value, otherwise one has an " inverse hierarchy" . 

Very important developments in the data have occurred in 2012. The value of the 
mixing angle #13 has been proven to be non vanishing and its value is by now known with 
fair precision. Several experiments have been involved in the #13 measurement and their 
results are reported in table 2. 



Quantity 


sin 2 2# 13 


sin 2 6>i3 


T2K [252] 
MINOS [253] 
DC [254] 

DYB [255] 
RENO [256] 


o.n±8S (o.i4i&S) 

0.041+^ (0.079^3) 
0.086 ±0.041 ±0.030 
0.092 ±0.016 ±0.005 
0.113 ±0.013 ±0.019 


0.0281^ (0.036lffi) 

o.oioiS (o.o2o«:P 
o.o22+8;8^ 

0.024 ±0.005 
0.029 ±0.006 



Table 2: The reactor angle measurements from the recent experiments T2K [252] , MINOS [253] , DOUBLE 
CHOOZ [254], Daya Bay [255] and RENO [256], for the normal (inverse) hierarchy. 



The present data are summarized in table 3 [240]. 



Am 2 un (10- 5 eV 2 ) 


7.5418;! 


Aml tm (10- 3 eV 2 ) 


2.431H8 (2.4218;?!) 


sin 2 #12 


0.30718:81^ 


sin 2 #23 


o.386i8:8^ (o.392i8;81) 


sin 2 6»i3 


o.o24ii8;88l (0.024418:881) 


(p/n 


1.08I8I? (1.0918:1) 



Table 3: Fits to neutrino oscillation data from Ref. [240]. The results for both the normal and the inverse 
(in the brackets) hierarchies are shown. 

Neutrino mixing is important because it could in principle provide new clues for the 
understanding of the flavour problem. Even more so since neutrino mixing angles show a 
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pattern that is completely different than that of quark mixing: for quarks all mixing angles 
are small, for neutrinos two angles are large (one is still compatible with the maximal value) 
and only the third one is small. In reality it is frustrating that no real illumination was 
sparked on the problem of flavour. We can reproduce in models the data on neutrino 
mixing, in a wide range of dynamical setups that goes from anarchy to discrete flavour 
symmetries (for reviews and references see, for example refs. [257-262]) but we have not 
yet been able to single out a unique and convincing baseline for the understanding of 
fermion masses and mixings. In spite of many interesting ideas and the formulation of 
many elegant models the mysteries of the flavour structure of the three generations of 
fermions have not been much unveiled. 



3.8 Quantization and renormalization of the Electroweak Theory 

The Higgs mechanism gives masses to the Z, the W and to fermions while the la- 
grangian density is still symmetric. In particular the gauge Ward identities and the sym- 
metric form of the gauge currents are preserved. The validity of these relations is an 
essential ingredient for renormalizability. In the previous Sections we have specified the 
Feynman vertices in the "unitary" gauge where only physical particles appear. However, 
as discussed in Chapter 1, in this gauge the massive gauge boson propagator would have a 
bad ultraviolet behaviour: g ^ 

WV = — (314) 

<f-rnfo 

A formulation of the standard EW theory with good apparent ultraviolet behaviour can be 
obtained by introducing the renormalizable or gauges [21], in analogy with the abelian 
case discussed in detail in Chapter 1. One parametrizes the Higgs doublet as: 



^3 + 204 / V V + 



V2 



(315) 



and similarly for <ffl, where w~ appears. The scalar fields and z are the pseudo Goldstone 
bosons associated with the longitudinal modes of the physical vector bosons W and Z. 
The i?£ gauge fixing lagrangian has the form: 

AC GF = - l -\d»W, - im w w\ 2 - ^(0% - Vmzzf - ^(^AJ 2 ■ (316) 

The W and Z propagators, as well as those of the scalars w and z, have exactly the same 
general forms as for the abelian case in eqs. (67)-(69) of Chapter 1, with parameters £ and 
?7, respectively (and the pseudo Goldstone bosons and z have masses ^my/ and rjmz)- 
In general, a set of associated ghost fields must be added, again in direct analogy with the 
treatment of gauges in the abelian case of Chapter 1. The complete Feynman rules for 
the standard EW theory can be found in a number of textbooks (see, for example, [263]). 
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The pseudo Goldstone bosons w ± and z are directly related to the longitudinal helic- 
ity states of the corresponding massive vector bosons W and Z. This correspondence 
materializes in a very interesting "equivalence theorem": at high energies of order E the 
amplitude for the emission of one or more longitudinal gauge bosons Vl (with V = W, Z) 
becomes equal (apart from terms down by powers of my / E) to the amplitude where each 
longitudinal gauge boson is replaced by the corresponding Goldstone field w ± or z [264]. 
For example, consider top decay with a longitudinal W in the final state: t — > bW^. The 
equivalence theorem asserts that we can compute the dominant contribution to this rate 
from the simpler t — > bw + matrix element: 

r(t ->■ bW£) = T(t -> bw + )[l + o(m 2 w /m 2 t )} . (317) 

In fact one finds: 

T(t bw+) = JU = f^f , (318) 
62ir 87TV2 

where h t = m t /v is the Yukawa coupling of the top quark (numerically very close to 1), 
and we used 1/v 2 = 2y/2Gp (see Eq.(279)). If we compare with Eq.(259), we see that this 
expression coincides with the total top width (i.e. including all polarizations for the W in 
the final state), computed at tree level, apart from terms down by powers of o{m 2 v /rn 2 ). 
In fact, the longitudinal W is dominant in the final state because hf » g 2 . Similarly 
the equivalence theorem can be applied to find the dominant terms at large yfs for the 
cross-section e + e" — >■ W£W£ , or the leading contribution, in the limit m H » m v , to the 
width for the decay T(H -> VV). 

The formalism of the gauges is also very useful in proving that spontaneously broken 
gauge theories are renormalizable. In fact, the non singular behaviour of propagators at 
large momenta is very suggestive of the result. Nevertheless it is by far not a simple matter 
to prove this statement. The fundamental theorem that in general a gauge theory with 
spontaneous symmetry breaking and the Higgs mechanism is renormalizable was proven 
by 't Hooft and Veltman [27], [28]. 

For a chiral theory like the SM an additional complication arises from the existence of 
chiral anomalies. But this problem is avoided in the SM because the quantum numbers of 
the quarks and leptons in each generation imply a remarkable (and, from the point of view 
of the SM, mysterious) cancellation of the anomaly, as originally observed in ref. [265]. 
In quantum field theory one encounters an anomaly when a symmetry of the classical 
lagrangian is broken by the process of quantization, regularization and renormalization of 
the theory. Of direct relevance for the EW theory is the Adler-Bell-Jackiw (ABJ) chiral 
anomaly [47]. The classical lagrangian of a theory with massless fermions is invariant 
under a U(l) chiral transformations ipf = e 1150 ^. The associated axial Noether current is 
conserved at the classical level. But, at the quantum level, chiral symmetry is broken due 
to the ABJ anomaly and the current is not conserved. The chiral breaking is produced by 
a clash between chiral symmetry, gauge invariance and the regularization procedure. 
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Figure 46: Triangle diagram that generates the ABJ anomaly [47]. 



The anomaly is generated by triangular fermion loops with one axial and two vector 
vertices (Fig. 46). For example, for the Z the axial coupling is proportional to the 3rd com- 
ponent of weak isospin t 3 , while the vector coupling is proportional to a linear combination 
of t 3 and the electric charge Q. Thus in order for the chiral anomaly to vanish all traces 
of the form tr{t 3 QQ}, tr{t 3 t 3 Q}, tr{t 3 t 3 t 3 } (and also tr{t + t_t 3 } when charged currents 
are also included) must vanish, where the trace is extended over all fermions in the theory 
that can circulate in the loop. Now all these traces happen to vanish for each fermion 
family separately. For example take tr{t 3 QQ}. In one family there are, with £3 = +1/2, 
three colours of up quarks with charge Q = +2/3 and one neutrino with Q = and, with 
t 3 = —1/2, three colours of down quarks with charge Q = —1/3 and one l~ with Q — — 1. 
Thus we obtain tr{t 3 QQ} = 1/2-3-4/9 - 1/2-31/9 - 1/21 = 0. This impressive can- 
cellation suggests an interplay among weak isospin, charge and colour quantum numbers 
which appears as a miracle from the point of view of the low energy theory but is in fact 
understandable from the point of view of the high energy theory. For example, in Grand 
Unified Theories (GUTs) (for reviews, see, for example, [266]) there are similar relations 
where charge quantization and colour are related: in the 5 of SU(5) we have the content 
(d, d, d, e + , v) and the charge generator has a vanishing trace in each SU(5) representation 
(the condition of unit determinant, represented by the letter S in the SU(5) group name, 
translates into zero trace for the generators). Thus the charge of d quarks is -1/3 of the 
positron charge because there are three colours. A whole family fits perfectly in one 16 
of SO(10) which is anomaly free. So GUTs can naturally explain the cancellation of the 
chiral anomaly. 

An important implication of chiral anomalies together with the topological properties of 
the vacuum in non abelian gauge theories is that the conservation of the charges associated 
to baryon (B) and lepton (L) numbers is broken by the anomaly [267], so that B and L 
conservation is actually violated in the standard electroweak theory (but B-L remains 
conserved). B and L are conserved to all orders in the perturbative expansion but the 
violation occurs via non perturbative instanton effects [268] (the amplitude is proportional 
to the typical non perturbative factor exp —c/g 2 , with c a constant and g the 577(2) gauge 
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coupling). The corresponding effect is totally negligible at zero temperature T, but becomes 
relevant at temperatures close to the electroweak symmetry breaking scale, precisely at 
T ~ o(TeV). The non conservation of B+L and the conservation of B-L near the weak 
scale plays a role in the theory of baryogenesis that quantitatively aims at explaining the 
observed matter antimatter asymmetry in the Universe (for reviews and references, see, 
for example, [269]). 

3.9 QED Tests: Lepton Anomalous Magnetic Moments 

The most precise tests of the electroweak theory apply to the QED sector. Here we 
discuss the anomalous magnetic moments of the electron and of the muon that are among 
the most precise measurements in the whole of physics. The magnetic moment jl and 
the spin S are related by jl = —geS /2m, where g is the gyromagnetic ratio (g = 2 for a 
pointlike Dirac particle). The quantity a = (g — 2)/2 measures the anomalous magnetic 
moment of the particle. Recently there have been new precise measurements of a e and a M 
for the electron [270] and the muon [271]: 

a e e xp = 11596521807.3(2.8) 10~ 13 , a e * p = 11659208.9(6.3) lO" 10 . (319) 

The theoretical calculations in general contain a pure QED part plus the sum of hadronic 
and weak contribution terms: 

a = a QED + a hadTonic + a weak = V C i (-) i + a hadronic + a weak . (320) 

i 

The QED part has been computed analytically for % = 1,2,3, while for % = 4 there is a 
numerical calculation with an error (see, for example, [273] and refs therein). The complete 
numerical evaluation of i = 5 for the muon case has been published in 2012 [274] as a new 
impressive achievement of Kinoshita and his group. The hadronic contribution is from 
vacuum polarization insertions and from light by light scattering diagrams (see Fig. 47). 
The weak contribution is from W or Z exchange. 



Figure 47: The hadronic contributions to the anomalous magnetic moment: vacuum polarization (left) 
and light by light scattering (right). 

For the electron case the weak contribution is essentially negligible and the hadronic 
term ^ a hadromc ^ (16.82 ± 0.19)10 -13 ) does not introduce an important uncertainty. As a 
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Figure 48: Compilation of recently published results for a M (in units of 10 -11 ) [279]: JN [280], DHMZ [276], 
HLMNT [277]. 

result this measurement can be used to obtain the most precise determination of the fine 
structure constant [274]: 



a 



137.0359991657(340) 



(321) 



In the muon case the experimental precision is less by about 3 orders of magnitude, but 
the sensitivity to new physics effects is typically increased by a factor (m M /m e ) 2 ~ 410 4 
(one mass factor arises because the effective operator needs a chirality flip and the second 
one is because, by definition, one must factor out the Bohr magneton e/2m). From the 
theory side, the QED term (using the value of a from a e in Eq.(321)), and the weak 
contribution [275] are affected by small errors and are given by (all theory numbers given 
here are taken from ref. [274]) 



oed = (H6584718.853 ± 0.037) 10' 



ii 



.uieafc 



(154 ±2.0)10 



-ii 



(322) 



The dominant ambiguities arise from the hadronic term. The lowest order (LO) vacuum 
polarization contribution can be evaluated from the measured cross sections in e + e~ — > 
hadrons at low energy via dispersion relations (the largest contribution is from the tttt final 
state) [276,277], with the result a^°10" n = 6949 ± 43. The higher order (HO) vacuum 
polarization contribution (from 2-loop diagrams containing an hadronic insertion) is given 
by: aj°10- n = -98.4 ±0.7 [277]. The contribution of the light by light (LbL) scattering 
diagrams is estimated to be: a^ bL 10~ n = 116 ±40 [278]. Adding the above contributions 
up the total hadronic result is reported as: 



hadronic 
> 



(6967 ±59)10 



-ii 



(323) 
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At face value this would lead to a 2.9a deviation from the experimental value a e * p in 
Eq.(319): 

a exp _ a th(e+e-) = ± 87 )- 10 -H. ( 324 ) 

However, the error estimate on the LbL term, mainly a theoretical uncertainty, is not 
compelling, and it could well be somewhat larger (although probably not by as much 
as to make the discrepancy to completely disappear). A minor puzzle is the fact that, 
using the conservation of the vector current (CVC) and isospin invariance, which are well 
established tools at low energy, can also be evaluated from r decays. But the results on 
the hadronic contribution from e + e~ and from r decay, nominally of comparable accuracy, 
are still somewhat apart (although the two are now closer than in the past), and the (g-2) 
discrepancy would be attenuated if one takes the r result (see Fig. 48 which refers to 
the most recent results). Since it is difficult to find a theoretical reason for the e + e~ vs 
t difference, one must conclude that perhaps there is something which is not understood 
either in the data or in the assessment of theoretical errors. The prevailing view is to take 
the e + e~ determination as the most directly reliable, which leads to Eq.(324), but some 
doubts remain. Finally, we note that, given the great accuracy of the a M measurement and 
the relative importance of the non QED contributions, it is not unreasonable that a first 
signal of new physics would appear in this quantity. 

3.10 Large Radiative Corrections to Electroweak Processes 

Since the SM theory is renormalizable higher order perturbative corrections can be 
reliably computed. Radiative corrections are very important for precision EW tests. The 
SM inherits all successes of the old V-A theory of charged currents and of QED. Modern 
tests have focussed on neutral current processes, the W mass and the measurement of 
triple gauge vertices. For Z physics and the W mass the state of the art computation 
of radiative corrections include the complete one loop diagrams and selected dominant 
multi loop corrections. In addition some resummation techniques are also implemented, 
like Dyson resummation of vacuum polarization functions and important renormalization 
group improvements for large QED and QCD logarithms. We now discuss in more detail 
sets of large radiative corrections which are particularly significant (for reviews of radiative 
corrections for LEP1 physics, see, for example: [281]; for a more pedagogical description of 
LEP physics see [282]). 

Even leaving aside QCD corrections, a set of important quantitative contributions to the 
radiative corrections arise from large logarithms [e.g. terms of the form (a/n In (m^/m/^)) n 
where fu is a light fermion]. The sequences of leading and close-to- leading logarithms 
are fixed by well-known and consolidated techniques (/3 functions, anomalous dimensions, 
penguin-like diagrams, etc.). For example, large logarithms from pure QED effects domi- 
nate the running of a from m e , the electron mass, up to mz- Similarly large logarithms of 
the form [a/n In (mz//j)] n also enter, for example, in the relation between sin 2 9w at the 
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scales m z (LEP, SLC) and \i (e.g. the scale of low-energy neutral-current experiments). 
Also, large logs from initial state radiation dramatically distort the line shape of the Z 
resonance as observed at LEP1 and SLC and this effect was accurately taken into account 
for the measurement of the Z mass and total width. The experimental accuracy on mz 
obtained at LEP1 is 5m z = ±2.1 MeV. Similarly, a measurement of the total width to an 
accuracy 5T = ±2.3 MeV has been achieved. The prediction of the Z line-shape in the SM 
to such an accuracy has posed a formidable challenge to theory, which has been successfully 
met. For the inclusive process e + e~ — > ffX, with / ^ e (for a concise discussion, we leave 
Bhabha scattering aside) and X including 7's and gluons, the physical cross-section can 
be written in the form of a convolution [281]: 

<t(s)= I dz a(zs)G(z,s) , (325) 

Jzo 

where & is the reduced cross-section, and G(z, s) is the radiator function that describes 
the effect of initial-state radiation; a includes the purely weak corrections, the effect of 
final-state radiation (of both 7's and gluons), and also non-factorizable terms (initial- and 
final-state radiation interferences, boxes, etc.) which, being small, can be treated in lowest 
order and effectively absorbed in a modified a. The radiator G(z, s) has an expansion of 
the form 

G(z,s) = 5(1 - z) ± a/ir(a 11 L ± a w ) ± (a / n) 2 (a 2 2L 2 ± a xx L + a 2 o) ± ... ± 

n 

+ (a/n) n Ys a mL l , (326) 

i=0 

where L = In s/m 2 ~ 24.2 for ~ mz- All first- and second-order terms are known 
exactly. The sequence of leading and next-to-leading logs can be exponentiated (closely 
following the formalism of structure functions in QCD). For mz ~ 91 GeV, the convolution 
displaces the peak by ±110 MeV, and reduces it by a factor of about 0.74. The exponen- 
tiation is important in that it amounts to an additional shift of about 14 MeV in the peak 
position with respect to the 1 loop radiative correction. 

Among the one loop EW radiative corrections, a very remarkable class of contributions 
are those terms that increase quadratically with the top mass. The sensitivity of radiative 
corrections to m t arises from the existence of these terms. The quadratic dependence 
on m t (and on other possible widely broken isospin multiplets from new physics) arises 
because, in spontaneously broken gauge theories, heavy virtual particles do not decouple. 
On the contrary, in QED or QCD, the running of a and a s at a scale Q is not affected 
by heavy quarks with mass M ^> Q. According to an intuitive decoupling theorem [283], 
diagrams with heavy virtual particles of mass M can be ignored at Q C M provided 
that the couplings do not grow with M and that the theory with no heavy particles is 
still renormalizable. In the spontaneously broken EW gauge theories both requirements 
are violated. First, one important difference with respect to unbroken gauge theories is in 
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the longitudinal modes of weak gauge bosons. These modes are generated by the Higgs 
mechanism, and their couplings grow with masses (as is also the case for the physical 
Higgs couplings). Second, the theory without the top quark is no more renormalizable 
because the gauge symmetry is broken as the (t,b) doublet would not be complete (also 
the chiral anomaly would not be completely cancelled). With the observed value of m t 
the quantitative importance of the terms of order G pw? t / An 2 \/2 is substantial but not 
dominant (they are enhanced by a factor mj/m^ ~ 5 with respect to ordinary terms). 
Both the large logarithms and the Gftu 2 terms have a simple structure and are to a large 
extent universal, i.e. common to a wide class of processes. In particular the Gpm\ terms 
appear in vacuum polarization diagrams which are universal (virtual loops inserted in 
gauge boson internal lines are independent of the nature of the vertices on each side of the 
propagator) and in the Z — > bb vertex which is not. This vertex is specifically sensitive 
to the top quark which, being the partner of the b quark in a doublet, runs in the loop. 
Instead all types of heavy particles could in principle contribute to vacuum polarization 
diagrams. The study of universal vacuum polarization contributions, also called "oblique" 
corrections, and of top enhanced terms is important for an understanding of the pattern 
of radiative corrections. More in general, the important consequence of non decoupling is 
that precision tests of the electroweak theory may apriori be sensitive to new physics even 
if the new particles are too heavy for their direct production, but aposteriori no signal of 
deviation has clearly emerged. 

While radiative corrections are quite sensitive to the top mass, they are unfortunately 
much less dependent on the Higgs mass. In fact, the dependence of one loop diagrams on 
run is only logarithmic: ~ Gftv^ \og{m 2 H / m^y) . Quadratic terms ~ G 2 F m 2 H only appear 
at two loops [284] and are too small to be detectable. The difference with the top case 
is that the splitting m 2 — m 2 is a direct breaking of the gauge symmetry that already 
affects the 1- loop corrections, while the Higgs couplings are "custodial" SU(2) symmetric 
in lowest order. 

3.11 Electroweak Precision Tests 

For the analysis of electroweak data in the SM one starts from the input parameters: 
as is the case in any renormalizable theory, masses and couplings have to be specified from 
outside. One can trade one parameter for another and this freedom is used to select the 
best measured ones as input parameters. Some of them, a, Gf and mz, are very precisely 
known, as we have seen, some other ones, nif Ught , mt and a s (mz) are less well determined 
while rriH, before the LHC, was largely unknown. In this Section we discuss the EW fit 
without the new input on m# from the LHC, in order to compare the limits so derived on 
mil with the LHC data. The discussion of the LHC results will follow in the next Sections. 
Among the light fermions, the quark masses are badly known, but fortunately, for the 
calculation of radiative corrections, they can be replaced by a(mz), the value of the QED 
running coupling at the Z mass scale. The value of the hadronic contribution to the running, 
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embodied in the value of A«^(m|) (see Fig. 49, [213] ) is obtained through dispersion 
relations from the data on e + e~ — > hadrons at moderate centre-of-mass energies. From the 
input parameters one computes the radiative corrections to a sufficient precision to match 
the experimental accuracy. Then one compares the theoretical predictions with the data 
for the numerous observables which have been measured [285], checks the consistency of 
the theory and derives constraints on m t , a s (mz) and m#. 
The basic tree level relations: 

g2 -° F g 2 sm 2 9 w = e 2 = 4ira (327) 



8m 2 w y/2' 



can be combined into 



sm2 ^7lS< (328) 



Always at tree level, a different definition of sin 9 W is from the gauge boson masses: 



m W „ 1 -2 a 1 m W 



, -po = l sin 2 ^ = l--^ (329) 

m z cos z t>w m z 

where po — 1 assuming that there are only Higgs doublets. The last two relations can be 
put into the convenient form 

m 2 w m 2 w = na 
m\' m 2 z ^2G F m 2 z 

Beyond tree level, these relations are modified by radiative corrections: 

vr? w Txa(m z ) 1 

m\ ] m\ y/2G F m 2 z l- &r w 

Apm (331) 



m w 



m 2 z cos 2 9y/ 



The Z and W masses are to be precisely defined, for example, in terms of the pole position in 
the respective propagators. Then, in the first relation the replacement of a with the running 
coupling at the Z mass a(mz) makes Arw completely determined at 1-loop by purely weak 
corrections {Gp is protected from logarithmic running as an indirect consequence of (V-A) 
current conservation in the massless theory). This relation defines Ar^ unambiguously, 
once the meaning of m wz and of a(mz) is specified (for example, MS). On the contrary, 
in the second relation Ap m depends on the definition of sin 2 9w beyond the tree level. For 
LEP physics sin 2 9w is usually defined from the Z — > p + p~ effective vertex. At the tree 
level the vector and axial- vector couplings gy and g\ are given in eqs.(254). Beyond the 
tree level a corrected vertex can be written down in terms of modified effective couplings. 
Then sin 2 9\y = sin 2 9 e ff is in general defined through the muon vertex: 



g v jg\ = l-4sin 2 0, 



eff 



119 



sin O eff = (1 + Ak)s , s c - 
1 



00 V2G F m% 

gf = \(1 + Ap) (332) 

We see that Sg and c 2 , are "improved" Born approximations (by including the running of 
a) for sm 2 9 e ff and cos 2 O e ff- Actually, since in the SM lepton universality is only broken 
by masses and is in agreement with experiment within the present accuracy, in practice 
the muon channel can be replaced with the average over charged leptons. 

We can write a symbolic equation that summarizes the status of what has been com- 
puted up to now for the radiative corrections (we list some recent work on each item from 
where older references can be retrieved) Arw [286], Ap [287] and Ak [288]: 

2 4 fi 

Ar w , Ap, Ak = g\l + a s ) + g 2 ^-(a 2 s + a s s ) + g 4 + g 4 ^a s + g e ^... (333) 

The meaning of this relation is that the one loop terms of order g 2 are completely known, 
together with their first order QCD corrections; the second and third order QCD corrections 
are only known for the g 2 terms enhanced by mf/m^); the two loop terms of order g 4 are 
completely known, while, only for Ap, the terms g 4 a s enhanced by the ratio mj/m^y and 
the terms g 6 ^- are also computed. 

m w 

In the SM the quantities Ar^y, Ap, Ak, for sufficiently large m t , are all dominated 
by quadratic terms in m t of order Gpm 2 . The quantity Ap m is not independent and can 
expressed in terms of them. As new physics can more easily be disentangled if not masked 
by large conventional m t effects, it is convenient to keep Ap while trading Arw and Ak 
for two quantities with no contributions of order Gpm 2 . One thus introduces the following 
linear combinations (epsilon parameters) [289]: 

ei = Ap, 

e 2 = c^Ap + S ° 2 A ^ W 2 - 2s 2 Ak, 
c o s o 

e 3 = c 2 Ap + (c 2 -s 2 )Ak. (334) 

The quantities e 2 and e 3 no longer contain terms of order Gpw? t but only logarithmic terms 
in m t . The leading terms for large Higgs mass, which are logarithmic, are contained in e\ 
and £3. To complete the set of top-enhanced radiative corrections one adds defined from 
the loop corrections to the Zbb vertex. One modifies gy and g\ as follows: 

9 b A = -i(l + ^)(l + <%), 
^ l-4/3sin 2 fl e// + 6 b 

9 b A l + e 6 ' 1 ] 

€b can be measured from = T(Z — > bb)/Y(Z — > hadrons) (see Fig. 49). This is clearly 
not the most general deviation from the SM in the Z — > bb vertex but is the quantity 
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where the large m t corrections are located in the SM. Thus, summarizing, in the SM one 
has the following "large" asymptotic contributions: 



3G F m 2 3G F m l w 2 m H 

ex = =■ — pr-tan t) w ln +...., 

8ir 2 V2 4ir 2 V2 rnz 

G F m 2 w m t 

6 2 = — prill + ...., 

2n 2 V2 rn z 
G F m 2 w m H G F m 2 w m t 

£s — 7= In — — In — , 

rnz Qn 2 V2 m z 
Cb = ,336, 

The 6i parameters vanish in the limit where only tree level SM effects are kept plus pure 
QED and/or QCD corrections. So they describe the effects of quantum corrections (i.e. 
loops) from weak interactions. A similar set of parameters are the S, T, U parameters [290]: 
the shifts induced by new physics on S, T and U are proportional to those induced on e 3 , 
e\ and e 2 , respectively. In principle, with no model dependence, one can measure the four 
€i from the basic observables of LEP physics T(Z — > /i + /i~), A FB and R b on the Z peak 
plus my/- With increasing model dependence, one can include other measurements in the 
fit for the ej. For example, use lepton universality to average the /i with the e and r final 
states, or include all lepton asymmetries and so on. The present experimental values of 
the €i, obtained from a fit of all LEP1-SLD measurements plus mw, are given by [213]: 

ei 10 3 = 5.4 ±1.0, e 2 10 3 = -8.9 ± 1.2, 

e 3 10 3 = 5.34 ±0.94, t b 'lO 3 = -5.0 ± 1.6. (337) 

Note that the e parameters are of order a few in 10 -3 and are known with an accuracy 
in the range 15 — 30%. As discussed in the next Section, these values are in agreement 
with the SM with a light Higgs. All models of new physics must be compared with these 
findings and pass this difficult test. 

3.12 Results of the SM Analysis of Precision Tests 

The electroweak Z pole measurements, combining the results of all the experiments, 
plus the W mass and width and the top mass m t , are summarised in Fig 49, as of March 
2012 [213]. The primary rates are given by the pole cross sections for the various final 
states a ; ratios thereof correspond to ratios of partial decay widths: 

m z L z o l In J- h 

Here Tu is the partial decay width for a pair of massless charged leptons. The partial 
decay width for a given fermion species contains information about the effective vector and 
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Figure 49: Summary of electroweak precision measurements at high Q 2 [213]. The first block shows the 
Z-pole measurements. The second block shows additional results from other experiments: the mass and 
the width of the W boson measured at the Tevatron and at LEP-2, the mass of the top quark measured 
at the Tevatron, and the contribution to a of the hadronic vacuum polarization. The SM fit results are 
also shown with the corresponding pulls (differences data - fits in units of standard deviations). 



axial- vector coupling constants of the neutral weak current: 

T fG F m 3 



ff 



C 6vrf ^ afCM + g2vfC ^ + Aew/QCD ' (339) 



where is the QCD colour factor, C{A,v}f are final-state QCD/QED correction factors 
also absorbing imaginary contributions to the effective coupling constants, g a f and g v f are 
the real parts of the effective couplings, and A contains non-factorisable mixed corrections. 

Besides total cross sections, various types of asymmetries have been measured. The 
results of all asymmetry measurements are quoted in terms of the asymmetry parameter 
Af, defined in terms of the real parts of the effective coupling constants, g a f and g v f, as: 

a _ o 9vf9af _ 9vf/ 9af A 0J _ 3 

9 v f + 9af i - + {9vf/g a f) 4 
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The measurements are: the forward-backward asymmetry (A° BB ), the tau polarization 
(A T ) and its forward backward asymmetry (A e ) measured at LEP, as well as the left-right 
and left-right forward-backward asymmetry measured at SLC (A e and Af, respectively). 
Hence the set of partial width and asymmetry results allows the extraction of the effective 
coupling constants. 



,0,1 
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Figure 50: Summary of sin 2 9 e ff precision measurements at high Q 2 [213]. 

The various asymmetries determine the effective electroweak mixing angle for leptons 
with highest sensitivity (see Fig. 50). The weighted average of these results, including 
small correlations, is: 

sin 2 9 eff = 0.23153 ± 0.00016, (341) 

Note, however, that this average has a x 2 of 11.8 for 5 degrees of freedom, corresponding 
to a probability of a few %. The \ 2 is pushed up by the two most precise measurements 
of sin 2 # e //, namely those derived from the measurements of A\ by SLD, dominated by 
the left-right asymmetry A° LR , and of the forward-backward asymmetry measured in bb 
production at LEP, A F b B , which differ by about 3crs. 

We now expand the discussion on the SM fit of the data. One can think of different types 
of fit, depending on which experimental results are included or which answers one wants 
to obtain. For example, in Table 4 we present in column 1 a fit of all Z pole data plus rriw 
and Vw (this is interesting as it shows the value of m t obtained indirectly from radiative 
corrections, to be compared with the value of m t measured in production experiments), in 
column 2 a fit of all Z pole data plus m t (here it is mw which is indirectly determined), 
and, finally, in column 3 a fit of all the data listed in Fig. 49 (which is the most relevant 
fit for constraining m#). From the fit in column 1 we see that the extracted value of m t 
is in good agreement with the direct measurement (see Fig 49). Similarly we see that the 
experimental measurement of m w is larger by about one standard deviation with respect 
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to the value from the fit in column 2. We have seen that quantum corrections depend 
only logarithmically on m#. In spite of this small sensitivity, the measurements are precise 
enough that one still obtains a quantitative indication of the mass range. From the fit in 
column 3 we obtain: log 10 tuh (GeV) = 1.97 ± 0.12 (or m# = 94^4 GeV). This result 
on the Higgs mass is particularly remarkable. The value of log 10 m#(GeV) is compatible 
with the small window between ~ 2 and ~ 3 which is allowed, on the one side, by the 
direct search limit (m# > 114 GeV from LEP-2 [213]), and, on the other side, by the 
theoretical upper limit on the Higgs mass in the minimal SM, m# < 600 — 800 GeV [291] 
to be discussed in Sect. 3.13. 



Fit 
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m t , mw, IV 


m t (GeV) 
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173.26 ±0.89 


m H (GeV) 
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12211? 


OA+29 


log [m H (GeV)} 


2.17 ±+0.38 




2.09 ±0.17 


1.97 ±0.12 


a s (m z ) 


0.1190 ±0.0025 




0.1191 ±0.0027 


0.1185 ±0.0026 


m w (MeV) 


80381 ± 13 


80363 ± 20 


80377 ± 12 



Table 4: Standard Model fits of electroweak data [213]. All fits use the Z pole results and Aa h ^ d (m 2 z ) as 
listed in Fig. 49. In addition, the measurements listed on top of each column are included in that case. 
The fitted W mass is also shown [213] (the directly measured value is mw = 80385 ± 15 MeV). 

Thus the whole picture of a perturbative theory with a fundamental Higgs is well sup- 
ported by the data on radiative corrections. It is important that there is a clear indication 
for a particularly light Higgs: at 95% c.l. m H < 152 GeV (which becomes m H < 171 GeV 
including the input from the LEP2 direct search result). This was quite encouraging for 
the LHC search for the Higgs particle. More in general, if the Higgs couplings are removed 
from the Lagrangian the resulting theory is non renormalizable. A cutoff A must be in- 
troduced. In the quantum corrections logm# is then replaced by log A plus a constant. 
The precise determination of the associated finite terms would be lost (that is, the value 
of the mass in the denominator in the argument of the logarithm). A heavy Higgs would 
need some unfortunate accident: the finite terms, different in the new theory from those 
of the SM, should by chance compensate for the heavy Higgs in a few key parameters of 
the radiative corrections (mainly e\ and e 3 , see, for example, [289]). Alternatively, addi- 
tional new physics, for example in the form of effective contact terms added to the minimal 
SM lagrangian, should accidentally do the compensation, which again needs some sort of 
conspiracy. 

To the list of precision tests of the SM one should add the results on low energy tests 
obtained from neutrino and antineutrino deep inelastic scattering (NuTeV [292]), parity 
violation in Cs atoms (APV [293]) and the recent measurement of the parity- violating 
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Figure 51: The data for sin 2 0^ are plotted vs m#. The theoretical prediction for the measured value 
of m t is also shown. For presentation purposes the measured points are shown each at the win value that 
would ideally correspond to it, given the central value of mt (updated from [295]) 

asymmetry in Moller scattering [294] . When these experimental results are compared with 
the SM predictions the agreement is good except for the NuTeV result that shows a devi- 
ation by three standard deviations. The NuTeV measurement is quoted as a measurement 
of sin 2 6w = 1 — mym| from the ratio of neutral to charged current deep inelastic cross- 
sections from and using the Fermilab beams. But it has been argued and it is now 
generally accepted that the NuTeV anomaly probably simply arises from an underestima- 
tion of the theoretical uncertainty in the QCD analysis needed to extract sin 2 9yy. In fact, 
the lowest order QCD parton formalism on which the analysis has been based is too crude 
to match the experimental accuracy. 

When confronted with these results, on the whole the SM performs rather well, so 
that it is fair to say that no clear indication for new physics emerges from the data. 
However, as already mentioned, one problem is that the two most precise measurements 
of sin 2 9 e g from A LR and A b FB differ by about 3cr's. In general, there appears to be a 
discrepancy between sin 2 9 e g measured from leptonic asymmetries ((sin 2 # c ff)z) an d from 
hadronic asymmetries ((sin 2 9 e g)h). In fact, the result from Alr is in good agreement 
with the leptonic asymmetries measured at LEP, while all hadronic asymmetries, though 
their errors are large, are better compatible with the result of A b FB . These two results for 
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sin 2 9 e fi are shown in Fig. 51 [295]. Each of them is plotted at the m# value that would 
correspond to it given the central value of m t . Of course, the value for m# indicated by 
each sin 2 # e ff has an horizontal ambiguity determined by the measurement error and the 
width of the ±ler band for m t . Even taking this spread into account it is clear that the 
implications on uih are sizably different. One might imagine that some new physics effect 
could be hidden in the Zbb vertex. Like for the top quark mass there could be other 
non decoupling effects from new heavy states or a mixing of the b quark with some other 
heavy quark. However, it is well known that this discrepancy is not easily explained in 
terms of some new physics effect in the Zbb vertex. A rather large change with respect 
to the SM of the b-quark right handed coupling to the Z is needed in order to reproduce 
the measured discrepancy (precisely a ~ 30% change in the right-handed coupling), an 
effect too large to be a loop effect but which could be produced at the tree level, e.g., by 
mixing of the b quark with a new heavy vectorlike quark [296]), or some mixing of the Z 
with ad hoc heavy states [297]. But then this effect should normally also appear in the 
direct measurement of Aj, performed at SLD using the left-right polarized b asymmetry, 
even within the moderate precision of this result. The measurements of neither At, at SLD 
nor Rb confirm the need of such a large effect (recently a numerical calculation of NLO 
corrections to Rb [298] gave a rather large result that worsens the agreement of theory and 
experiment for this quantity [193], but the corresponding coupling constant shift would still 
be negligible with respect to that implied by A FB ). Alternatively, the observed discrepancy 
could be simply due to a large statistical fluctuation or an unknown experimental problem. 
As a consequence of this problem, the ambiguity in the measured value of sin 2 9 e g is in 
practice larger than the nominal error, reported in Eq. 341, obtained from averaging all 
the existing determinations, and the interpretation of precision tests is less sharp than it 
would otherwise be. 

We have already observed that the experimental value of my/ (with good agreement 
between LEP and the Tevatron) is a bit high compared to the SM prediction (see Fig. 52). 
The value of itlh indicated by my/ is on the low side, just in the same interval as for 
sin 2 measured from leptonic asymmetries. 

In conclusion, the experimental information on the Higgs sector, obtained from EW 
precision tests at LEP-1,2 and the Tevatron can be summarized as follows. First, the rela- 
tion Myj = M§ cos 2 9\y, Eq.(280), modified by small, computable radiative corrections, has 
been experimentally proven. This relation means that the effective Higgs (be it fundamen- 
tal or composite) is indeed a weak isospin doublet. The direct lower limit uih > 114.5 GeV 
(at 95% c.l.) was obtained from searches at LEP-2. The radiative corrections computed 
in the SM when compared to the data on precision EW tests lead to a clear indication for 
a light Higgs, not too far from the present lower bound. The upper limit for itlh in the 
SM from the EW tests depends on the value of the top quark mass m t . The CDF and 
DO combined value after Run II is at present [213] m t = 173.2 ± 0.9 GeV. As a conse- 
quence the limit on m# from the LEP and Tevatron measurements is rather stringent [213]: 
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Figure 52: The data for my/ are plotted vs m t [213]. 

mil < 171 GeV (at 95% c.L, after including the information from the 114.5 GeV direct 
bound). 

3.13 The Search for the SM Higgs 

The Higgs problem is really central in particle physics today. On the one hand, the 
experimental verification of the Standard Model (SM) cannot be considered complete until 
the structure of the Higgs sector is not established by experiment. On the other hand, 
the Higgs is also related to most of the major problems of particle physics, like the flavour 
problem and the hierarchy problem, the latter strongly suggesting the need for new physics 
near the weak scale (that so far was not found) . In turn the discovery of new physics could 
clarify the dark matter identity. It was already clear before the LHC that some sort of 
Higgs mechanism is at work. The W or the Z with longitudinal polarization that we observe 
are not present in an unbroken gauge theory (massless spin-1 particles, like the photon, are 
transversely polarized): the longitudinal degrees of freedom for the W or the Z are borrowed 
from the Higgs sector and are an evidence for it. Also, at LEP it has been precisely 
established that the gauge symmetry is unbroken in the vertices of the theory: all currents 
and charges are indeed symmetric. Yet there is obvious evidence that the symmetry is 
instead badly broken in the masses. Not only the W and the Z have large masses, but the 
large splitting of, for example, the t-b doublet shows that even a global weak SU(2) is not 
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at all respected by the fermion spectrum. This is a clear signal of spontaneous symmetry 
breaking and the implementation of spontaneous symmetry breaking in a gauge theory is 
via the Higgs mechanism. The big questions are about the nature and the properties of the 
Higgs particle(s). The search for the Higgs boson and for possible new physics that could 
accompany it has been the main goal of the LHC from the start. On the Higgs the LHC 
should answer the following questions: do some Higgs particles exist? Which ones: a single 
doublet, more doublets, additional singlets? SM Higgs or SUSY Higgses? Fundamental 
or composite (of fermions, of WW...)? Pseudo-Goldstone boson of an enlarged symmetry? 
A manifestation of large extra dimensions (5th component of a gauge boson, an effect of 
orbifolding or of boundary conditions...)? Or some combination of the above or something 
so far unthought of? By now we have a candidate Higgs boson that really looks like the 
simplest realization of the Higgs mechanism, as described by the minimal SM Higgs. In 
the following we first consider the apriori expectations for the Higgs sector and then the 
profile of the Higgs candidate discovered at the LHC. 



3.14 Theoretical bounds on the SM Higgs Mass 

A strong argument indicating that the solution of the Higgs problem could not be too far 
away (that is, either discovering the Higgs or finding the new physics that complicates the 
picture) is the fact that, in the absence of a Higgs particle or of an alternative mechanism, 
violations of unitarity appear in some scattering amplitudes at energies in the few TeV 
range [299]. In particular, amplitudes involving longitudinal gauge bosons (those most 
directly related to the Higgs sector) are affected. For example, at tree level in the absence 
of Higgs exchange, for s » m 2 z one obtains: 

A(W+W£ ->• Z L Z L ) no mggs ~ i± (342) 

In the SM this unacceptable large energy behaviour is quenched by the Higgs exchange 
diagram contribution: 

A(W+W£ -> Z L Z L ) mggs ~ -i f (343) 
Thus the total result in the SM is: 

2 

S777 

A{W+W£ Z l Zl)sm - ~ V( a -<) (344) 

which at large energies saturates at a constant value. To be compatible with unitarity 
bounds one needs m 2 H < 4it\/2/Gf or mg < 1.5 TeV. This is an important theorem 
that guarantees that either the Higgs boson(s) or new physics or both must be present in 
the few TeV energy range. 

It is well known, as described in [300] and references therein, that in the SM with only 
one Higgs doublet an upper bound on m H (with mild dependence on m t and the QCD 
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coupling a 8 ) is obtained from the requirement that the perturbative description of the 
theory remains valid up to a large energy scale A where the SM model breaks down and 
new physics appears. Similarly a lower limit on m# can be derived from the requirement of 
vacuum stability [301], [302], [303] (or, in milder form, of a moderate instability, compatible 
with the lifetime of the Universe [304,305]). The Higgs mass enters because it fixes the 
initial value of the quartic Higgs coupling A in its running up to the large scale A. We now 
briefly recall the derivation of these limits. 

The upper limit on the Higgs mass in the SM is clearly important for an apriori assess- 
ment of the chances of success for the LHC as an accelerator designed to solve the Higgs 
problem. One way to estimate the upper limit [300] is to require that the Landau pole 
associated with the non asymptotically free behaviour of the A0 4 theory does not occur 
below the scale A. The running of A (A) at one loop is given by: 

d\ 3 

— = — ^-[A 2 + 3Xh 2 — 9/i 4 + small gauge and Yukawa terms] , (345) 
dt Aiii 1 

with the normalization such that at t = 0, A = Ao = m 2 H /2v 2 (from the minimum condition 
in Eq.(288)) and the top Yukawa coupling is given by h® = m t /v. The initial value of A 
at the weak scale increases with m# and the derivative is positive at large A (because of 
the positive A 2 term - the Xip 4 theory is not asymptotically free - which overwhelms the 
negative top- Yukawa term). Thus, if tuh is too large, the point where A computed from the 
perturbative beta function becomes infinite (the Landau pole) occurs at too low an energy. 
Of course in the vicinity of the Landau pole the 2-loop evaluation of the beta function is 
not reliable. Indeed the limit indicates the frontier of the domain where the theory is well 
described by the perturbative expansion. Thus the quantitative evaluation of the limit is 
only indicative, although it has been to some extent supported by simulations of the Higgs 
sector of the EW theory on the lattice. For the upper limit on itlh one finds [300] 

m H < 180 GeV for A ~ M GUT - M Pl 

ttih < 0.5 - 0.8 TeV for A ~ 1 TeV. (346) 

As for a lower limit on the SM Higgs mass, a possible instability of the Higgs potential 
V[4>] is generated by the quantum loop corrections to the classical expression of V[4>}. 
At large the derivative V'[<j)} could become negative and the potential would become 
unbound from below. The one- loop corrections to V[<f>] in the SM are well known and 
change the dominant term at large <fi according to \<p A — > (A + 7 log 2 /A 2 )0 4 . This 
one-loop approximation is not enough in this case, because it fails at large enough 0, 
when 7 log </> 2 /A 2 becomes of order 1. The renormalization group improved version of the 
corrected potential leads to the replacement A0 4 — > A(A)0 /4 (A) where A (A) is the running 
coupling and (p'(^) = 4> ex P f l{t')dt', with j(t) being an anomalous dimension function 
and t = logA/f (v is the vacuum expectation value v = (2\/2GV) _1 / 2 )- As a result, the 
positivity condition for the potential amounts to the requirement that the running coupling 
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A (A) never becomes negative. A more precise calculation, which also takes into account 
the quadratic term in the potential, confirms that the requirements of positive A(A) leads 
to the correct bound down to scales A as low as ~ 1 TeV. We see that, for mu small and 
mt fixed at its measured value, A decreases with t and can become negative. If one requires 
that A remains positive up to A = 10 16 -10 19 GeV, then the resulting bound on mu in the 
SM with only one Higgs doublet, obtained from a recent state of the art calculation [304] 
is given by: 



m H (GeV) > 129.4 + 1.4 



m t (GeV) - 173.1 
07 



0.5 



a s (m z ) - 0.1184 
0.0007 



(347) 



The estimate of the ambiguity associated with m t can be questioned: is the definition of 
mass as measured at the Tevatron relevant for this calculation [306]? Note that this limit is 
evaded in models with more Higgs doublets. In this case the limit applies to some average 
mass but the lightest Higgs particle can well be below, as it is the case in the minimal 
SUSY extension of the SM (MSSM). 

In conclusion, for mt ~ 173 GeV, only a small range of values for mu is allowed, 
130 < rriH < ~ 180 GeV, if the SM holds and the vacuum is absolutely stable up to an 
energy scale A ~ Mqvt or Mpi. For Higgs masses below this range one can still have a 
domain where the SM is viable because the vacuum can be unstable but with a lifetime 
longer than the age of the Universe. We shall come back on that later (see Fig. 55, [304]). 



3.15 SM Higgs Decays 

The total width and the branching ratios for the SM Higgs as function of mu are given 
in Fig. 53 [313]. Since the couplings of the Higgs particle are in proportion to masses, when 
mil increases the Higgs particle becomes strongly coupled. This is reflected in the sharp 
rise of the total width with mu- For mu in the range 114-130 GeV, the width is below 5 
MeV, much less than the widths of the W or the Z which have a comparable mass. The 
dominant channel for such a Higgs is H — > bb. In Born approximation the partial width 
into a fermion pair is given by [313], [314]: 

T(H ff) = N c £^m H m 2 ^ (348) 

where /?/ = (1 — Am 2 /m 2 H ) 1 / 2 . The factor of (3 3 appears because parity requires that the 
fermion pair must be in a p-state of orbital angular momentum for a scalar Higgs (with 
parity P=+l), (this factor would be (3 for a pseudoscalar Higgs boson). We see that the 
width is suppressed by a factor m 2 /m 2 H (the Higgs coupling is proportional to the fermion 
mass) with respect to the natural size Gf^h for the width of a particle of mass mu 
decaying through a diagram with only one weak vertex. 

A glance to the branching ratios shows that the branching ratio into r pairs is larger 
by more than a factor of 2 with respect to the cc channel. This is at first sight surprising 
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Figure 53: Left: The total width of the SM Higgs boson as function of the mass. Right: The branching 
ratios of the SM Higgs boson as function of the mass (solid: fermions, dashed: bosons) [313] 



because the colour factor Nq favours the quark channels and the masses of r's and of D 
mesons are quite similar. This is due to the fact that the QCD corrections replace the 
charm mass at the scale of charm with the charm mass at the scale m#, which is lower 
by about a factor of 2.5. The masses run logarithmically in QCD, similar to the coupling 
constant. The corresponding logs are already present in the 1-loop QCD correction that 
amounts to the replacement m 2 — > m 2 q [l + 2a s /7r(log m 2 q /m 2 H + 3/2)] ~ m 2 {m 2 H ). 

The Higgs width sharply increases as the WW threshold is approached. For decay into 
a real pair of V's, with V = W, Z , one obtains in Born approximation [313], [314]: 



T(H -> VV) = ^^5 v p v (l - 4x + 12x 2 ) 
IGttv^ 



(349) 



where f3y = a/1 — 4x with x = my/m 2 H and 8w = 2, 5z = 1. Much above threshold the 
VV channels are dominant and the total width, given approximately by: 

IV ~ 0.5 TeV(^) 3 (350) 

becomes very large, signalling that the Higgs sector is becoming strongly interacting (recall 
the upper limit on the SM Higgs mass in eq.(346)). The VV dominates over the tt because 
of the (3 threshold factors that disfavour the fermion channel and, at large m#, by the 
cubic versus linear behaviour with itlh of the partial widths for VV versus tt. Below the 
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VV threshold the decays into virtual V particles is important: VV* and V*V*. Note in 
particular the dip of the ZZ branching ratio just below the ZZ threshold: this is due to 
the fact that the W is lighter than the Z and the opening of its threshold depletes all other 
branching ratios. When the ZZ threshold is also passed then the ZZ branching fraction 
comes back to the ratio of approximately 1:2 with the WW channel (just the number of 
degrees of freedom: two hermitian fields for the W, one for the Z). 




Figure 54: Typical one-loop diagrams for Higgs decay into 77, Z7 and, for only the quark loop, to gg. 



The decay channels into 77, Z7 and gg proceed through loop diagrams, with the 
contributions from W (only for 77 and Z7 ) and from fermion loops (for all) (Fig. 54). 
We reproduce here the results for T(H — > 77) and T(H — > gg) [313], [314]: 



m 77) = * 3 7= \AwM + }Z N cQ}MTf)\ 2 (351) 



Gpot rn H 

A f 

^»>=1^'5> (T ' >IJ « 352 » 



where 75 = m^/4m^ and: 



Mr) = 4[r+(r-l)/(r)] 



r 2 

A w (t) = -l[2r 2 + 3r + 3(2r-l)/(r)] (353) 



with: 



f{r) = arcsin \fr for r < 1 



/(r) = -^[log^^=|=-m] 2 forr>l (354) 

For H — > 77 (as well as for H — > Z*y) the W loop is the dominant contribution at small 
and moderate m.#. We recall that the 77 mode is a possible channel for Higgs discovery 
only for tuh near its lower bound (i.e for 114 < itlh < 150 GeV). In this domain of m# we 
have F(H — > 77) ~ 6 — 23 KeV. For example, in the limit win « 2mj, or r — > 0, we have 
^4vf(0) = —7 and Af(0) = 4/3. The two contributions become comparable only for m H ~ 
650 GeV where the two amplitudes, still of opposite sign, nearly cancel. The top loop is 
dominant among fermions (lighter fermions are suppressed by m 2 f/m 2 H modulo logs) and, 
as we have seen, it approaches a constant for large mt. Thus the fermion loop amplitude for 
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the Higgs would be sensitive to effects from very heavy fermions, in particular the H — )■ gg 
effective vertex would be sensitive to all possible very heavy coloured quarks (of course 
there is no W-loop in this case and the top quark gives the dominant contribution in the 
loop). As discussed in the QCD Chapter, the gg — >• H vertex provides one of the main 
production channels for the Higgs boson at hadron colliders (another important channel 
at present is the WH associate production). 

3.16 The Higgs Discovery at the LHC 

On July 4th, 2012 at CERN the ATLAS and CMS Collaborations [12, 13] announced 
the observation of a particle with mass around 126 GeV that, within the present accuracy, 
indeed looks like the SM Higgs boson. This is a great breakthrough that, by itself, already 
makes an adequate return for the LHC investment. One main impact of the Higgs discovery 
is the fact that it represents the only known example in physics of a fundamental, weakly 
coupled, scalar field with VEV. Another extremely important result from the LHC at 7 
and 8 TeV center of mass energy is that no new physics signals have been seen so far. This 
negative result is certainly less exciting than a positive discovery, but it is a crucial new 
input that, if confirmed, will be instrumental in re-directing our perspective of the field. 
In this Section we summarize the relevant data on the Higgs signal as they are known at 
present while the analysis of the data from the 2012 LHC run is still in progress. 

A particle compatible with a Higgs boson of mass ran ~ 126 GeV has been observed 
at the LHC (at present with 7 cr in ATLAS, 7 a in CMS and also, at the Tevatron, about 
3a from CDF and DO combined). From a naive combination the overall confidence level is 
about 10 a. The decay channels where the signal has been observed are 77, ZZ*, WW*, bb 
and r + r~. The ATLAS and CMS combined value for the mass is ran = 125.7 ± 0.4. This 
light Higgs is what one expects from a direct interpretation of EW precision tests: no need 
of a fancy conspiracy with new physics to fake a light Higgs while the real one is heavy 
(this well matches with the fact that no conspirators have been spotted: no new physics). 
As shown in the previous Section the observed value of ran is a bit too low for the SM 
to be valid up to the Planck mass with an absolutely stable vacuum (see eq. 347) but it 
corresponds to a metastable value with a lifetime longer than the age of the universe, so 
that the SM can well be valid up to the Planck mass (if one is ready to accept the immense 
fine tuning that this option implies, as discussed in Sect. 3.17). This is shown in Fig. 
55 where the stability domains as functions of ra t and tuh are shown, as obtained from a 
recent state-of-the-art evaluation of the relevant boundaries [304]. 

In order to be sure that this is the SM Higgs boson one must confirm that the spin- 
parity is + and that the couplings are as predicted by the theory. As for the spin (see, 
for example, [315]), the existence of the H — > 77 mode proves that spin cannot be 1 and 
must be either or 2, in the assumption of an s-wave decay. The bb and t + t~ modes 
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Figure 55: Vacuum stability domains in the SM for the observed values of m t and m# [304]. On the right 
an expanded view of the most relevant domain in the mt-m# plane. The dotted contour-lines show the 
scale A in GeV where the instability sets in, for a s {mz) = 0.1184. 



are compatible with both possibilities. With large enough statistics the spin-parity can 
be determined from the distributions of H — > ZZ* — > 4 leptons, or WW* — > 4 leptons. 
Information can also be obtained from the HZ invariant mass distributions in the associated 
production [315]. The existing data appear to favour a scalar versus a pseudoscalar Higgs. 
We do not expect surprises on the spin-parity assignment because, if different, then all the 
lagrangian vertices would be changed and the profile of the SM Higgs particle would be 
completely altered. The tree level couplings of the Higgs are in proportion to masses and, 
as a consequence, are very hierarchical. The loop effective vertices to 77 and gg, g being 
the gluon, are also completely specified in the SM, where no new heavy states exist and 
contribute in the loop. As a consequence the SM Higgs couplings are predicted to exhibit a 
very special and very pronounced pattern (see Fig. 56) which would be extremely difficult 
to fake by a random particle (only a dilaton, particle coupled to the energy-momentum 
tensor, could come close to simulate a Higgs particle, although in general there would be 
a common proportionality factor in the couplings). The hierarchy of couplings is reflected 
in the branching ratios and the rates of production channels, as seen in Figs. 57. 

Within the present (February '13) limited, accuracy the measured Higgs couplings are 
in reasonable agreement with the sharp predictions of the SM, as can be seen from Fig. 58. 
Great interest was excited by a hint of an enhanced Higgs signal in 77 but the evidence 
is far from being established. All together if it is not the SM Higgs it must be a very 
close relative! Still it would be really astonishing if the H couplings would exactly be 
those of the minimal SM, meaning that no new physics distortion reaches an appreciable 
contribution level. Thus, it becomes a firm priority to establish a roadmap for measuring 
the H couplings as precisely as possible. The planning of new machines beyond the LHC 
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Figure 57: The branching ratios of the SM Higgs boson in the mass range m# = 120 130GeV (left) and 
its production cross sections at the LHC for various cm. energies (right) [307]. 
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has already started. Meanwhile many theory groups are studying the best strategies for 
analysing the already available data. For example, one can introduce a universal factor 
multiplying all H^i/i couplings to fermions, denoted by c, and another factor a multiplying 
the HWW and HZZ vertices. Both a and c are 1 in the SM limit. All existing data on 
production times branching ratios are compared with the a- and c-distorted formulae to 
obtain the best fit values of these parameters [308-310]. Now this fit is performed directly 
by the experimental Collaborations as shown in Fig. 59 [311,312]. In conclusion it really 
appears that the Higgs sector of the minimal SM, with good approximation, is realized in 
nature, while many of us just considered it as a temporary toy model to be soon replaced 
by a more complicated reality. 



3.17 Limitations of the Standard Model 

No signal of new physics has been found neither by direct production of new particles 
at the LHC nor in electroweak precision tests nor in flavour physics. Given the success of 
the SM why are we not satisfied with this theory? Once the Higgs particle has been found, 
why we do not declare that particle physics is closed? The reason is that there are both 
conceptual problems and phenomenological indications for physics beyond the SM. On the 
conceptual side the most obvious problems are that quantum gravity is not included in the 
SM and the hierarchy problem. Among the main phenomenological hints for new physics 
we can list coupling unification, dark matter, neutrino masses (discussed in Sect. (3.7)), 
baryogenesis and the cosmological vacuum energy. At accelerator experiments one possible 
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departure from the SM is the muon anomalous magnetic moment that, as discussed in Sect. 
3.9 shows a deviation by about 3 a, but some words of caution should be kept in mind 
as a large fraction of the error is the theoretical uncertainty, in particular on the hadronic 
contribution to light-light scattering [279]. 

The computed evolution with energy of the effective SM gauge couplings clearly points 
towards the unification of the electro-weak and strong forces (GUT's) at scales of energy 
Mqut ~ 10 15 — 10 16 GeV [266] which are close to the scale of quantum gravity, Mpi ~ 
10 19 GeV. One is led to imagine a unified theory of all interactions also including gravity 
(at present superstrings [316] provide the best attempt at such a theory). Thus GUT's and 
the realm of quantum gravity set a very distant energy horizon that modern particle theory 
cannot ignore. Can the SM without new physics be valid up to such large energies? One can 
imagine that some obvious problems of the SM could be postponed to the more fundamental 
theory at the Planck mass. For example, the explanation of the three generations of 
fermions and the understanding of fermion masses and mixing angles can be postponed. 
But other problems must find their solution in the low energy theory. In particular, the 
structure of the SM could not naturally explain the relative smallness of the weak scale of 
mass, set by the Higgs mechanism at fi ~ 1 /y/G~F ~ 250 GeV with Gf being the Fermi 
coupling constant. This so-called hierarchy problem [317] is due to the instability of the SM 
with respect to quantum corrections. This is related to the presence of fundamental scalar 
fields in the theory with quadratic mass divergences and no protective extra symmetry at 
= 0. For fermion masses, first, the divergences are logarithmic and, second, they are 
forbidden by the SU (2) (^) U{1) gauge symmetry plus the fact that at m = an additional 
symmetry, i.e. chiral symmetry, is restored. Here, when talking of divergences, we are not 



137 



worried of actual infinities. The theory is renormalizable and finite once the dependence 
on the cut off A is absorbed in a redefinition of masses and couplings. Rather the hierarchy 
problem is one of naturalness. We can look at the cut off as a parameterization of our 
ignorance on the new physics that will modify the theory at large energy scales. Then it 
is relevant to look at the dependence of physical quantities on the cut off and to demand 
that no unexplained enormously accurate cancellations arise. 

The hierarchy problem can be put in very practical terms (the "little hierarchy prob- 
lem"): loop corrections to the Higgs mass squared are quadratic in A. The most pressing 
problem is from the top loop. With m\ = m\ are + 5m 2 h the top loop gives 

Sml ltop ~ -^^ 2 A 2 ~ -(0.2A) 2 (355) 

If we demand that the correction does not exceed the light Higgs mass indicated by ex- 
periment, A must be close, A ~ o(l TeV). Similar constraints arise from the quadratic A 
dependence of loops with gauge bosons and scalars, which, however, lead to less pressing 
bounds. So the hierarchy problem demands new physics to be very close (in particular 
the mechanism that quenches the top loop). Actually, this new physics must be rather 
special, because it must be very close, yet its effects are not already clearly visible neither 
in precision electroweak tests (the "LEP Paradox" [318]) nor in flavour changing processes 
and CP violation. Examples of proposed classes of solutions for the hierarchy problem 
are: 1) Supersymmetry [319]. In the limit of exact boson-fermion symmetry the quadratic 
divergences of bosons cancel so that only log divergences remain. However, exact SUSY 
is clearly unrealistic. For approximate SUSY (with soft breaking terms and R-parity con- 
ervation), which is the basis for most practical models, A 2 is essentially replaced by the 
splitting of SUSY multiplets, A 2 ~ m2 susY ~ m ord- m particular, the top loop is quenched 
by partial cancellation with s-top exchange, so the s-top cannot be too heavy. 2) Techni- 
color [320]. The Higgs system is a condensate of new fermions. There is no fundamental 
scalar Higgs sector, hence no quadratic divergences associated to the \j? mass in the scalar 
potential. This mechanism needs a very strong binding force, Ktc ~ 10 3 Aqcd- It is diffi- 
cult to arrange that such nearby strong force is not showing up in precision tests. Hence 
this class of models has been disfavoured by LEP, although some special class of models 
have been devised aposteriori, like walking TC, top-color assisted TC etc [321] (for reviews, 
see, for example, [322]). 3) Extra dimensions [323,324] (for pedagogical introductions, see, 
for example, [325]). The idea is that Mpi appears very large, or equivalently that gravity 
appears very weak, because we are fooled by hidden extra dimensions so that either the 
real gravity scale is reduced down to a lower scale, even possibly down to o(l TeV) or 
the intensity of gravity is red shifted away by an exponential warping factor [324]. This 
possibility is very exciting in itself and it is really remarkable that it is compatible with 
experiment. It provides a very rich framework with many different scenarios. 4) "Lit- 
tle Higgs" models [326]. In these models the Higgs is a pseudo-Goldstone boson and extra 
symmetries allow m h ^ only at two-loop level, so that A can be as large as o(10 TeV) 
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with the Higgs within present bounds (the top loop is quenched by exchange of heavy 
vectorlike new quarks with charge 2/3). 5) The anthropic evasion of the problem. The 
observed value of the cosmological constant A also poses a tremendous, unsolved natural- 
ness problem [327]. Yet the value of A is close to the Weinberg upper bound for galaxy 
formation [328]. Possibly our Universe is just one of infinitely many bubbles (Multiverse) 
continuously created from the vacuum by quantum fluctuations. Different physics takes 
place in different Universes according to the multitude of string theory solutions [329] 
(~ 10 500 ). Perhaps we live in a very unlikely Universe but the only one that allows our 
existence [330], [331]. I find applying the anthropic principle to the SM hierarchy problem 
somewhat excessive. After all one can find plenty of models that easily reduce the fine 
tuning from 10 14 to 10 2 : why make our Universe so terribly unlikely? If to the SM we add, 
say, supersymmetry, does the Universe become less fit for our existence? In the Multiverse 
there should be plenty of less fine tuned Universes where more natural solutions are realized 
and yet are suitable for our living. By comparison the case of the cosmological constant 
is a lot different: the context is not as fully specified as the for the SM (quantum gravity, 
string cosmology, branes in extra dimensions, wormholes through different Universes....). 
It is true that the data impose on us a substantial amount of apparent fine tuning and 
certainly the criterium of naturalness has so far failed, so that we are now lacking a reli- 
able argument on where precisely the new physics threshold is located, but still we remain 
confident that some new physics will appear not too far from the weak scale. While I 
remain skeptical I would like here to sketch one possibility on how the SM can be extended 
in agreement with the anthropic idea. If we ignore completely the hierarchy problem and 
only want to reproduce the most compelling data that demand new physics beyond the 
SM, a possible scenario is the following one. The SM spectrum is completed by the just 
discovered light Higgs and no other new physics is in the LHC range (how sad!) except 
perhaps a Z' , for example a Z' B _ L . In particular there is no SUSY in this model. At the 
GUT scale of Mqvt > 10 16 GeV the unifying group is 5*0(10), broken at an interme- 
diate scale, typically M int ~ 10 10 — 10 12 down to a subgroup like the Pati-Salam group 
SU(4)®SU(2) L <g)SU(2) R or SU(3) (g) 17(1) ®SU(2) L (g) SU(2) R [332]. Note that, in 
general, unification in 577(5) would not work because we need a group of rank larger than 
4 in order to allow for a two step (at least) breaking needed, in the absence of SUSY, to re- 
store coupling unification and to avoid a too fast proton decay. An alternative is to assume 
some ad hoc intermediate threshold to modify the evolution towards unification [333] . The 
Dark Matter problem is one of the strongest evidences for new physics. In this model it 
should be solved by axions [52-54]. It must be said that axions have the problem that their 
mass should be ad hoc fixed in order to reproduce the observed amount of Dark Matter. In 
this respect the WIMP (Weakly Interacting Massive Particle) solution, like the neutralinos 
in SUSY models, is much more attractive. Lepton number violation, Majorana neutrinos 
and the see-saw mechanism give rise to neutrino mass and mixing. Baryogenesis occurs 
through leptogenesis [269]. One should one day observe proton decay and neutrino-less 
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beta decay. None of the alleged indications for new physics at colliders should survive (in 
particular even the claimed muon (g-2) [271] discrepancy should be attributed, if not to an 
experimental problem, to an underestimate of the theoretical errors or, otherwise, to some 
specific addition to the above model [334]). This model is in line with the non observation 
of /i — > e7 at MEG [233] , of the electric dipole moment of the neutron [335] etc. It is a very 
important challenge to experiment to falsify this scenario by establishing a firm evidence 
of new physics at the LHC or at another "low energy" experiment. 
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